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Abstract 



We review some theoretical and statistical aspects of the origin of the large-scale 
structure in the Universe, in view of the two most widely known and accepted sce- 
narios: the inflaton scenario (primordial curvature perturbation ( generated by the 
quantum fluctuations of the light scalar field (p that drives inflation, named the in- 
flaton), and the curvaton scenario (C generated by the quantum fluctuations of a 
weakly coupled light scalar field a that does not drive inflation, named the curva- 
ton). Among the theoretical aspects, we point out the impossibility of having a low 
inflationary energy scale in the simplest curvaton model. A couple of modifications 
to the simplest setup are explored, corresponding to the implementation of a second 
(thermal) inflationary period whose end makes the curvaton field 'heavy', triggering 
either its oscillations or immediate decay. Low scale inflation is then possible to at- 
tain with (the Hubble parameter a few Hubble times after horizon exit) being as 
low as 1 TeV. Among the statistical aspects, we study the bispectrum -B^(fci, k2, k^) 
of C whose normalisation /nl gives information about the level of non-gaussianity 
in the primordial curvature perturbation. In connection with /nl, several conserved 
and/or gauge invariant quantities described as the second-order curvature perturba- 
tion have been given in the literature. Wc review each of these quantities showing 
how to interpret one in terms of the others, and analyze the respective expected /nl 
in both the inflaton and the curvaton scenarios as well as in other less known mod- 
els for the generation of primordial perturbations and/or non-gaussianities. The 5N 
formalism turns out to be a powerful technique to compute /nl in multi-component 
slow-roll inflation, as the knowledge of the evolution of some family of unperturbed 
universes is the only requirement. We present for the first time this formalism and 
apply it to selected examples. 
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Chapter 1 
Introduction 



The Friedmann-Robertson- Walker (FRW) cosmological model (known also as the 
standard or Big-Bang cosmological model) [721 H^^t 11871 1188| l22Uj is the successful 
framework that describes the observed properties of the Universe: homogeneity and 
isotropy at large scales, Hubble expansion, and almost 14 billion years of evolution 
in agreement with globular clusters and radioactive isotopes dating. The additional 
predictions of the cosmic background radiation, confirmed by Penzias and Wilson's 
discovery in 1965 [1^1 1168j . and the relative abundances of light elements [01111113 
IH71 1165| I219| I221j in full agreement with observation, have estabhshed a solid base 
for the study of the Universe throughout its history, turning the old speculative 
cosmology into well established and experimentally supported science ^E^. 

The introduction of a period of exponential expansion (called inflationary) [3 
I7^lll7j . prior to the Big-Bang, brought an elegant solution to the horizon, flatness, 
and unwanted relics problems that were present in the original standard cosmo- 
logical model [HI IZni linn I117| I185j . The horizon problem, or why is the Cosmic 
Microwave Background radiation (CMB) temperature highly isotropic?, was solved 
as the accelerated expansion blows up a region initially in thermal equilibrium to a 
much bigger size, making the horizon at the end of the matter dominated era be still 
inside that region; as a result all regions in the sky appear today at the same back- 
ground temperature. The same accelerated expansion makes the comoving horizon 
shrink so rapidly that our local patch of the universe becomes extremely flat despite 
its actual topology, solving this way the flatness problem, or why is our observable 
Universe almost perfectly flat?. The huge dilution of the abundances of unwanted 
relics (e.g. topological defects: magnetic monopoles, cosmic strings, and domain 
walls) caused by the exponential grow of the size of the Universe during inflation 
and the huge production of entropy by the decay of the scalar field that drives in- 
flation, gave solution to the unwanted relics problem, or where are the topological 
defects (and some other troublesome stuff) predicted by the standard cosmology? ^. 

^Be aware however that a minimum of inflationary expansion (at least 70 e- folds which corre- 
spond to a minimum of 10~^^ seconds of inflation) is required to successfully solve the horizon, 
flatness, and unwanted relics problems, assuming standard evolution. 
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In spite of its success at solving the above mentioned problems, the inflationary 
period became perhaps more important because of its ability to stretch the quantum 
fluctuations of the fields living in the FRW spacetime [121 1^ IHi EHl CIS CHSl 
I21()j . making them classical [H ESI IHH I122[ I124j and almost constant soon after 
horizon exit. They correspond to small inhomogeneities in the energy density and 
are responsible, via gravitational attraction, of the large-scale structure seen today 
in the Universe (see Figs. 11.11 and ll.2|) . If this scenario turned to be correct, the 
energy density inhomogeneities should have left their trace in the CMB released 
at the time of recombination. Indeed, the Cosmic Background Explorer (COBE) 
in 1992 |2()(ij found small anisotropies in the CMB temperature of the order of 1 
part in 10^ (with average temperature Tq = 2.725 ± 0.002 K [2^), on scales of order 
thousands of Megaparsecs. With 30 times better angular resolution and sensitivity 
than COBE, the Wilkinson Microwave Anisotropy Probe (WMAP) |223j confirmed 
this picture in 2003 (see Fig. I1.3|l . measuring in turn the cosmological parameters 
with a 1% order precision |i207ij on scales of order tens of Megaparsecs. The PLANCK 
satellite [171], due to be launched in 2007, will be able to refine these observations 
(see Fig. II. 4j) . With 10 times better angular resolution and sensitivity than WMAP, 
PLANCK promises to determine the temperature anisotropies with a resolution of 
the order of 1 part in 10^, and the cosmological parameters with a 0.1% order 
precision. 

The anisotropies in the CMB temperature^ ST/Tq are directly related to the 
perturbations in the energy density Sp/po at the time of recombination (Sachs- 
Wolfe effect), whose primarily origin is the stretched quantum fluctuations of one or 
several scalar fields (pi that fill the Universe during inflation |112| 1191] ^: 



The perturbations in the energy density at the time of recombination can in turn 
be quantified by the gauge-invariant primordial curvature perturbation ( [112. ll9lj : 



which, on fiat slices*^, can be expressed in terms of only the fluctuations in the fields 
(pi. For instance, in the case of only one scalar field (f present during inflation, 

^From now on, and unless otherwise stated, the perturbation 6y in any quantity y will be 
regarded as first-order in cosmological perturbation theory. Unperturbed quantities will be denoted 
by a subscript unless otherwise stated. 

■^In Eqs. Ijl.tll and (|1.2(l all the quantities are evaluated at time of last scattering, being a the 
global expansion parameter and H = a/a the global Hubble parameter. A dot means derivative 
with respect to the cosmic time. The subscripts k stand for the fourier modes with comoving 
wavenumber k. 

choice of coordinates defines a threading of spacetime into lines of fixed spatial coordinates 
Xi, and a slicing into constant time t hypcrsurfaces. The flat slices are defined such that the 
intrinsic spatial curvature ^^^R vanishes in those hypcrsurfaces. 




(LI) 




(1.2) 
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Figure 1.1: The 2-degree Field Galaxy Redsliift Survey (2dFGRS) map ^A^. The 
2dFGRS obtained spectra for 245591 objects, mainly galaxies. Reliable redshifts 
were obtained for 221414 galaxies. The galaxies cover an area of approximately 
1500 square degrees selected from the extended APM Galaxy Survey in three re- 
gions: a north galactic pole strip, a south galactic pole strip and random fields 
scattered around the south galactic pole strip. The figure shows the map of the 
galaxy distribution produced from the completed survey. The filamentary structure 
can be qualitatively compared with that predicted by the Cold Dark Matter model 
with dark energy (ACDM model) (see Fig. II. 2j) (Courtesy of the Anglo-Australian 
Observatory's 2dFGRS group ^Lj)- 



C dH Ca una CnSl is given by 



C = -i/inf— , (1.3) 



where Hi^f is the global Hubble parameter during inflation. The curvature pertur- 
bation is a convenient quantity to describe the primordial perturbations since it 
is conserved on superhorizon scales {k -C aHi^i), as long as the pressure is a unique 
function of the energy density |127^ 1179^ I222j , and it is well defined even after the 
scalar fields (pi have decayed. We will define C in a rigorous way in Chapter |21 

The scalar fields (pi we have been talking about might or might not have dom- 
inated the energy density during inflation and, therefore, driven the inflationary 
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Redshift z = {t = 13.6 Gyr) 




Redshift z = 18.3 {t = 0.21 Gyr) 



Figure 1.2: These slices show the dark matter distribution in the Universe, based 
on the Millennium Simulation |208j : the largest N-body simulation carried out thus 
far (more than 10^° particles). By zooming in on a massive cluster of galaxies, the 
slices highlight the morphology of the structure on different scales, and the large 
dynamic range of the simulation (10^ per dimension in 3D). The zoom extends 
from scales of a few Gpc down to resolved substructures as small as 100 Mpc. 
The slices through the density field are all 15 Mpc/h thick, with h = 0.71 ± 0.03 
|207j . For each redshift, three panels are displayed. Subsequent panels zoom in 
by a factor of four with respect to the previous ones. This filamentary structure 
can be qualitatively compared with that observed by the 2-degree Field Galaxy 
Redshift Survey (see Fig. II. 1|) (Courtesy of the Virgo Consortium for cosmological 
supercomputer simulations [218J). 




Figure 1.3: CMB temperature anisotropies as seen by the WMAP satellite 26 . The 
oval shape is a projection to display the whole sky. The temperature anisotropies 
are found to be of the order of 1 part in 10^. The background temperature is 
To = 2.725 ± 0.002 K; regions at that temperature are in very light blue. The 
hottest regions (in red) correspond to AT = 200/iK. The coldest regions (in very 
dark blue) correspond to AT = -200;uK (Courtesy of the NASA/WMAP Science 
Team [223 )• 



stage prior to the Big-Bang^. Our adopted definition for the inflaton field in this 
thesis will be the light field that dominates the energy density during inflation^ 
and drives the exponential expansion. This field in most cases parameterises the 
distance along the inflationary trajectories. Until recently, the most widely known 
and accepted scenario for the origin of the density perturbations identified the in- 
flaton with the scalar field whose fluctuations were responsible for the primordial 
density perturbations. This scenario, called the inflaton scenario^ 11121 [TT71ll30j . 
describes very well the properties of C, leading to an almost scale invariant power 

^Although the waterfall field in hybrid inflation dominates the energy density, it does not drive 
inflation. This is not in contradiction with our previous statement since, in the strict sense, 
the waterfall field fluctuations are suppressed as the waterfall mass is much bigger than iJ* (the 
star '*' denoting the global Hubble parameter evaluated a few Hubble times after horizon exit) 

mainniinniiniiinni. 

^By light field we mean a field (f) whose mass is much less than H^,. 

^For simplicity we will just consider single-component inflationary models of the slow-roll variety 
in the present introduction and in the following chapter. The multi-component case, which contains 
one infiaton and one or more 'light non-infiaton fields', will be considered in Chapters El and El 
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Figure 1.4: Simulation of the CMB temperature anisotropies as seen by the 
PLANCK satelhte. PLANCK will provide a map of the CMB field at all angu- 
lar resolutions greater than 10 arcminutes and with a temperature resolution of the 
order of 1 part in 10*^ (ten times better than WMAP) (Courtesy of ESA's PLANCK 
mission jl71j ). 



spectrum 



which is defined by the statistical average 

27r2 

(Ck.Ck.) = ^5'(ki+k2)Pc(fc), (1.5) 

calculated over an ensemble of universes. The amplitude and spectral index 
in Eq. (0|) ^: 

^ - J^* , (1-6) 

y/QSTTmp 



rit; = 2r]^-Qe, (1.7) 
are functions of the slow-roll parameters 



1.^ 



-ninf</2o 



^In this thesis we wih use natural units such that c = fi = ks = ^, and Newton's gravitational 
constant given by SttG = Wp^, with nip — 2.436 x 10^^ GeV being the reduced Planck mass. 
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that characterize the inflationary behaviour and satisfy the slow-roll conditions 
e < 1 and |?7<^| < 1 pTTl El IHDl CHSl- The scale of inflation, given by H^, is 
not completely determined in the inflaton scenario by the already measured ampli- 
tude 1^4^ I |2()7j . since the slow-roll parameter e is unknown (except for the upper 
bound e < 0.01 flHT] ). 

The next statistical significant quantity after V^ik), the bispectrum Bi^{ki, k2, k^) 
defined by the statistical average 

(CkiCk.Cka) = {2ny'/'B^{k,, k2, ks)6^{k, + k, + ks) , (1.10) 

is also well described in the inflaton scenario as its normalisation /nl, deflned by 



6 

B(;{ki, k2, ks) = -^/nl 



27r2 27r2 

— 3-P^(A;i)— -3-P^(A;2) + cyclic permutations 
ki ^2 



'1.11^ 



is suppressed by e and r/^ |142j describing a highly gaussian set of perturbations: 

/nl = ^[2r/^ - 6£ - 2e/(fci, k^, k,)] . (1.12) 

The value of the scale dependent function / in the previous expression lies in the 
range < / < 5/6, being precisely determined by the respective wavevector conflg- 
uration |142j . 

We make a couple of observations [3 EHl llll| I112| I19()j which are valid not only 
for the inflaton scenario but also for the curvaton one discussed below. First, the 
amplitude At of the power spectrum of gravitational waves 



Vrik) ^ Al (-^ ] , (1.13) 



which is deflned by the statistical average 

E(4^k.) = ^^'(ki + ^2)VT{k) , (1.14) 

where hij is the tensor perturbation in the perturbed metric tensor, depends only 
on the scale of inflation (given by H^,): 

At ~ . (1-15) 

vrmp 

Second, its tilt ut is a function only of e: 

nT = -2e. (1.16) 

In view of our previous discussion about V(^{k), we can write a consistency relation 
involving A^, At, and ut, which presents itself as a nice prediction of the inflaton 
scenario: ^ 

rTC = ^ = -8nT. (1.17) 
13 



Its confirmation, as well as a negative detection of non-gaussianity by both WMAP 
[TIM ITHHl 1223 and PLANCK EH satellites, would give strong support to the 
inflaton scenario as the correct framework to understand the origin of the large-scale 
structure in the Universe. We will discuss in more detail the inflaton scenario in 
Chapter 121 

An alternative to the inflaton scenario is when the weakly coupled light field a, 
whose fluctuations are responsible for the primordial density perturbations, does not 
dominate the energy density, and therefore does not drive inflation. This scenario 
is called the curvaton scenario |1H8^ I159j (see also Refs. [HH l \1W\ and 
the field a receives the name curvaton. Introduced in 2002, this scenario describes 
also very well the properties of (, with an almost scale invariant power spectrum 
P^(/c) = Aj{k / aHiai)"''^ whose spectral index written as 



n. = 2r]^-2e, {l.lt 



is function of e and 



being V the scalar potential, and whose amplitude Ac_ is given by if,,, the unper- 
turbed component of the curvaton field during inflation a*, and the fractional global 
curvaton energy density VL^ec = Po-o/Ptotaio just before the curvaton decay: 

^^^ILSldec^ (1.20) 

The scale of inflation in the curvaton scenario, given by H^, does not depend on 
e nor on rj^j, instead it depends on the free parameters cr* and f^dec- A conse- 
quence from the previous expression is that there is no analogous to the consis- 
tency relation [c.f. Eq. ()1.17|) ] for the curvaton scenario. However, there are dis- 
tinctive non-gaussian signatures that can allow us to distinguish this model from 
other scenarios for the origin of the large-scale structure in the Universe (see below) 
)14 ^ IT7 | Il2(i| \1'.M\ llH8j . The curvaton scenario will be studied in detail in Chapter |21 
Perhaps the main motivation for the introduction of the curvaton scenario is 
that it liberates the inflaton field from the generation of the primordial perturba- 
tions 11601 1161j . This is particularly good from the particle physics point of 
view since the intrinsic difficulty at embedding the inflaton scenario in a particle 
physics model is greatly alleviated. Indeed, in the inflaton scenario, the energy 
scale of inflation is likely to be quite high^ {H^, < 10^^ GeV) [123 in order to pro- 
duce the required level of primordial perturbations ^ 5 x 10^^ [2Q2j^°- This 

^The upper bound on can be obtained from Eq. p.6|l and the current bound |2()7| on the 
slow-roll parameter e < 0.01. 

^''Nevertheless there are examples where an e parameter of order ^ 10"^'' is naturally obtained, 
so that the right amount of primordial perturbations is generated for if* ~ 10"^ GeV (see e.g. 

Ref. ^rn\). 
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makes very difficult the identification of tlie infiaton field with one of the scalar 
fields present in the supersymmetry (SUSY) breaking sector or with one of the Min- 
imal Supersymmetric Standard Model (MSSM) fiat directions. This is because the 
characteristic SUSY energy scale, which depends on the symmetry breaking scheme, 
goes typically from ~ lO'^ GeV to ~ 10^ GeV |I^|I2E!. In contrast, in the 
curvaton scenario, the energy scale has to be much lower, satisfying < 10^^ GeV 
|138| 11391 1159] , so that the infiaton field does not contribute to the curvature per- 
turbation [c.f. Eq. (^H)]^^. The last bound can be taken as an anti-smoking gun 
for the curvaton scenario because a positive gravitational wave signal would require 

> 10^2 GeV ISSl EH ESI mH ma mni, mhng out the curvaton as the source of 
primordial density fluctuations. Some people can see this as a bad feature of the 
scenario in question as it sends an unpromising message to all those who are making 
big efforts to detect gravitational waves. However, as described in Ref. |17Uj . the 
curvaton scenario may well be consistent with detectable gravitational waves as far 
as the infiaton is a 'heavy' field (m<^ > H^), violating this way the slow-roll condition 
\rj^\ -C 1 and suppressing the amplitude of the curvature perturbation produced 
by the infiaton itself. 

One of our main concerns in this thesis is the possibility to accommodate low- 
scale inflation in the curvaton scenario. Unfortunately, even when the energy scale 
may in principle be greatly reduced compared to the infiaton scenario, the sim- 
plest curvaton model still requires a quite high inflationary energy scale satisfying 
H^, > 10^ GeV |126j . This of course makes impossible the embedding of the infiaton 
field within the framework of a SUSY particle physics model. In Chapters El and 0] we 
explore two modifications to the curvaton model which can instead allow inflation 
at a low scale [SH I149[ I175[ I189j . In the ffist modification [HIj the end of a second 
(thermal) inflationary stage |1()8| I137j , driven by the rolling of a flaton field x 
coupled to the curvaton, makes the curvaton mass rria- increase suddenly at some 
moment after the end of inflation but before the onset of the curvaton oscillations. 
This proposal can work but not in a completely natural way. Nevertheless, we show 
that inflation with as low as 1 TeV or lower is possible to attain. In the second 
modification jl75^ I189j the increment in rricr at the end of the thermal inflation era is 
so huge that the decay rate overtakes the Hubble parameter and the curvaton field 
decays immediately. The advantage of this second modification is that low scale 
inflation is achieved for more natural values in the relevant parameter space. 

Aside from the previously introduced theoretical aspects of the origin of the 
large-scale structure in the Universe, we also study some of its statistical aspects. 

"'^^An interesting scenario is when both the infiaton and the curvaton fields contribute to C 
j7UlllU5j . This is the case where the curvaton starts oscillating around the minimum of its potential 
when it already contributes significantly to the total energy density ptotai- The upper bound 
iJ* < 10^^ GeV is, in this case, therefore relaxed. In this thesis we will consider only the standard 
curvaton scenario, where the infiaton field does not contribute to C and the curvaton oscillations 
begin when the curvaton energy density pa is still subdominant. 

"'^^Of course, inflation in this case is not of the slow- roll variety. Some possibilities are fast-roll 
| 115| or hilltop inflation |31|. 
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Conversely to the inflaton scenario [c.f. Eq. ()1.12j) ]. ( in the simplest curvaton model 
may present a sizable non-gaussian component if the curvaton does not dominate 
the energy density before decaying [HI [13 I13H I138j . More specifically, according 
to the expression 

which is valid in the curvaton scenario, |/nl| ^ 1 is obtained if r is very small. In 
the previous expression r is defined by 

- - (1.22) 



4pro + 3p, 



and is evaluated just before the curvaton decay (being pr the global radiation en- 
ergy density). This is of extreme importance since the next WMAP data release 
Unaiinnil^, or in its defect the future PLANCK satellite data [EnillZI], will ei- 
ther detect non-gaussianity or put strong constraints on /nl, offering the possibility 
of successfully discriminating among the different inflaton and curvaton models. The 
current constraint on /nl, according to the first-year WMAP data, is |/nl| < 10^ 
|l()2j . The next data release is expected to lower this bound by one order of magni- 
tude or so |l()3j . 

The non-linearity parameter /nl, if independent of position, is closely related to 
the second-order curvature perturbation ^2 defined by 

C(X) = Cl(x) + ic2(x) 

= C.(x)-^/nl(C,^(x)-(0), (1.23) 

where (g is gaussian with {(g) = 0. In connection with this issue we point out 
that several conserved and/or gauge invariant quantities described as the second- 
order curvature perturbation have been given in the literature [HI I14UI I142[ I145j . 
In Chapter El we revisit various scenarios for the generation of second-order non- 
gaussianity in (, employing for the first time a unified notation and focusing on /nl 
|131j . When ( first appears a few Hubble times after horizon exit, |/nl| is much 
less than 1 and is, therefore, negligible. Thereafter ( (and hence /nl) is conserved 
as long as the pressure is a unique function of the energy density (adiabatic pres- 
sure) |127| I179| I222j . Non-adiabatic pressure comes presumably only from the effect 
of fields, other than the one pointing along the inflationary trajectory, which are 
light during inflation (light non-infiaton fields) [SHIIZZ!- Our expectation is that, 
although during single- component inflation /nl is constant, multi-component infla- 
tion might generate |/nl| ~ 1 or bigger. We mention some recent proposals where 
non-gaussianity can be generated during the preheating stage following inflation 
jnni inn in2] , and conjecture that preheating can affect /nl only in atypical scenarios 
where it involves light non- inflaton fields |8| I24 | l25 | l62 | 110(1] . We also study the cur- 
vaton scenario and derive Eq. p. 211) . showing that the simplest model typically gives 
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/nl <^ — 1 or /nl = +5/4. The inhomogeneous reheating scenario [^21 EH EZI (see 
also Refs. (USj 11481 ll5Ul I15H I216| I224j ) . where ( is generated by the inhomogeneities 
in the inflaton decay rate during reheating, is quickly reviewed showing that it can 
give a wide range of values for /nl |216[ I224j . One important conclusion from this 
chapter is that it will be crucial to calculate the precise observational limit on /nl 
using second order theory in case that, unless there is a detection, observation could 



eventually provide a limit |/i 



NL 



< 



I unni. 



A new and interesting proposal is the extension to second order of the 6N for- 
malism |19(ij (see also Refs. |21()| 1211] ). used initially to calculate ( at first order 
and the spectral index in multi-component slow-roll models of inflation [196]^^. In 
this formalism C{t,x) is expressed as the perturbation in the amount of expansion 



iV(t,x) = In 



a(t, x) 
a(tin) 



;i.24) 



from an initially fiat slice at time tin to a final slice of uniform energy density at 
time t: 

Cit,x.) = 6N = N{t,x) - No{t) . (1.25) 

Here 

a{t) 



No{t) = In 



a{ti-_ 



(1.26) 



is the unperturbed amount of expansion, and a(t,x) is the local expansion parame- 
ter. Thus, ( up to second order is given by 



at, x) = 5: Ni{t)5<p, +-y: N,,{t)s<pMj , 

i ij 



;i.27) 



where 



Ni 



ON 
d(f)io 

^2 



;i.28) 
;i.29) 



and the fields (pi, evaluated a few Hubble times after horizon exit, are those relevant 
for the generation of (. In view of Eq. ()1.2Hj) . this presents as a powerful method 
to calculate /nl in any multi-component slow-roll inflationary model for the gen- 
eration of ( |132j . In Chapter ini we give for the first time this formalism, which 
allows us to extract all the stochastic properties of ( if the initial field perturbations 
are gaussian^^. The elegance and power of this method lies in the fact that the 



^^For an extension to multi-component non slow-roll models see Ref. |11U| . 

^^In the case that the initial field perturbations arc non-gaussian, there is an additional contri- 
bution to Eq. which is strongly wavevector dependent |2()()j . This contribution is in any 
case very small ,141, being < 10~^. 
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calculation requires only the knowledge of the evolution of some family of unper- 
turbed universes. The following formula is given for /nl in terms of iV(t,0jg) and 
its derivatives |32 t I132j : 



where is a typical cosmological scale and L is the size of the region within which 
the stochastic properties are specified, so that ln(A;L) ~ 1. We apply the above 
formula to the Kadota and Stewart modular inflation model IHOl IHI] , the curvaton 
scenario jlH8niH9^ ll59j . and the multi-component 'hybrid' inflation model of Enqvist 
and Vaihkonen [6^. The relation of this formula to cosmological perturbation theory 
is also explained. 

The conclusions of this thesis are drawn in Chapter [7| 





(1.30) 
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Chapter 2 



Two mechanisms for the origin of 
the large-scale structure 

2.1 Introduction 

In the standard inflationary scenario j^ lll2t[TT7|ll3Uj a single scalar field, named the 
inflaton, is responsible for the solution of the horizon, flatness, and unwanted relics 
problems, as well as for the origin of the large scale structure seen in the observable 
Universe. This double mission for the inflaton field imposes strong constraints on 
the parameters of the inflationary models, leading to big intrinsic difficulties at 
building successful and realistic models of inflation |112[ ll3Uj . To rescue the well 
motivated inflationary models that fail at generating the required level of primordial 
perturbations jlHj, the inflaton field is left in charge of driving inflation only. The 
other task, the generation of the primordial perturbations, is assigned to a weakly 
coupled light field a different from the inflaton. This is the curvaton scenario |1H8[ 
I139tll59j (see also Refs. [^^11161 fl56j ). where the original curvature perturbation 
associated to and produced by a during inflation, is gradually transformed into the 
total curvature perturbation C ^. The conversion process starts during the radiation 
dominated epoch that follows the reheating stage produced by the inflaton decay^. 

This chapter will give some preliminary definitions and basic facts about the 
inflaton and the curvaton scenarios, such as the first-order perturbations in the 
metric, the precise definition of (, the slow-roll conditions, the characteristics of the 
spectrum of ( in both the inflaton and the curvaton scenarios, and the spectrum 
of gravitational waves and the consequences derived from its possible detection. 

^During inflation and until the start of the radiation dominated epoch just after reheating, C,a 
is actually an isocurvature perturbation. The reason of the name isocurvature is because Co- during 
that time does not contribute at all to the total curvature perturbation C. 

^The cause of the conversion process is the relative redshifting between the radiation and the 
curvaton fluid energy densities. The curvaton at this stage is considered a matter fluid since the 
period of its oscillations around the minimum of its potential is much less than the characteristic 
expansion time scale |213| . This makes the average curvaton pressure essentially zero. 
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Having this information at hand, it will be easier to follow the main discussions of 
this thesis that are exposed in Chapters El to IHl 

2.2 Metric perturbations and the primordial cur- 
vature perturbation ( 

Our observable expanding Universe is well described as homogeneous and isotropic, 
down to scales of order of tens of Megaparsecs. The FRW line element 

ds^ = -de + a\t) ( + rHd^ + sin^ ddA , (2.1) 
y 1 — Kr^ j 

describes such an Universe |72l I186| I187| I188t l22Uj , assuming that the galaxies are 
always at the same comoving spherical coordinates r, ^, and 0, and that the inter- 
galactic space is continually increasing with time due to the expansion parameter 
(scale factor) a{t\ The proper time, once synchronized in all the galaxies, runs at 
the same rate everywhere, so the metric is well described by a global cosmic time t. 
Finally, the topology of the Universe can be classified by the parameter k as closed 
(k > 0), open (k = 0), or flat (k < 0), according to whether the space is finite but 
unbounded, infinite with certain curvature, or infinite but strictly flat a la Minkowski 
^. The inflationary stage however shrinks the comoving horizon so rapidly that the 
present Universe looks so extremely fiat. This makes the parameter k be completely 
irrelevant when going back in time, so that we can safely discard it when study- 
ing the inflationary and early Universe processes. Thus, and switching to cartesian 
coordinates and conformal time t] defined by dr] = dt/a, the FRW line element 
looks 

ds^ = a^{r]) {-drf + 5,jdx'dx^) . (2.2) 

2.2.1 First-order perturbations in the FRW line element: 
classification and number of degrees of freedom 

The homogeneity and isotropy assumptions make a good description of the Universe, 
but they are not in any case perfect conditions. The observed galaxy distribution 
and the temperature fluctuations in the CMB (see Figs. II. 11 - 11. 4|l reveal the impor- 
tance of modifying the FRW line element so that it accounts for the small deviations 
from homogeneity and isotropy required to correctly describe such large-scale struc- 
tures. The introduction of small inhomogeneities in the FRW line element, and the 
truncation at first-order of the perturbed Einstein equations, are well justified since 
the galaxy density contrast 5p/po and the CMB temperature anisotropies ST/Tq 
are observed to be of order 10~^ jSHl I1U4[ I169j . In consequence, considering only 

■^The actual value of k is conventional since it depends on the chosen value for the present 
expansion parameter Otoday 
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scalar degrees of freedom, the most generic first-order perturbed line element reads 
jTTl inm ITiHl ITT)^ ITHHl IT^ ITTIH] 



ds'^ = a^{ri) [-(1 + 2(t)G)d7f + 2 diB drjdx' + ((1 - 2^)5^^ + DijE)dx'dx^\ , (2.3) 

where Dij = didj — (l/3)(5jjV^, and 0g, -B, '?/', and are scalar quantities. 

The peculiar form of the perturbed line element in Eq. ()2.3|) . as well as the right 
number of scalar degrees of freedom in that expression, are justified as follows: 

• The full perturbed metric tensor is generically described by scalar, vector, and 
tensor perturbations. 

• The total number of degrees of freedom for a symmetric tensor in an (n + 1)- 
dimensional spacetime is (n + 2)(n + l)/2. 

• The qqo entry can be written as —(1 + 20g), where 0g corresponds to a purely 
scalar perturbation (1 scalar degree of freedom). 

^ — > 

• The Qqi entries can be parameterised by Bi = diB + vi, where V ■ "U = 0. 
The term d^B corresponds to a purely scalar perturbation (1 scalar degree of 
freedom). The term Vi corresponds to a purely vector perturbation (n — 1 
vector degrees of freedom because of the V • i; = constraint). 

• The Qij entries can be written as (1 — 2ilj)6ij + Uij, where ip is a. purely scalar 
perturbation (1 scalar degree of freedom) that accounts for the trace of Qij, 
and Uij is a symmetric traceless tensor. 

• Uij = Ufj + UYj + Ufj is expressed in terms of a purely scalar perturbation 
(encoded in H^), a purely vector perturbation (encoded in H]^), and a purely 
tensor perturbation (encoded in H^). 

• n^- can be written as DijE because it is a symmetric traceless tensor corre- 
sponding to the purely scalar perturbation E (1 scalar degree of freedom). 

• UYj can be written as {diUj + djIli)/2, where Ilj is a purely vector perturbation 
satisfying V ■ 11 = (n — 1 vector degrees of freedom). 

• Ufj corresponds to a purely tensor perturbation such that diUfj = 0. Con- 
sidering the previous three items, and the fact that the number of degrees of 
freedom coming from a traceless symmetric tensor in an n-dimensional space 
is (n + l)n/2 — 1, the number of tensor degrees of freedom is (n — 2)(n + l)/2. 

• The total number of scalar degrees of freedom (4), vector degrees of freedom 
(2(n— 1)), and tensor degrees of freedom ((n— 2)(n+l)/2), add up to reproduce 
the total number of degrees of freedom of the metric tensor g^j^^, given in the 
second item: (n + 2)(n + l)/2. 
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Two very important facts follow from the first-order perturbed Einstein equations 
fn[ 1^ I163j . First, at first order the vector perturbations are decoupled from 
the scalar perturbations being anyway usually neglected due to their rapid decrease 
with time. Second, the tensor perturbations at first order also decouple from their 
scalar counterparts. These two facts, along with the parameterisation of the metric 
tensor discussed in the above items, justify the use of the Eq. ()2.3|) as the most 
generic first-order metric perturbed line element that describes the energy density 
(scalar) perturbations in the FRW spacetime. 



2.2.2 The curvature perturbation tjj and its non- invar iance 
under infinitesimal coordinate transformations 

In the perturbed line element of Eq. ()2.3p . ip represents the intrinsic spatial curvature 
^^'^R on hypersurfaces of constant conformal time r] 

^^^R = 4 V'^ , (2.4) 

where the operator is the comoving Laplacian operator. For this reason the 
quantity tp is usually referred to as the curvature perturbation. The curvature per- 
turbation ip, as well as the other scalar perturbations (pc, B, and E, are however 
not invariant under a coordinate transformation. Indeed, from the most general 
infinitesimal coordinate transformation^ 

T] ^ 7] + ax^), (2.5) 
X, ^ Xi + diPix'') , (2.6) 

the scalar perturbations (pc, B, tp, and E transform as JT] 

(pG ^ (pG-c-m, (2.7) 

B ^ B + ^-P', (2.8) 

^ ^ i^ + l'^'p + m, (2.9) 

E E-2p, (2.10) 

where a prime means derivation with respect to the conformal time r/, and Ti. is the 
global conformal Hubble parameter defined by 7i = a' /a. The above transformation 
rules are under the physical proviso that the perturbed line element ds'^ in Eq. ()2.3|) 
should remain invariant. 

The parameters ^ and (3 that give account of the infinitesimal coordinate transfor- 
mation in Eqs. ()2.5p and ()2.6|) can be adjusted (fixing the slicing and the threading) 

^The vector component that should appear in the transformation of Xi has been discarded. This 
is because the infinitesimal first-order vector shifts contribute to the transformation rules of the 
vector perturbations only. 
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so that two of the four scalar perturbations in Eqs. ()2.7|) to ()2.10j) vanish. The 
longitudinal gauge, for instance, corresponds to choose B and E equal to zero. In 
this specific gauge the gravitational potential 0g becomes equal to the curvature 
perturbation ip up to first order IHEl I143[ I162[ \W^\ I185j as long as the considered 
fluid is described by a perfect isotropic stress, examples of such a fluid being the 
inflaton and/or the curvaton fields. All of this leads us to say with confidence that 
the curvature perturbation ip really represents the effect of the inhomogeneities in 
the FRW spacetime and it is, therefore, the quantity to study. 



2.2.3 The gauge-invariant curvature perturbation ( 

To parameterise adequately the inhomogeneities in the FRW spacetime, we need a 
quantity invariant under the coordinate (gauge) transformations in Eqs. ()2.5|) and 
()2.fij) . This is not completely possible to do, but we may define a quantity which is 
invariant under transformations in time only. This is done taking into account that, 
for any scalar quantity f{x^) different to 0g, B, ip, and E, the transformation law 
in the associated perturbation Sf{x^) is given by 

5f^5f-fli. (2.11) 

The first gauge invariant quantity that we may define is that which represents the 
curvature perturbation ip in the comoving slices, defined them as the constant time 
hypersurfaces where there is no flux of energy. Considering the inflaton field 93, the 
comoving slices coincide with those where ip is uniform {dipcom = 0) pT) tll27^fl(iH^I222j 
so that the time translation ^com(3^'^) required to go from a generic slice to the 
comoving slice is: 

1 S(p 

Ccom = -t[^^~ ^^com] = — T ■ (2-12) 

Therefore the comoving curvature perturbation ipcom, denoted from now on as 
is written in terms of ip and (p in the generic slice as 

C = -i'- ^^com 

= -^-H^. (2.13) 

The overall minus sign is just a convention, chosen in this thesis to match the agreed 
definition of ( by most of the authors. 

The second gauge invariant quantity that we may define represents the curvature 
perturbation ip in the slices of uniform energy density. Again, and following the 
same steps as before, we have to consider the infinitesimal time translation ^um{x^) 
required to go from a generic slice to the uniform energy density slice where 5puni = 0: 

euni=^[5p-5puni] = ^. (2.14) 

Po Po 
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Thus, the curvature perturbation in the uniform density shce V'uni; denoted from 
now on as —C, is given by 

= -i;-n^-^. (2.15) 

Po 

We have chosen to denote both the comoving and the uniform density shce curvature 
perturbations with the same letter ( because they are equivalent on superhorizon 
scales {k -C a-ffinf) ^I]- The superhorizon scale region is of great importance because 
( is conserved in that region if the adiabatic condition, described below, is satisfied. 
In addition, at superhorizon scales ( becomes truly gauge-invariant because the 
contribution of in Eq. ()2.9|) gets suppressed by the spatial derivatives, so that ip 
does not depend on the changes in the threading^. 

The conservation of ( is guaranteed as long as the pressure P is a unique function 
of the energy density p (the adiabatic condition). The last statement follows from 
the local energy conservation equation in the uniform density slicing at large scales 
[17112221: 

p(t) = -3{H + CMt)+Pit,^)]. (2.16) 

If P satisfies the adiabatic condition then it becomes spatially uniform = (Q), 
and so does C- As a result 

C = 0, (2.17) 

because 

C=(C) = |(C) = 0. (2.18) 

If P does not satisfy the adiabatic condition, i.e. if it has a non-adiabatic component 
Pnad, Eq- ()2.1()|1 implies 

C = ^5Pnad. (2.19) 

The curvaton scenario is one example where the existence of a non-adiabatic pressure 
perturbation makes ( evolve from the negligible curvature perturbation produced by 
the inflaton to the right value observed today. The non-adiabatic pressure pertur- 
bation is, in this case, the result of the presence of two weakly interacting fluids, the 
curvaton matter fluid and the radiation fluid, in the period that follows the inflaton 
decay and reheating. 



^Notice that if (/Sq ^ or Pq ^ the curvature perturbations in Eqs. (|2.13|l and H2.15|l blow 
up. To avoid such a disaster we need very small values for Ti. so as to generate the observed value 
for (. If for some reason ip'g = or p'^ = 0, ( as given in Eq. H2.13() or H2.15|l becomes ill defined 
and we would need to look for a better well defined gauge invariant quantity that represents the 
intrinsic curvature perturbation ip. 
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2.3 Inflation and its effect on the spectrum of 
perturbations of a non- dominating massless 
scalar fleld during a de Sitter stage 

Any period in the history of the Universe during which the expansion is accelerated is 
denominated as inflationary |5tl79 |ITT2|lll7[IT3n] . The inflationary stage prior to the 
Hot Big-Bang has the nice property to stretch the quantum fluctuations of the scalar 
fields living in the FRW spacetime [El IHOl IHH EHl 11121 CESl EID] , so that they 
become classical |31 [7S1 El I122t I124j and almost constant, sourcing the primordial 
density inhomogeneities responsible for the presently observed large-scale structure. 
The amplitude of the spectrum of the classical field perturbations (for light fields) 
is generically the same for all kinds of quasi de Sitter models. However, the spectral 
index is written down in a certain way for fields that dominate the energy density 
and in another way for fields that do not. In this section we will describe inflation, 
paying special attention to the constraint imposed by it on the slow-roll parameter 
e, and review the properties of the power spectrum of perturbations of a massless 
scalar field during a de Sitter stage, characterized the latter by a constant Hubble 
parameter ifjnf = -ff* during inflation. 



2.3.1 Inflation 

Inflation can be rigorously defined as the period when the global expansion param- 
eter a satisfies the condition 

a>0. (2.20) 

Such a condition translates into a definite requirement on the global Hubble param- 
eter during inflation H^^f. From the definition of i^inf in terms of a, which can be 
written alternatively as the following evolution equation: 



a(t) = ainicxp ^* H,^dt^ , (2.21) 



being Oini the expansion parameter at some time tmi, the inflationary condition in 
Eq. ^TM is satisfied while 

H,^,>-Hl,. (2.22) 

The last expression reduces to an upper bound on the slow-roll parameter e defined 
in Eq. (fn?|) : 

..-(|l)_<l. ,2.23) 

Inflation is held as long as e < 1, being one possibility when H^^f is constant in 
time. This scenario is commonly recognized as the de Sitter stage, and it will be 
useful at studying the spectrum of perturbations of a massless scalar field. Other 
possibilities correspond to nonvanishing values for -ffinf, being the quasi de Sitter 
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case {e <^ 1) the most popular. Indeed, the quasi de Sitter stage supplemented by 
the condition |?7(^| <^ 1, being rj^p the slow-roll parameter defined in Eq. ()1.9|) . is 
what is known as slow-roll inflation. This variety of infiation will be discussed in 
Subsection 12.5.11 Meanwhile in the following subsection, we will study the fluctua- 
tions of a non- dominating massless scalar field during a de Sitter stage. Later on we 
will generalise these results to the fiuctuations of non-dominating and dominating 
massive scalar fields in a quasi de Sitter stage. 



2.3.2 Spectrum of perturbations of a non- dominating mass- 
less scalar field during a de Sitter stage 

In this subsection we will consider the effects of a de Sitter inflationary stage on 
the fluctuations of a massless scalar fleld that does not dominate the energy density. 
Let's call this fleld (p. During inflation, and before horizon exit, the fluctuations in 
(j) can still be regarded as quantum operators. If we further assume that is almost 
a non-interacting fleld, we can write down the fleld perturbation operator (50(x, t) 
in terms of the usual creation and annihilation operators dj^ and dk: 

50(x, t)= J j^-^ expitk ■ x)<50k(t) , (2.24) 

where 

6Mt) = uJkitW + ^lit)a\ . (2.25) 

As it was discussed in the introduction of this thesis, the properties of the primor- 
dial curvature perturbation ( are specifled by the spectrum Vc_{k), which is deflned 
by the statistical average (CkiCk2) over an ensemble of universes. The same deflnition 
can be applied for the perturbations in 0, so that the only information we need to 
know to calculate Vs^{k) is the quantum state of the Universe during inflation, being 
the most reasonable choice the vacuum state®. Since the universes in the ensemble 
are all in the vacuum state during inflation, the statistical average (50ki'^0k2) is now 
very easy to calculate corresponding to the expectation value (O|(50ki'^0k2 10)- It is 
straightforward then to recognize that the spectrum Vs,j,{k) of the S(j) perturbations, 
deflned by 

(5</'k,50k2) = -^^'(ki + ^2)Vs^ik) , (2.26) 
is given by the simple formula 

Vs^ik) = ^Icukl^ . (2.27) 

To calculate \(jJk\ we need to solve the Klein-Gordon equation for the fluctuations 
in (j). The usual Klein-Gordon equation 

^2 

Uk + 3H,iUk + = , (2.28) 



^The vacuum state guarantees the homogeneity and isotropy in the whole space. 
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is only applicable if the inflaton field (that which dominates the energy density) is 
unable to generate significant primordial perturbations. That guarantees that the 
line element in Eq. ()2.2|) is not modified by any scalar perturbation as is assumed 
not to dominate the energy density during infiation. As a consequence, the usual 
structure of the Klein-Gordon equation is unmodified. 

Eq. ()2.28|) is easily solved by making the following change of variables: 

uJk = —, (2.29) 
a 

and going to conformal time where the expansion parameter in a de Sitter stage 
(i^inf being constant) 

a{t) = Oini exp[H^{t - tini)] , (2.30) 

is given by 

a{v) = -^, (2.31) 

with conformal time rj taking negative values. Thus, the Klein-Gordon equation in 
Eq. reduces to 



2 



+ ( ^ - 1 = ' (2-32) 



^^^expH^ (2.33) 




whose exact solution is 



2k 

The reader might worry about the fact that Eq. ()2.33|) is valid up to a multi- 
plicative constant. To solve this ambiguity we note that in the sub horizon regime 
{k 3> ct-ffinf) Eq. ()2.32|1 reproduces the Klein-Gordon equation for a massless scalar 
field in Minkowski spacetime. The solution for in Eq. ()2.33|1 satisfies the Bunch 
and Davies normalisation EZj on subhorizon scales, so that no integration con- 
stant should amplify it. 

The subhorizon regime is useful in the sense that we can find the adequate 
solution and normalisation for Afc. However, the superhorizon regime [k -C aHiai) is 
much more interesting because C defined in that regime is truly gauge invariant and 
conserved as long as the pressure satisfies the adiabatic condition. From Eqs. ()2.29j] . 
()2.31|) . and ()2.33|) . the magnitude of the mode function for superhorizon scales 
turns out to be constant and it is given by 

|u;fc|^^=. (2.34) 



The latter expression allows us to obtain the spectrum of perturbations in 0, given 
byEq. (jTTfjl : 



V^{k)^[^] . (2.35) 
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As we can see, the spectrum Vs^{k) of the perturbations of a massless scalar field 
in a de Sitter stage is constant and scale invariant. We will see later on that the 
previous amplitude is generically the same for any kind of quasi de Sitter model 
as far as we deal with light fields. A possible scale dependence will arise if the 
field under consideration has a finite mass and/or if the inflationary stage is not 
completely de Sitter (i^inf 7^ 0). 

2.4 The curvaton scenario: an example of a non- 
dominating light scalar field during a quasi de 
Sitter stage 

In the previous section we described the spectrum of perturbations of a non-domina- 
ting massless scalar field during an inflationary period with ifjnf being constant (de 
Sitter stage). Now we move a step ahead considering a non-dominating scalar field 
cr whose mass mg- satisfies m„ <^ (light field), during an inflationary period 
where the Hubble parameter H[^f is not constant but evolves slowly satisfying the 
condition —Hi^f/H^^^ <^ 1 so that e <C 1 (quasi de Sitter stage). We will see that the 
introduction of a slowly varying ifinf and a small mass rricr for the field a results in 
a small scale dependence for Vsa{k) which is parameterised by e and rj„ = m^/3if^. 
The amplitude of Vsa{k) will turn out to be the same as that for a non-dominating 
massless scalar field in a de Sitter stage, which was the case described in Subsection 
12.3.21 This example will help us to understand the properties of the curvature 
perturbation produced in one of the most satisfying models proposed to explain the 
origin of the large-scale structure in the Universe: the curvaton scenario. 

2.4.1 Spectrum of perturbations of a non- dominating light 
scalar field during a quasi de Sitter stage 

The calculation of the spectrum of perturbations of a non- dominating light scalar 
field a during a quasi de Sitter stage closely resembles that done in Subsection 12.3.21 
We will now have to consider the mass mg- of a and the running of ifinf given by 
the slow-roll parameter e. Since the inflaton field is supposed to produce negligible 
curvature perturbation, and a does not dominate the energy density during inflation, 
the line element is still given by Eq. ()2.2j) . This means that the Klein-Gordon 
equation for the mode functions Uk does not change compared with that for the 
mode functions of the background component of a: 



In the quasi de Sitter stage the expansion is almost exponential; nevertheless it 




(2.36) 
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is better described by the following evolution equation: 

a{t) = ai„i[l + H,^{{ti^i)e{t - ti^i)]'/" 



(2.37) 



where -?/inf(iini) is the Hubble parameter at the time tmi- The last expression can be 
written down using the conformal time 7] in an easier way: 



a{7]) 



(2.38) 



having in mind that 7] takes negative values. Going to conformal time and making 
the change of variables 

uJk = —, (2.39) 



the following equation of motion for Afc is obtained: 



K + 



1 



1 



A. 



(2.40) 



where 



1/2 



(2.41) 



The solution for such an equation is given in terms of the Hankel's functions of the 
first and second kind i/W and H^f ^351 IT^ ITTIH] : 



-v[c^ik)Hi'Ji-kr^) + C2ik)Hi'J{-kv)] 



(2.42) 



where ci{k) and C2{k) are integration constants that are determined going to the 
subhorizon regime [k ^ aifjnf), which corresponds to —kr] ^ 1, and normalising the 
solution according to Bunch and Davies [^niEZl^- Indeed, taking into account that in 
the subhorizon regime the Hankel's functions are well approximated by [S^ll85||l95j 



Hil\-kr^:>l] 
Hi'Ji-kr^:^!] 



Txkrj 



Tikr] 



exp 



exp 



7r 



-i(k7] + -Va 



i[kri + -Va 



37r 

+ T 

T 



the Bunch and Davies normalisation 



A. 



exp{—ikr]) 



2k 



(2.43) 
(2.44) 

(2.45) 



^The requirement of a solution for normalised a la Bunch and Davis is justified because 
Eq. I|2.40|l on subhorizon scales reduces to the Klein-Gordon equation for a massless scalar field in 
Minkowski spacetime. 
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is obtained in the subhorizon regime by choosing the following values for the inte- 
gration constants: 



TT 



TT 



C2{k) = 0. 

Thus, the exact solution in Eq. ()2.42j) is rewritten as 



exp 







tt" 












2. 



(2.46) 
(2.47) 



(2.48) 



To find out the spectrum of 5a on superhorizon scales {k <^ aHi^f), corresponding 
to —ki] <^ 1, we need to consider the behaviour of in that regime. This is 
given by the following approximation for the Hankel's function of the first kind on 
superhorizon scales 1185^ |195j 



Hi'Ji-kr^^l] 



exp 



r(3/2) 



(2.49) 



The magnitude of the mode function on superhorizon scales is then almost con- 
stant and approximately given by 



[2{1 - e)r-Hl - e] 



inf 



r(3/2) V«^inf. 



ViT-e 



(2.50) 



where we have used the relation involving v^, e, and established in Eq. ()2.41|1 . 

Making use of the expression in Eq. ()2.27|1 . the spectrum of perturbations in a 
is finally written down as 



V2 



TT 



inf , 



whose spectral index n^^ is given by 

nsa = 21]^ - 2e . 



(2.51) 



(2.52) 



Comparing the previous result with that found for the case of a non-dominating 
massless scalar field during a de Sitter stage [c.f. Eq. ()2.35p ]. we see that the 
amplitude of Vsa{k) is exactly the same for both cases, but now we have a small 
scale dependence parameterised by the slow evolution of i^inf and the mass of 
the light field a, characterized by the smallness of the parameters e and rjcr. 
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2.4.2 The curvaton scenario 



In the basic curvaton setup the Hubble parameter during inflation is assumed 
to be slowly varying (e -C 1), the curvaton energy density pa during inflation is 
assumed to be negligible, and the inflaton fleld is supposed to produce a negligible 
curvature perturbation (r !138| I139| I159j which is imprinted to the radiation fluid 
while the inflaton decays. After the end of inflation the Universe is composed by the 
almost unperturbed radiation fluid that originates from the reheating process fol- 
lowing the inflaton decay, and the weakly coupled llZj light curvaton fleld a ® whose 
unperturbed component is kept frozen at a value ex* until the Hubble parameter 
H becomes of the order of the curvaton mass rrio- during the radiation dominated 
epoch. Once a is unfrozen, it begins oscillating around the minimum of its potential, 
which is taken to be quadratic, with an oscillation period which rapidly becomes 
much less than the characteristic expansion time scale. This ensures that the aver- 
age curvaton pressure vanishes and, therefore, a may be considered as a matter fluid 
|213j . During the oscillatory period, the curvaton energy density decreases with 
time according to p„ oc a~^, while the radiation energy density pr decreases with 
time faster than p„ according to pr oc a~^. Eventually the curvaton will decay, but 
by that time the contribution of p^ to the total energy density will be big enough 
for the original isocurvature perturbation (cr, generated by a during inflation and 
which is not negligible, to become the total curvature perturbation C. 

Since the oscillations of a around the minimum of its potential are so fast, we 
can approximate its energy density by 



where (7a(t,x) is the amplitude of the oscillations. Notice that, under these circum- 
stances, the expression for the curvaton energy density in Eq. ()2.53|) corresponds 
also to the expression for the curvaton potential V{a). The no appearance of quar- 
tic or higher order terms in the potential is essential for the success of the model 
because, otherwise, the density ratio Pa/Pr would not increase with time jTfl . 

Making use of the curvature perturbation deflnition in Eq. ()2.15|) . we can write 
down the total ( in the curvaton scenario as: 



where the total energy density ptotai is simply the addition of the curvaton and radi- 
ation energy densities Po- and pr, that deflne the conserved curvaton and radiation 

*The curvaton field must be weakly coupled to avoid premature thermalisation, which would 
erase all the information about the generatated curvature perturbation |T7j . 




(2.53) 




(2.54) 
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curvature perturbations Co- and C 



9. 



. + (2.55) 

Pao -J Pa-Q 

C . -^-H^J^ = -i, + \^->^. (2.56) 

Pro ^ Pro 

In the above expressions we have employed the background continuity equation 

[zn cnii EE! 

po + 3Hipo + Po) = , (2.57) 

where Pq = for a matter fluid, and Pq = Po/3 for a radiation fluid. Combin- 
ing Eqs. ()2.54|) . ()2.55p . and ()2.56|) . the total curvature perturbation ( can then be 
written down as the weighted sum 

C = (l-r)C. + rC, (2.58) 

with modulation factor 

r = -^ _ 

4pro + ^Pao 

Notice that just at the beginning of the radiation dominated epoch that follows 
the reheating stage produced by the inflaton decay, r is almost zero since p^ is 
negligible by that time; therefore ( ^ (r which is negligible. However, r grows 
in time due to the relative redshifting between po- and pr until when eventually a 
decays. In view of Eqs. ()2.58|1 and ()2.59|1 . the total ( grows then in time approaching 
more and more to the curvaton curvature perturbation One extreme example 
is when a has dominated the energy density before decaying; in that case r ^ 1 
and therefore ( ^ (f^. When a decays, ( is imprinted in remaining radiation fluid 
starting this way the gravitational instability process that ends up with the presently 
observed large-scale structure. 

As we have already pointed out, one of the requirements of the curvaton scenario 
is that the curvature perturbation produced by the inflaton during inflation (r is 
completely negligible compared with that produced by the curvaton during the 
same period: Cr Co-. Under that assumption, the expression for the total C after 
a decays comes from Eq. ()2.58j) as 

C ~ <<x , (2.60) 



in the sudden decay approximation |139j . For a model that goes beyond this ap- 
proximation the expression for the total C in terms of Ca is only obtained by means 
of numerical calculations, the result being in that case |138[ IH 



C ~ ^decC , (2.61) 



^The curvature perturbations Co- and Cr, associated to the curvaton and the radiation fluids 
respectively, are conserved since the adiabatic condition is satisfied separately being both fluids 
non- interacting. 
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where fidec is the fractional global curvaton energy density just before the curvaton 
decay: 

f^dcc = f . (2.62) 

\ Ptotalo / dec 

As we mentioned in Chapter ^ and will explain in Chapter El the normalisation /nl 
of the bispectrum _B^(/ci, /c2, ^3) of C in the curvaton scenario is directly related to 
n^ec if the latter is not so close to 1 [HI El rTTmrHIK] : 

The parameter /nl gives information about the level of non-gaussianity present in 
(, and the actual bound on it, coming from WMAP data |in2j . is |/nl| ^ 10^. This 
bound translates into a lower bound for Qdec, that combined with the obvious energy 
density condition fl^ec < 1, gives the allowed range 

0.01 < fidec < 1 . (2.64) 

The present lower bound on Q^ec is likely to be increased |il03j by the next WMAP 
data release or the future PLANCK satellite data if non-gaussianity effects are not 
detected. 

Once we have studied how the curvature perturbation is produced in the curvaton 
scenario, we proceed now to study the spectrum V(^{k) of (. In view of Eq. ()2.53j) . 
and having in mind that the equation of motion for Saa is the same as that for the 
background field aa^, throughout inflation and during the post-inflationary period, 
as long as the non-gauge invariant curvature perturbation ip is fixed to be zero^°, we 
can relate the contrast in the energy density of a at any time t with the contrast in 
a some time after horizon exit but before the onset of the curvaton oscillations: 

^^2f^U2^. (2.65) 

From Eqs. ()2.55p . ()2.61|) . and ()2.65p . ( is expressed in terms of the perturbations in 
a a few Hubble times after horizon exit: 

C ~ l^dcc— , (2.66) 

and in consequence the spectrum V({k) is given by 



n(^)-i^dec^. (2.67) 



"'^"That is indeed the case while the curvature perturbation in the radiation fluid is taken to 
be neghgible j2t)4| . which is one of the key assumptions in the curvaton scenario. 
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The curvaton field a is a light field whose energy density is negligible during infiation; 
therefore the discussion and results of Subsection 12.4.11 apply [c.f. Eqs. ()2.5ip and 
()2.52|l ]. giving as a result 



The spectral index nc_ is in good agreement with observation, which requires an 
almost-scale invariant power spectrum j2()7j : 

- 0.048 < rif < 0.016. (2.69) 

Unfortunately the Hubble parameter a few Hubble times after horizon exit H^, 
which gives information about the infiationary energy scale, is not predicted by the 
amplitude of V^lk) since a* is an unknown and unbounded parameter. Nevertheless, 
a lower bound for will be found in Chapter IHl by taking into consideration other 
effects that give a different relation between H^, and a^,, and the amplitude A,^ in 
Eq. (ITHHIl once the WMAP normalisation (|v4^| ^ 5 x 10"^) ^U7\ is taken into 
account: 

a, ^ (I.Stt X IQ-^y^ndccH^. (2.70) 

What is important however to emphasise at this point is that the biggest possible 
value for in the curvaton scenario is for sure 10^^ GeV. Otherwise the curvature 
perturbation (r, produced by the inflaton field during infiation, would contribute 
significantly to the total (, spoiling the main motivation for the proposal of the 
curvaton scenario^^. The justification of this assertion will be given in the following 
section. 



2.5 The inflaton scenario 

In this section we will discuss the main facts about the infiaton scenario where 
infiation is assumed to be of the slow-roll variety [Sl lll2l[TT71ll30j . Slow- roll infiation 
corresponds to the case where the infiaton field cp slowly-roll down towards the 
minimum of its potential. We will specify the slow-roll conditions and see what 
their consequences are on the shape of the infiaton potential as well as on the value 
and structural form of the slow- roll parameters e and rj^. Being ip in the infiaton 
scenario the responsible of driving infiation and also of generating the curvature 
perturbation (, the power spectrum V^{k) of ( presents definite signatures that are 
expressed in terms of e and r]^. We will calculate (fc) and see what the constraints 
on the infiaton potential are in order to produce enough primordial perturbations. 
The Hubble parameter during infiation H^, will turn out to be likely quite high 

^^The only way to have iJ* > 10^^ GeV in the curvaton scenario while making (r negligible is 
by requiring the inflaton field not to be light during inflation {m^ > i?*) |17()j . A non slow-roll 
inflationary model is in that case compulsory. 
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(if* < 10^"^ GeV) for the inflaton scenario to be consistent with the amphtude of 
perturbations observed by WMAP. The scale of inflation is, therefore, likely high 
enough to impose severe constraints on concrete inflation models |I67] . 



2.5.1 The slow- roll conditions 

We begin by considering the Friedmann and continuity equations, derived from the 
background Einstein equations for the FRW cosmological model jTTJ IIUH I112j , that 
relate the Hubble parameter at any time with the global energy density and pressure 
of the fluid that fills the Universe: 

= (2.71) 
3mp 

Po = -3ii(po + i^o) • (2.72) 

A direct consequence of both equations is that the second derivative of the global 
expansion parameter a with respect to the cosmic time is given by a simple relation 
involving po and Pq\ 

^ = _^i±|^ . (2.73) 

This expression tells us that, to satisfy the inflationary condition d > 0, the pressure 
of the fluid that fills the Universe must be negative satisfying 

Po + 3Po < . (2.74) 

As an application of the above formula we may study the dynamics of the inflaton 
field (fi knowing that, from the energy momentum tensor for a homogeneous scalar 
field |..112j . the unperturbed energy density and pressure of ip are given by 

P^o = l^l + Vi^o), (2.75) 

P^o = ^00-n^o). (2.76) 
The inflationary condition in Eq. ()2.74j) is then satisfied provided that 

00 < ^(^o) . (2.77) 

The most popular type of inflationary models assume that the kinetic energy 
density of the inflaton field is much less than the potential energy density: 

« Vivo) , (2.78) 

which corresponds intuitively to a very flat potential along which the inflaton field 
if slowly-roU down during inflation towards the minimum of its potential |112| 1130] . 
If the expression in Eq. ()2.78j) is supplemented by the condition 

|<^o| < \mnm\ , (2.79) 
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the inflaton field satisfies what is known as the slow-roll conditions As we 

will see, these conditions can be expressed in terms of the parameters e and rj^ that 
parameterise the spectral index and amplitude of V({k) in the inflaton scenario. 
Notice that, under the slow-roll conditions, the background field (po follows the 
slow-roll equation of motion 

dV 

3//infy3o ^ , (2.80) 

which corresponds to the background Klein-Gordon equation under the condition 
given by Eq. (jTYHl) . 

As discussed in Subsection 12.3. H the requirement to have a period of accelerated 
expansion is easily expressed as an upper bound on the slow-roll parameter e that 
describes the rate of change of the Hubble parameter a few Hubble times after 
horizon exit: 

..-(^)^<l. (2.81) 

The true reason why e is called a slow-roll parameter is because it is constrained 
to be much less than 1 under the slow-roll conditions in Eqs. ()2.78|1 and ()2.79|1 . 
being easily expressed in terms of the unperturbed inflaton potential ^(v^o) and its 
derivative with respect to ip 

The flatness condition on the potential V{(p) required for ip to slowly-roU during 
inflation is evident from the above expression. 

Two slow-roll conditions (Eqs. ()2.78|) and ()2.79|l ) require constraints on two 
slow-roll parameters. One of them is that given in Eq. ()2.82|) in terms of e; the 
other one is given in terms of the parameter tj^, already defined in Eq. 

= ^ - • (2.83) 

The respective constraint on t]^ and its relation with V{(p) are obtained once we 
take into consideration the slow-roll conditions in Eqs. ()2.78|) and ()2.79|) [lllj : 



ml, 



V d^i 



< 1 . (2.84) 



This relation again shows how fiat the potential of the inflaton field ought to be to 
drive inflation. This is particularly good in order to generate enough inflation as 
p} spends a lot of time rolling along the flat part of its potential, which is in turn 
perhaps the main motivation to have an inflationary slow-roll model. 
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From the practical point of view, inflation is said to start when V{ip) satisfies 
both Eqs. ()2.82|) and ()2.84|) . and ends when any of them is violated. Let's however 
remember that, in any case, the slow-roll conditions are sufficient but not necessary 
to drive inflation. Strictly speaking inflation may end some time after the slow-roll 
conditions are violated, but this time is very small compared with the 70 e-folds or 
so of inflation required to solve the horizon, flatness, and unwanted relics problems, 
under standard evolution. 



2.5.2 Spectrum of perturbations of a dominating light scalar 
field during a quasi de Sitter stage 

When we consider a scalar field (p that dominates the energy density during inflation 
and whose perturbations are sizable enough, the spacetime stops being perfectly 
smooth so that we have to leave the comfortable unperturbed metric line element 
in Eq. ()2.2p to adopt the perturbed line element described in Eq. fj2.3p . As a result 
the Klein-Gordon equation for the Fourier modes Uk of the perturbations in ip is 
modified to take into account the backreaction of the metric \W2\ Il(i3| \lHf)\ : 

ujk + SHiniuJk + -T + TTT = -20Gfc7^ h 00^00 + ^ipk^o " k'^Bk^po . (2.85) 

The above equation looks quite difficult to manage but, fortunately, we can eliminate 
some of the scalar degrees of freedom by fixing the gauge and using the perturbed 
Einstein equations for the inflaton field (p \W2\ Il(i3| 118,^^ . For instance, going to 
the longitudinal gauge, we can fix the scalar perturbations B and E to be zero in the 
metric line element of Eq. ()2.3|) . whereas the non-diagonal part of the ij component 
of the perturbed Einstein equations requires 0g = ^ being the stress associated to 
ifo completely isotropic. The modified Klein-Gordon equation reduces in this case 
to 

Ldk + ^HiniiUk + [^ + TTT = -SV'fcT^ h 4?/'fc0o • (2.86) 

To solve the previous equation we still require to know the behaviour of ip. To 
that aim we take advantage of the 00, Oi, and the diagonal part of the ij components 
of the perturbed Einstein equations in the longitudinal gauge I162t |l63t I185j : 

V'^ib 1 ( dv \ 

- 3iJinf (z^ + i/infV') + ^ = ^ ( <^o^0 - 00^ + ^^^A ' (2-87) 

ip + Hi^fij = -^<po6^, (2.88) 
2^ + ij2f)v> + 4ifi„f^ + V^ = -1^ - - ^5<^^ , (2.89) 
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which combined give the following equation for ^jj^ in terms of the slow-roll param- 
eters e and rj: 



(2.90) 



A quick look at the previous expression reveals that, on superhorizon scales, ip 
behaves as ipk ~ 2{2e — ?7)ifinfV'fc so that |4'?/'/c0o| ^ \2ipkdV/dipo\, whereas on 
subhorizon scales V'fc ~ 0. In view of this, and by making use of the Oi component of 
the perturbed Einstein equations [c.f. Eq. ()2.88j) ]. which may also be written down 
as 

ij + H,niilj = sHl,^ , (2.91) 



we conclude that, on superhorizon scales. 



00 



so that the equation of motion for uj^ is finally given by 







for k ^ ai/inf 



6eH?^fU}k , for -C aHi 



inf ' 



(2.92) 



(2.93) 



This kind of differential equation is much more familiar to us, and we know that 
it can be solved going to conformal time and making the usual change of variables 



Ak = . 



(2.94) 



The resultant equation of motion for is then 



(2.95) 



with v^p defined by 



5e 



1/2 



(2.96) 



The solution for this equation is immediate, based on the results found in Subsec- 
tion 12.4.11 The magnitude of the mode function ook on superhorizon scales is then 
almost constant and given by 



[2(1 -£)]"--§ (1-e^ ^(^'^ 



inf 



k 



r(3/2) V«^mf. 



k 



(2.97) 
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which is used to calculate the spectrum Vs^{k) of perturbations in the infiaton field 
ip by means of Eq. ()2.27|) : 




(2.98) 



The spectrum of perturbations in (p is almost scale-invariant with spectral index n^^ 
given in terms of the slow-roll parameters e and r^,^: 

ns^ = 27]^ - 6e . (2.99) 

Comparing the spectrum obtained [c.f. Eqs. ()2.98|1 and ()2.99j) ] with that for a 
non-dominating scalar field [c.f. Eqs. 1)2.511) and ()2.52j) ]. we see that the backreaction 
of the metric only affects the spectral index of the spectrum of perturbations. The 
amplitude remains the same either the respective field dominates the energy density 
or not. 



2.5.3 The spectrum of ( in the infiaton scenario 

Now we are in position to calculate the spectrum of the curvature perturbation ( 
in the infiaton scenario, based on the results found in the previous subsection. We 
first begin by invoking the definition of ( given in Eq. ()2.13|) in terms of the ip field: 



c 



-ijj - iiinf — 



which, on superhorizon scales, reduces to 



:i+^)//inf- 



--ninf— , 

'^0 



(2.100) 



(2.101) 



where the expression in Eq. ()2.92|1 has been used. 

The spectrum V(^{k) of ( is, in view of the latter, given in terms of the spectrum 
'Psipik) of the perturbations in ip: 



Vs^k) 



(2.102) 



which, according to Eq. ()2.98|) . gives the final expression 



V^k) ^ Al 



k 



inf , 



k 



2 / 



Seirmp 



inf , 



(2.103) 



where, in the second line, the amplitude Ac_ is written down in terms of e. 
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In the inflaton scenario tp is responsible of driving inflation and also of generating 
the required level of primordial perturbations measured by WMAP (I^J ^ 5 x 10~^ 
|207j ). That imposes the following constraint on the Hubble parameter during in- 
flation in terms of e: 

ff* ^ lO^^v^ GeV, (2.104) 

that combined with the present bound e < 0.01 coming from spectral index and 
gravitational waves constraints |2()7j requires < 10^^ GeV [123J. Despite the fact 
that the inflationary energy scale, directly related to the value of H^,, is in the inflaton 
scenario regulated by the parameter e, low-scale inflation may well be obtained but 
only at the expense of a very small e, which in turn requires a high level of fine- 
tuning jl67j (see however Ref. |177j ). As a consequence serious problems appear 
when trying to build successful particle physics inflationary models [ 112| llHOj . The 
relevance of finding such a kind of low-scale inflation models is evident since the 
inflaton field could be identified with one of the MSSM flat directions or one of the 
scalar fields in the SUSY breaking sector (see for example Refs. j^Zl IHll UHl I125t 

Going back to the curvaton scenario, we said that one of the assumptions of 
the model was a negligible curvature perturbation generated by the inflaton. This 
assumption may be quantifled by requiring (r to be, say, at most 1% of the total (, 
which means, from Eq. ()2.103|) . that H^, in the curvaton scenario must satisfy 

< lO^^v^ GeV < 10^2 GeV. (2.105) 

In the following section we will show how such an upper bound on H^, makes the 
detection of gravitational waves an anti-smoking gun for the curvaton scenario. 

2.6 Gravitational waves 

Primordial tensor-type perturbations in spacetime are regarded as gravitational 
waves, being unsourced during inflation and susceptible to be decomposed in a 
polarisation tensor basis. They also propagate in a way that each subhorizon mode 
function follows a harmonic wave equation of motion, i.e. the Klein-Gordon equation 
associated to a massless scalar fleld in Minkowski spacetime. 

The relevant quantity to discuss now is 11^-, which is the tensor component of the 
full perturbed metric tensor, and that, as shown in Subsection l2.2.11 is characterized 
by two degrees of freedom in our three dimensional space as well as by satisfying 
the transversality condition c^ill^- = 0. Hereafter we will write H^- as 2hij as it is 
done in most of the literature. Displaying only the tensor perturbation, the metric 
line element in conformal time reads then 

ds^ = a^{r])[-dri^ + {5ij + 2hij{ri,x))dx'dx^] . (2.106) 

Tensor perturbations are decoupled from their scalar and vector counterparts. The 
line element in Eq. fl2.10(ij) shows then that hij is a gauge invariant quantity. 
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The Einstein-Hilbert action involving hy, and given in a general way by 

SE^-'^jdH-gf'R. (2.107) 

being g the determinant of the g^j_i, metric tensor and R the Ricci scalar, is given as 
a function of the kinetic term associated to hij as obtained from Eq. ()2.106|) : 

Se = -'^J d'x{-gy/'^d,h,,d^h,, . (2.108) 

Notice that no more terms have been added to Eqs. (j2.1()7|l and ()2.1()8|) because no 
tensor-type contributions to the energy-momentum tensor exist during the inflation- 
ary period. The primordial tensor perturbations hij are in consequence unsourced 
so that they propagate freely throughout space following a harmonic (on subhorizon 
scales) wave propagation pattern. 

To clearly show how the hij perturbations propagate, we apply to them the same 
kind of treatment we do to the scalar perturbations in previous sections. We begin 
by decomposing hij in canonically normalised Fourier modes h^ : 

hijiji, x) = — - y exp(ik ■ x) ^ Sijij), \i)hl{r]) + h.c. , (2.109) 

where p = +, x are the two degrees of freedom (polarisation states), and the factors 
Eij {p, k) are the polarisation tensors that satisfy 

Y^kiEi.ip^k) = 0, (2.110) 

i 

^4.(p,k)£,,(p',k) = 26,,,, (2.111) 
Y^e^^%{+,k)e,i{x,k) = -Y^e^^%{x,k)e,,{+,k) = 2^ (2.112) 

ijl ijl I I 

Y,e'^'%{+,k)6ji{+,k) = -J2e^^%{x,k)e,i{x,k) = 0, (2.113) 

ijl ijl 

according to the transversality condition and the properties of the rotational trans- 
formations |129|I134] ^^. Next, we recognize that the mode functions h^., which satisfy 
the Klein-Gordon equation of motion derived from the Einstein action in Eq. ()2.108|) 

hi + 3H,^ihl + ^hl = 0, (2.114) 
are better handled if we rescale them as 



K^-^. (2.115) 



a 



^^hi the expressions of Eqs. (|2.112|l and p. 113(1 . e^^^ is the totally antisymmetric Levi-Civita 
tensor. 
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Thus, the equation of motion for is the same as that for a massless scalar field^^: 



4 + 



Zk = ^. 



(2.116) 



and reduces in the subhorizon hmit {—krj ^ 1) to the Klein-Gordon equation in 
Minkowski spacetime |S2|- We point out that, in deriving the previous expression, 
we have worked in conformal time during a quasi de Sitter stage. The expansion 
parameter aijf) is in this case given by 

«(^) = -TTTTTr^ ' ^2.117) 
where the conformal time f] takes negative values, and the parameter Vh is 



1 2-e 

4^ fl-£)2 



+ e. 



1/2 



(2.118) 



The solution to Eq. ()2.116j) is well known from previous sections (see specifically 
Eq. fl2.5()|l in Subsection 12.4. Ij) . and reduces to the almost time- independent value 



[2{1 - er^~Hl - e] 



k 



T{Vh) i^inf 

r(3/2) [aHi^f^ 



k 



(2.119) 



for the magnitude of the mode function h^. in the superhorizon regime. 

We are interested in the statistical properties of the gravitational waves which 
are well described by the spectrum Vrik) defined by the statistical average 



Y.{KK)^^^\^i + ^2)VT{k) 



(2.120) 



over an ensemble of universes. Here /i^ stands for 

/^k^ = Sijip, k)hl + h.c. . 

mp y 



(2.121) 



To calculate the statistical average during inflation, we must promote the gravita- 
tional wave amplitude to an operator hij by introducing the creation and annihila- 
tion operators and that depend on the polarisation p and wave vector k, and 
satisfy the commutation relation I129j 



pp 



(2.122) 
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That property reflects the masslessness of the graviton (the gravity messenger particle). 
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The gravitational amplitude operator hij{ri,x.) that generalizes Eq. ()2.109|) is then 
given by 

kAv,^) = J j0j-,eM^.^)h^{v), (2.123) 

with 

^k (^) = ^ E hfe mivK + sUp, -m'iv)a'X] ■ (2.124) 

Being all the universes in the ensemble in the vacuum state during inflation, the 
statistical average {h^^_^h^^^) is easily identified with the expectation value {0\h^^_^h^^jO) . 
The spectrum of gravitational perturbations Vrik), defined by Eq. ()2.12()|1 . is then 

VAk) = -£^EE^^.(P,k)4(p',k)/z^(r/)/.r'(^) 

^ y pp 



Ak 



3 



hlW (2.125) 



where one of the properties of the polarisation tensors [c.f. Eq. p.lllll ] has been 
used. Now we can make use of the result in Eq. ()2.119|1 to finally arrive to a definite 
expression for Vrik) on superhorizon scales in terms of if* and e P|l?) ^ITTT|lll2|IT??I]j : 



aHi 



inf , 



irmp 



inf , 



(2.126) 



This nice result shows that the infiationary energy scale, given by H^,, can be 
known from a direct measurement of the amplitude At- Unfortunately, at the 
moment all the efforts to detect gravity waves have been fruitless, leaving only the 
upper bound H^: < 10^^ GeV |207j . In addition, technological restrictions impose 
the lower bound if^, > 10^^ GeV if gravity waves may some day be detected jSSl 
inn inni 12011 12()2[ 12051 . a positive detection would kill then the curvaton scenario 
because, as we had discussed in Subsection 12. 5. 3^ the inflationary energy scale in 
this scenario is required to satisfy < 10^^ GeV to make the inflaton fleld (f not to 
generate enough curvature perturbation. Only in non-slow roll inflationary models 
|170j (speciflcally if rj^ = m'^/SH'^ > 1), the energy scale during inflation could 
be high enough to allow the detection of gravity waves consistent with a negligible 
contribution of ip to (. 

We end up this section by reporting the existence of a consistency relation be- 
tween the curvature perturbation spectrum P^(/c) and the gravitational waves one 
Vrik) in the inflaton scenario pi IHHl ITT^ I190j . The ratio between the amplitudes 
of both spectra [c.f. Eqs. ()2.103|1 and ()2.126|1 ] is given by the slow-roll parameter e: 

'^Tc = ^ = 16e, (2.127) 
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which in turn gives information about the spectral index of Vrik) : 

nT = -2e. (2.128) 

The ratio r^-^ is then consistently related to the spectral index ut, this relation being 
given by the expression 

r-r^ = -8nr. (2.129) 

No consistency relation of this type is encountered in the curvaton scenario or in 
other scenarios for the origin of the large-scale structure in the Universe, although 
it is true that r^^ is always smaller in the multi-component inflationary case jl3n] . 
so that its future confirmation would mean good news for the inflaton scenario. 
Nevertheless, if the consistency relation turns out to be experimentally wrong, that 
does not mean necessarily that the inflationary paradigm is wrong, just that the 
single-field variant is not nature's choosing. Anyway, the non-gaussianity signatures 
associated to each model would add up to the consistency relation in Eq. ()2.129|) . to 
act as powerful discriminators for models that give origin to the primordial energy 
density perturbations (see Chapters and IHl) • 

2.7 Conclusions 

The curvature perturbation is a well defined quantity, gauge-invariant and con- 
served on large scales (if the adiabatic condition is satisfied), that allows us to 
quantify the primordial energy density inhomogeneities produced during inflation. 
The statistical properties of ( are given by the spectrum Vi^{k) whose amplitude and 
spectral index strongly depend on the specific mechanism of production of density 
inhomogeneities. This makes V(^{k) act as a discriminator for the different produc- 
tion mechanisms, at least as the spectral index and the possible relation of 
with the gravitational waves spectrum Vrik) are concerned |112[ 111^(1] . The infla- 
tionary energy scale, given by H^, is well determined by the amplitude of Vrik) 
so that the current upper bound on Aj- leads to < 10^"^ GeV |207j . The en- 
ergy scale may be well below lO^'' GeV, but only at the expense of a high level of 
fine tuning to make the slow- roll parameter e be extremely below 1. This reflects 
how constrained is the inflaton potential in the inflaton scenario, in order to pro- 
duce enough curvature perturbation while driving inflation. As a result the particle 
physics motivated inflationary models are quite unrealistic if we insist that the infla- 
ton field (p has to produce the energy density inhomogeneities |167j . It is here where 
the curvaton scenario comes to the rescue: by requiring cp just to drive inflation, 
the weakly coupled curvaton field a is in charge of giving origin to C. |1HH| lllffl^ . 
The inflationary energy scale is in this case easily lowered so as to possibly associate 
(p with one of the fields appearing in supersymmetric extensions of the Standard 
Model of particle physics. Gravitational waves are in this case so tiny to ever be 
detected, since the current detection technology restricts H^, to be above 10^^ GeV 
(HSl EH ESI I2()H I2()2[ I2()5j . Which scenario for the generation of ( is correct will 
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be determined by future observations. At the moment, we will just try to do our 
best to successfully integrate cosmology and particle physics, being the first step 
the determination of the lower bound for if* in the simplest curvaton model |126j . 
As we will see in the next chapter, in such a model is still quite high, being 
> 10^ GeV, so that a modification to the basic setup is urgently needed. Two 
modifications to the simplest curvaton model will be explored in Chapters 01 and |3] 
jSH 1189] , showing that low scale inflation with if* ~ 1 TeV or lower is possible to 
be obtained. 
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Chapter 3 



Low scale inflation and the 
curvaton mechanism 

3.1 Introduction 

The primordial curvature perturbation ( is generated presumably from the pertur- 
bation of some scalar field, which in turn is generated from the vacuum fluctuation 
during inflation. The scalar field responsible for the primordial curvature perturba- 
tion is traditionally supposed to be the inflaton field (p, i.e. the field responsible for 
the dynamics and, in particular, the end of inflation |112j . This 'inflaton hypothesis' 
is economical, but it is quite difficult to implement and, if many scalar fields exist, 
it presumably is not particularly likely. An alternative is that the curvature pertur- 
bation is generated by the weakly coupled curvaton field a, which could dominate 
(though not necessarily) the energy density before it decays |138| I159j (see also 
Refs. inni IHEl USED- According to this 'curvaton hypothesis', the contribution of 
the inflaton to the curvature perturbation is negligible. This is especially true if 
the energy scale of inflation is much lower than the scale of grand unification, the 
latter scale being the typical requirement of the traditional inflaton hypothesis^. In 
fact, one of the advantages of the curvaton scenario is the relaxation of the con- 
straints on the inflationary energy scale, which alleviates many tuning problems in 
inflation model-building and allows for the construction of realistic, theoretically 
well-motivated inflation models I160| I161j . 

In the simplest version of the curvaton model though, the scale of inflation is 
still required to be quite high corresponding to Hubble parameter i?^, > 10^ GeV 
|126j . The purpose of this chapter is to systematically explore a modification of the 
curvaton model which can instead allow inflation at an even lower scale To 
be specific, we aim for ~ 10^ GeV, which holds if the inflationary potential is 

^Although some exceptions exist (see for example Ref. EZ])- 

^Low scale inflation has also been studied '174' in the context of the 'inhomogeneous reheating 
scenario' [HSIEIEZI (see also Refs. |H3 1142. EHIl EISl EH), where the inhomogeneities in 
the inflaton decay rate during inflation give origin to 
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generated by some mechanism of gravity-mediated supersymmetry breaking which 
holds also in the vacuum. 

We begin by presenting some known bounds in a unified notation. Then we 
consider the possibility that the curvaton mass increases suddenly at some moment 
after the end of infiation but before the onset of the curvaton oscillations 

3.2 The bounds on the scale of inflation 

In this section we present four bounds on the scale of infiation, in terms of three 
parameters which encode possible modifications of the simplest curvaton scenario. 
These bounds have been presented at least implicitly in earlier works jl2()[ 1149^ I175j 
but not in the unified notation that we employ. The advantage of this notation is 
that it will allow us to compare the bounds in various situations, establishing with 
ease which is the most crucial. The three parameters are 

• The ratio e = ct^/cTosc, where a* is the global value of the curvaton field at 
horizon exit and cJosc is its global value when it starts to oscillate. 

• The ratio / = ifosc/^^o-, where ifosc is the Hubble parameter at the start of 
the oscillations and fh^ is the effective curvaton mass after the onset of the 
oscillations. 



• The ratio 6 = J Hose/ where if* is the Hubble parameter a few Hubble 



times after horizon exit. 
3.2.1 Curvaton physics considerations 

The observed value of the nearly scale invariant spectrum of curvature perturbations, 
parameterised by the amplitude A^, is |y4J fa 5x 10~^ |2U7j . In the curvaton scenario 
C is given by [ISHl lISi [cf. Eq. (EI)^ ] 



where Qdec < 1 is the density fraction of the global curvaton energy density over 
the global total energy density of the Universe ptotaio at the time of the decay of the 
curvaton: 



and Ccr is the curvature perturbation of the curvaton field a, which is jlT] [cf. 




C ~ ^dccC 



(3.1) 




(3.2) 



Eq. (I^3HD ] 




(3.3) 



where 'osc' denotes the time when the curvaton oscillations begin and 'dec' denotes 
the time of curvaton decay. 
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In all the cases which we consider, 

(") - (v) ■ 

\ / * \ / osc 

where '*' denotes the epoch when the cosmological scales exit the horizon during 
inflation. The above typically holds true because the curvaton (being a light field) 
is frozen during and after infiation until the onset of its oscillations. However, this 
does not mean that ex* ~ cxosc necessarily. Indeed, in the case of a pseudo Nambu- 
Goldstone boson (PNGB) curvaton with a varying order parameter v, the curvaton 
field is associated with the angular displacement 6 from the minimum of its potential 

as isniEii 

a = V2ve. (3.5) 

Therefore, even though after the end of infiation, 6 remains approximately frozen 
(the angular motion is over damped), we may have e -C 1, where 

(3.6) 



a. 



because [cf. Eq. ()3.5|) ] = eVosc fosc- However, in this case too, for the curvaton 
fractional perturbation we find 



(3.7) 



which agrees nicely with Eq. ()3.4|) . 

Now, for the perturbation of the curvaton we have the following value for the 
amplitude Asa, of the spectrum of perturbations [cf. Eq. 1)2.511) ] 

As.. - § . (3.8) 
Combining Eqs. ()3.fij) and ()3.8j) we find 

(3.9) 

which means that, if the order parameter of a PNGB curvaton grows, the curvaton 
perturbation is amplified by a factor [ST]. 
From Eqs. p.ip and ()3.3|) we have 

o"osc ~ ^doc — 7— = ^dcc— • (3.10) 



Using Eq. ()3.9j) . we can recast the above as 



H ^ dec / \ 

0"osc ~ 7—- (3.11) 
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3.2.2 The main bound on the scale of inflation 



For the density fraction at the onset of the curvaton oscillations we have: 

VPtotaio/o,, \rnpj 

where 

f^^, (3.13) 

and we used that (pao)osc - |"^a'^osc and (ptotaio)osc = SH^^^mj,. Here, denotes 
the effective mass of the curvaton after the onset of its oscillations. In the basic 
setup of the curvaton hypothesis this effective mass is the bare mass If this 
is the case then = m„ ~ Hose (i-e. / ^ 1). However, in the heavy curvaton 
scenario, the mass of the curvaton is supposed to be suddenly incremented at some 
time after the end of the inflationary epoch due to a coupling of the form Xx^cr^ with 
a field x which acquires a large vacuum expectation value (VEV) at some time after 
the end of inflation jSIl 1123 IIH3- In this case rhl = ml + X{x)^ ^ X{xf > H^^c 
(i.e. / « 1). 

Now, we need to consider separately the cases when the curvaton decays before 
it dominates the Universe {fldec <^ 1) or after it does so {^dec ~ !)• Note, that the 
WMAP constraints on non-gaussianity in the CMB [102J impose a lower bound on 
fidec, which allows the range ^138j [c.f. Eqs. (ITHHl) and 

0.01 < ^]dcc < 1 . (3.14) 

Because of the above bound we might require that the density ratio Po-/ptotai grows 
substantially after the end of inflation. Typically, in the curvaton scenario this does 
indeed take place after the curvaton begins oscillating, but only if the curvaton 
oscillates in a quadratic potential during the radiation era. As it was shown in 
Ref. if the curvaton oscillates in a quartic or even higher order potential, its 
density ratio does not increase with time (it may well decrease instead) and satisfying 
the bound in Eq. ()3.14j) might be very hard. Due to this fact, in the following, we 
assume that the period of oscillations occurs in the radiation era with a quadratic 
potential. Hence, we consider that ifosc < r^f, being Finf the infiaton decay rate. 

Suppose, at first, that the curvaton decays before dominating the density of the 
Universe so that fl^ec 1- Assuming that the curvaton oscillates in a quadratic po- 
tential, during the radiation epoch, its density fraction grows as Po-o/Ptotaio H^^^"^. 
Therefore, at curvaton decay we have 

fhlal 



^dec ~ J's/r 3/2 ' (3-15) 

J dec -Hose Tflp 



where we used Eq. fl3.12|l and also that (ptotaio)dec ~ T^eci with Tdec being the radi- 
ation temperature just after the curvaton decay. Using Eq. ()3.11|) the above can be 
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recast as 

(3.16) 

where we used that T^^^ ~ H^^c'^np- 

Now, suppose that the curvaton decays after it dominates the Universe so that 
^dec ~ 1- Since (po- /ptotai)dom — 1 by definition, using again that, during the radia- 
tion epoch, Po-o/Ptotaio ^ H~^/'^ and in view of Eq. ()3.12|) . we obtain 

i^dom-i^cj-'f — (3.17) 

\mp/ 

where 'dom' denotes the time of curvaton domination^. Employing again Eq. ()3.1H) . 
the above can be written as 

H.^7,eA^fmp(^)''\ (3.18) 

V -Hose / 

Combining Eqs. ()3.1fj|l and ()3.18j) we find that, in all cases 

r. neAJ^ ( ^^^^^^1^21^11^) . (3.19) 

V"dec \ "osc / 



This can be rewritten as 

1/5 



H. ~ (#^) ( °""'^-°-^-'"' ) (n.A,f)"Hn,,miy'^ , (3.20) 

V-Hosc/ \ -nBBN / 



or equivalently (using V*^^ ~ y/H^mp) 
V:'' - (#^)''" f """^^ l {^eA,fr\T^^^m%f'^ , (3.21) 

V-nosc/ \ -HbBN / 

where 'BBN' denotes the epoch of Big-Bang Nucleosynthesis (BBN) (Tbbn ~ 1 MeV). 
Now, according to Eq. ()3.14j) we have fidec < 1. Also, we require that the curva- 
ton decays before BBN, i.e. i^doc > -f^BBN. Hence, the above provides the following 
bounds 



> {TreAjY/\T^^^m%Yl^ ~ (e/)^/^ x 10^ GeV, 
V}'^ > (7reA^/)2/5(TBBNm|,)V5 ^ (e/)2/5 X 1012 GeV. 



(3.22) 



■^Here we define Hdom by Hdom = Hcq, where Hcq is the Hubble parameter at the time when 
the global curvaton energy density po-o makes equal to the global radiation energy density p^o ■ 
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In the standard setup of the curvaton scenario e = / = 1 and the above bounds do 
not allow inflation at low energy scales to take place |126j . However, we see that if 
either e or / are much smaller than unity the lower bound on the inflationary scale 
can be substantially relaxed and low scale inflation can be accommodated^. Still, 
though, there are more bounds to be considered. 



3.2.3 Other bounds related to curvaton decay 



Firstly, let us consider the bound coming from the fact that the decay rate of the 
curvaton field cannot be arbitrarily small. Indeed, in view of the fact that the 
curvaton interactions are at least of gravitational strength, we find the following 
decay rate for the curvaton 



< 



(3.23) 



where 7^ > 1- 

Suppose, at first, that the curvaton decays after the onset of its oscillations, as 
in the basic setup of the curvaton scenario. In this case, < Hose and i^dec = Tcr- 
Using the fact that maxjifdom. To-} > To-, Eq- ()3.23|) suggests 



max{i^dom,-f^dec} . r-ifm^\ 

7} >lJ — ) 

Hqsc \mpj 

Including the above into Eq. ()3.19|) the latter becomes 

mp /Ho, 



dec 



(3.24) 



(3.25) 



which results in the bounds 



> {neA^Yy/f 6^ nip ~ e^y/f 6^ x 10" GeV, 
K^/^ > neAJ^/^Smp ~ x 10^^ GeV, 



(3.26) 



where we have defined 



6 = 



I Ho 



H^ 



(3.27) 



which must be really small in order to reduce the bounds in Eq. ()3.22j) to satisfactory 
levels. In the case of a PNGB curvaton we see that the bounds in Eq. ()3.2(ij) are 
drastically reduced with e, compared with the bounds in Eq. ()3.22|1 . 

'^The relevance of a low e makes evident in Ref. where the scenario of the curvaton as 
a PNGB is studied. The specific explored model refers to a PNGB whose order parameter v is 
increased after the cosmological scales exit the horizon during inflation, but before the onset of the 
curvaton oscillations. That makes e very small. 
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Now, provided we demand that the curvaton field does not itself result in a 
period of infiation, we see that the curvaton cannot dominate the Universe before 
the onset of its oscillations. This results into the constraint 



^ Ptotalo 



dec 



{neA(^)mp 



(3.28) 



where we used Eqs. (jHHH), ()3.12|1 . ()3.13|) and ()3.27|1 . Inserting the above into 
Eq. ()3.25|) we obtain 

(3.29) 



a 



dec 



which results in the bounds 



H, > (neA^fd^mp ~ e^^^ x 10^ GeV, 
K'/' > {ireA^f/^^mp ~ e^/^S^ x lO^^ GeV. 



(3.30) 



A similar bound is reached with the use of the upper bound on rricr 



(3.31) 



which comes from V„ < H^sc and Eq. ()3.23|) . instead of the bound in Eq. ()3.28p . 
Inserting the above into Eq. ()3.25|) one finds [cf. Eq. ()3.29|) ] 



H,>l.{ixeA^f5'^ 



n 



(3.32) 



dec 



which, again, results in the bound in Eq. 



as it was suggested in Ref. |175j . 
= 1 and also i^osc — minjifpt, rrio-}. 



In the heavy curvaton scenario we have e 
where Hpt corresponds to the phase transition which increases the effective mass 
of the curvaton. Then, if 5 ^ 1, the bounds in Eq. ()3.30p are not possible to be 
relaxed below the standard case discussed in Ref. |126j despite the fact that we may 
have / -C 1 in Eqs. ()3.22|1 and ()3.2(ij) . Therefore, in the heavy curvaton scenario we 
require 5 -C 1, i.e. the onset of the curvaton oscillations has to occur much later 
than the end of inflation so that ^ i/osc > Ttr |149j . In this case, as can be seen 
in Eq. ()3.30p . it is easy to lower the bound on the infiationary scale even for a not- 
so-small value of 6. This is a very nice feature of this scenario. Note also, that in the 
case of a PNGB curvaton i^osc ~ '^a ^ and 6 is very small necessarily. Because, 
in this case, e <C 1, it is straightforward to see that the bounds in Eq. ()3.3Up are 
much weaker than the bounds in Eq. ()3.22|) . 

As it was pointed out in Ref. |175j . the sudden increment in the curvaton mass 
might lead to a growth in the curvaton decay rate enough for > Hpt. This 
would force the curvaton to decay immediately and we can write Hose ~ Hpt ~ H^ee- 
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Obviously, in this case we cannot have H^^c < -f^dom and there is no period when 
Pao/P totaio oc This means that (pao/Ptotaio)osc ~ ^dec- Using Eqs. and 

).12|) it is easy to find 

(3.33) 



~ lie A J 



mp 



which results in the following bounds 



dec 



H, > TieAJmp ef x 10^^ GeV 



V;^^^ > J'KeA(;fmp ~ (e/)i/2 x lO^^ GeV. 



(3.34) 



It is evident that the above bounds may challenge the WMAP constraint for the 
curvaton scenario jlH] leading to excessive curvature perturbations from the infiaton 
field if e and/or / are not much smaller than unity. 

The bounds in Eqs. ()3.22p . ()3.26|) . and ()3.30p provide the basis for our investiga- 
tion jni], leaving the fourth bound in Eq. ()3.34|1 to be considered in the next chapter 
jl89j . As a matter of completeness we have considered all the other possible bounds 
coming from the requirements that < rhcr and ifdec > -f^BBN- We have found that 
these bounds lead to consistent and/or weaker constraints than the above four. 



3.3 The case of a heavy curvaton 

In this section we are going to consider the so called 'heavy curvaton scenario' 
where an increment in the curvaton mass, at some moment after the end of infiation 
but before the onset of the curvaton oscillations, leads to a huge decrease of the 
infiationary scale through the attainment of a very small parameter 6 [cf. Eq. ()3.30|) ] . 
We will do so by the implementation of a second infiationary period following the 
idea first presented in Ref. jl49p . We identify this second infiationary period as 
the thermal infiation one which triggers the increment in the curvaton mass when 
the fiaton field, that responsible for the generation of the thermal infiation era, rolls 
down towards the minimum of the potential. 

3.3.1 The thermal inflation model 

Thermal infiation was introduced as a very nice mechanism to get rid of some un- 
wanted relics that the main infiationary epoch is not able to dilute, without affecting 
the density perturbations generated during ordinary inflation. As its name suggests, 

^Note however that any post-inflationary phase transition could serve the purpose of giving an 
efi^ective mass to the curvaton field. 
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thermal inflation relies on the finite-temperature effects on the flaton scalar poten- 
tial. A flaton field x could be defined as a field with mass and vacuum expec- 
tation value M 3> |136[ I137j . More specifically, a fiaton field is a MSSM fiat 
direction lifted by non-renormalisable terms. SUSY breaking provides soft terms 
which create a large vacuum expectation value because the absence of quartic terms 
in the potential. The possible candidates for a fiaton field within particle physics 
are either one of the many expected gauge singlets in string theory |172j or the GUT 
Higgs (which is a scalar field charged under the GUT gauge symmetry but neutral 
under the Standard Model one) with ~ 10^ GeV and M ~ 10^^ GeV [TT^ ^. 
After the period of reheating following the main infiationary epoch, the thermal 
background modifies the fiaton potential V trapping the fiaton field at the origin 
and preventing it to roll-down towards M [13, 108J. At this stage the total energy 
density ptotai and pressure Ptotai are 



r , 



Ptotai = V + p. 
Ptotai = -V + ^Pr, (3.35) 

making the condition for thermal infiation, ptotaio + SPtotaio < 0, valid when the 
global thermal energy density p^o falls below the height of the potential Vh, which 
corresponds to a temperature of roughly Vj^^'^. Thermal infiation ends when the 
finite temperature becomes ineffective at confining the field, at a temperature of 
order m^, so the number of e-folds this infiationary period lasts is 

/ «end \ , M start \ I 




N = \n{ = In ~ In ~ - In ~ 10 . (3.36) 

Vastart/ VTcnd/ \ "^X / ^ V"^x/ 



Here we have used the fact that, in a fiaton potential of the form 

oo 

V = V,-iml- gT')\x\' + E Kmf^lxl'''^' , (3-37) 

n=l 

where the nth term dominates: 

ml = 2{n + l)ml, (3.38) 

" ^ 'ml, (3.39) 



^2n+2^-2n _ [2(n + 1) (n + 2) A„ 

Vh = [2{n + 2)]-'mlM^ 



(3.40) 



Note that the gT"^ contribution to the effective mass of the flaton fleld stands for 
the effect of the thermal background, which changes the slope of the potential in the 
X direction and traps the flaton fleld at the origin of the potential [T3t llOSj . It is 
worthwhile to mention that the potential is stabilized by non-renormalisable terms, 

^Note, though, that in some GUT models there are additional Higgs fields with much smaller 

vevs EniEiiiiiiiinHi- 
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with dimensionless couplings A„ ~ 1 to make the theory vahd up to the Planck scale. 
Notice also that the A4|x|^ term is absent in the potential; otherwise, the vacuum 
expectation value M would not be much bigger than m^^, spoiling the suppression 
of the unwanted relics. 

Before embedding the thermal inflation epoch and the curvaton mechanism into a 
single model, we want to clarify some issues about the nature of the interactions that 
produce the thermal background. If the flaton is a GUT Higgs, it is coupled with 
those fields charged under the GUT gauge symmetry, in particular with those the 
inflaton field decays into. That collection of particles makes the thermal background, 
and its interaction with the flaton field produces the thermal correction. If the 
fiaton field is a gauge singlet it still can be coupled, via Yukawa coupling terms, 
with some other fields, possibly in a hidden sector, that the inflaton field decays 
into. Again, a thermal correction is generated. The actual interactions and decay 
rate are not important as the main objective of this chapter and the next one is to 
obtain some particle physics model-independent information about the possibility 
of reconciling low scale inflation with the curvaton mechanism, in a scenario that 
involves a second period of inflation (thermal inflation), without going into the 
details of the identification of all the relevant fields (inflaton, flaton, and curvaton) 
in the framework of a particle physics model (GUT theories, MSSM, etc.), which 
would make the results highly particle physics model-dependent. The flaton could 
be either a gauge singlet or the GUT Higgs. In the former case the flaton can 
be coupled with some other fields that the inflaton field decays into, via Yukawa 
coupling terms, and the specific interactions would be known once we choose what 
of the many gauge singlets expected in string theory is the flaton. In the latter case 
the interactions in the GUT models are already known. The specific interactions 
are important of course, but there are so many possibilities that the general result 
would be hidden behind the characteristics associated to any definite particle physics 
model. 

Having discussed the nature of the flaton interactions, and guided by the result 
in Ref. |149j . we proceed to implement a second inflationary stage into the curvaton 
scenario in order to lower the main inflationary energy scale. If this second epoch 
of inflation is the thermal inflation one devised in Refs. |1U8| I136| I137j we would be 
solving not only the issue of the ordinary inflation energy scale but also the moduli 
problem present in the standard cosmology [TUl IH^ WA 137]. 

In the curvaton model supplemented by a thermal inflation epoch two fields 
X and (J, which we identify as the flaton and the curvaton fields respectively, are 
embedded into the radiation background left by the inflaton decay. It is assumed 
that the curvaton field could be either a gauge singlet |172j . the Peccei-Quinn field 
im, a PNGB 001 Eni EE], or a MSSM flat direction |SHl EHl E3 EH US 1121 E21 

I154| 1158^ I173| 1176] ■ and has just a quadratic interaction with the flaton one so 
that the unperturbed component is frozen at some value a* until the time when the 
flaton field is released from the origin and rolls down towards the minimum of the 
potential. This in turn signals the end of the thermal inflation era and the beginning 
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of the oscillations of the curvaton field around the minimum of its quadratic potential 
|47j . The flaton field, in addition to the non renormalisable terms with A„ ~ 1 that 
stabilize the potential and make its slope in the x direction be very flat, presents 
a quadratic interaction with the curvaton field. The complete expression for the 
potential is 

oo 
n=l 

where ~ 10'^ GeV due to the soft SUSY contributions in a gravity mediated 
SUSY breaking scheme. Under these circumstances the condition for an inflationary 
period, ptotaio + SPtotaio < 0, is satisfied when the global thermal energy density p^o 
falls below Vh. Of course, this period of thermal inflation ends when the effect of the 
thermal background becomes unimportant, at a temperature T ~ m-^, liberating 
the fiaton field to roll down towards the minimum of the potential and letting it 
acquire a large vacuum expectation value M given by: 



1/2 



M ~ . (3.42) 



The evolution of the energy densities associated to the different fluids in this case 
are sketched in Fig. 13.11 

Let's assume that the usual inflation and its corresponding reheating have al- 
ready happened, so that the flaton and the curvaton fields are embedded into a 
radiation bath. Therefore, even when the minimum of the potential is located at 
Xo = ^x(^*) 7^ ctq = 0, X is trapped at the origin because of the finite- 
temperature effects and |(Jo| = a* 7^ because ra„ < H < H^. Thus, the value of 
the scalar potential at this stage is: 



V^(Xo = 0,ao = (T,) = V^h + mX, (3-43) 



with 



ml = 2(n + l)(m2 -AlaoH, (3.44) 
Mf+2^p2n ^ [{n + 2)X^]-\ml-X\ao\'), (3.45) 
Vh = [2(n + 2)]-i(mjM2) |,„=o • (3.46) 

When the thermal energy density falls below Vh + m'^al thermal inflation begins. 
This period lasts until the temperature is of the order the effective mass of the 
flaton field which is fh^ = (m^ — XalY^"^. Note that Aa^ < because otherwise 
there is no thermal inflation. Then, we obtain a first constraint on the value of the 
parameter A: 

^ ml 10-2 GeV' ,^ 
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Figure 3.1: Evolution of the energy densities in the thermal inflation model where 
the curvaton field a has some time to oscillate before decaying 51 . The continuous 
line corresponds to the global radiation energy density pro; the dashed dotted line 
corresponds to the global flaton energy density p;^^, and the dashed line corresponds 
to the global curvaton energy density Po-o- The horizontal axis represents the ex- 
pansion parameter a. From the left to A radiation dominates the energy density, 
although it decreases following p^^ oc a~^. At this stage the flaton and curvaton fields 
X and a are frozen at Xo = and ctq = cr* making their energy densities constants. 
When reaches Vh at A, thermal inflation begins. The thermal inflation period 
lasts until B when the temperature T becomes of the order of the flaton mass m^. 
Thermal inflation stage is portrade by the dashed region. After thermal inflation 
ends, the parametric resonance process transforms a substantial fraction of p-^ into 
pr |2Sl IHHl inni I212j . The flaton field is liberated by this time and begins oscillating 
around the minimum of its potential, behaving then as a matter fluid with p^^ oc a~^. 
The curvaton field increments suddenly its mass m^- at S as a result of the oscilla- 
tions of X around the vacuum expectation value M, which yields to a much bigger, 
but still subdominant, Po-. The increment is enough for the effective curvaton mass 
rhfj to overtake H^t (the Hubble parameter at B) so that ctq gets unfrozen and starts 
oscillating around (Jq = 0. The curvaton field behaves then as a matter fluid so that 
Po-p oc a~^. By the time C, x already dominates the energy density before decaying 
into radiation. The curvaton field continues to oscillate until D when it decays into 
radiation after having come to dominate (though not necessarily) the total energy 
density. The curvature perturbation is transfered to the radiation at this moment 
due to the decay of a. 
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where we have used the Eq. fl3.1H) and focused on ~ 10'^ GeV which comes from 
the gravity-mediated SUSY breaking contributions. 

When thermal inflation ends the thermal energy density is no longer dominant. 
The Hubble parameter at the end of thermal inflation is then associated to the 
energy density coming from the curvaton and the flaton fields: 

= 3^^ ^ ' ^^-^^^ 

so that 

i/osc-lO-^'M, (3.49) 
and therefore the parameter / [cf. Eq. ()3.13j) ] is 

Z^:^^ 10-16^1^ (3.50) 

where M = 1 0-0=0 is somewhere in the range 10'^ GeV <^ M < 10^^ GeV. 

With this so-low value for the Hubble parameter at the end of thermal inflation, 
the parameter 6 [cf. Eq. fl3.27|l ] is 




5-10-^/— , (3.51) 



so that the bounds in Eqs. ()3.22j] and ()3.2(j|l become^: 



> 10~'GeV^jj^, (3.52) 

/Vf3/4 

> lO-'GeV'/^ (3.53) 

The effective mass of the curvaton field after the end of thermal inflation, i.e., when 
X = and a = are the average over the oscillations of the flaton and the curvaton 
fields, is 

= {ml + XM^y/^ . (3.54) 

Note that we are focusing in the case of a final curvaton decay rate smaller 
than the Hubble parameter at the beginning of the oscillations Hose- This is to 
allow the curvaton field to decay after the flaton field so that we can keep working 
in the simplest curvaton scenario where the curvaton field oscillates in a radiation 
background ^[1331123 (see Fig. EH). 

Making use of the constraint in Eq. ()3.47|) and the expression in Eq. ()3.54j) . 
and taking into account that the bare curvaton mass rricr is smaller than the Hubble 



^The bound in Eq. (|3.30() is consistent with low scale inflation in view of AI < 10^^ GeV. Notice 
also that, in the heavy curvaton mechanism, e = 1 because there is no amplification of the curvaton 
perturbations. 
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parameter i^osc at the end of thermal inflation, we obtain an upper bound on the 
effective mass of the curvaton field: 

M 

< 10-^ GeV . (3.55) 

-'^ * i i dec 

This bound, when applied to Eq. ()3.52|) . is consistent with low scale inflation. 
When Eq. ()3.55|) is applied to Eq. ()3.53|1 . we obtain a lower bound for which is 
consistent too with low-energy scale inflation since M < 10^^ GeV: 

> 10-9 GeV^/''M^/^n^^J^ . (3.56) 

The last inequality is stronger than that of Eq. ()3.52|1 only while the effective mass 
of the curvaton field is 

> 10^ GeV''/''M'/''n%l' . (3.57) 
Otherwise, we still need to consider the expression in Eq. ()3.52|) . 



3.3.2 Required parameter space 

Once we have checked the viability of a low-energy scale infiation we proceed to 
investigate the required range of values for the parameters of the Lagrangian. Re- 
member that we are going to focus on the gravity-mediated SUSY breaking scheme 
where the Hubble parameter during infiation is ~ rns/2 ~ 10'^ GeV. After thermal 
infiation has ended, the fiaton and curvaton fields start to oscillate, eventually de- 
caying into thermalised radiation (see Fig. 13. Ij) . The decay process is distinguished 
by the decay rate. The field with the biggest decay rate will decay first. The fiaton 
and curvaton decay rates are given by 

r,^7x^ and (3.58) 

with 7x ^ 1 |in8| I136[ I137j and 7^ > 1. Since we like the curvaton mechanism not 
to be modified, the fiaton field must decay well before the curvaton decay. This 
requires 

o o ml 10^*^ GeV^ 
'^'«<W2 (3-59) 

Now, using the expression in Eq. (HTH^ . which is relevant for rh^ < 10^ GeV^^/^^M^/^fi^^c^ 
[cf. Eq. fl3.57p ]. we require 

> IQ-^^Mn-^/'^ , (3.60) 

in order to obtain low-energy scale infiation. Note that, combining the above with 
Eq. (Unni), we find 

/ < IO-^a/tW^ 1' (3-61) 
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as required by the heavy curvaton scenario. Similarly to the above, using the 
expression m Eq. (ITHHIl . which is relevant for > 10^ GeV^^^^^M^/^^n%l^ 
[cf. Eq. (I3.57|) ]. we require 

> 10~^°GeV~^M^fidcc • (3-62) 

Thus, for values of rhcr less than 10^ GeV^°^^^M^/^^fij^^^ the required range of 
values for rn-o- is^: 

IQ-^^M <m,< 10^ GeV^°/"Mi/i^ , (3.63) 

where the lower bound comes from Eq. ()3.60|) . The vacuum expectation value M is 
in the range 

10^2 GeV < M < 10^^ GeV , (3.64) 

where the lower bound comes from the solution to the moduli problem as we will see 
later, and the upper bound comes from Eq. ()3.63|) . On the other hand, for values 
of bigger than 10^ GeVi°/"Mi/"r]d^c^ the required range of values for is: 

max{102 GeVi°/"Mi/i\ 10"^° GeV-^M^} <m,< 10^^ GeV^/V^'^' , (3-65) 

where we have used Eqs. ()3.59p and ()3.62p . and M can be, a priori, up to rap. 
We have considered all the other possible constraints on rh^ and found they are 
irrelevant compared with those in Eq. ()3.63p and Eq. ()3.65|) . 

Fig. l3.2l shows the required parameter space A vs M (grey region) as a logarithmic 
plot. We have made use of the definition of the curvaton effective mass m„ in terms 
of the coupling constant A and the vacuum expectation value M [cf. Eq. ()3.54p ]: 

ml ^ , (3.66) 

and the required parameter space rhcr vs M studied before. Note that for values of 
M higher than ~ 10^^ GeV it is impossible to satisfy Eq. ()3.65|) . so our final range 
for M is 

10^2 GeV < M < 10^^ GeV . (3.67) 
The required values for A, according to Fig. 13.21 

10~22 < A < 10~i° , (3.68) 

are in agreement with the upper bound in the Eq. ()3.47j) : 

lO^GeV^ 4 
^<^f2 10 ' (3-69) 



and with the lower bound 



3mp 



^The bound in Eq. is weaker than rh„ < 10^ GeV^"^" Af^/" within the allowed range 

for M (see Fig. ESI- 
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log A 




12 125 13 13.5 14 14.5 15 



log (M /GeV) 

Figure 3.2: Required A — M parameter space (grey region) as a logarithmic plot. 
The two lines in the middle of the graph correspond to the limits in Eq. ()3.63|) 
which are valid up to the meeting point of the three lowest lines. The slanting lines 
correspond to the limits in Eq. ()3.65p . Note that it is impossible to satisfy the 
conditions in Eq. ()3.(i5j) beyond the meeting point of the highest and the lowest 
lines. Such a range of values for M suggests that the flaton field cannot be the GUT 
Higgs field studied in Ref. |136j . The flaton field should be, therefore, one of the 
gauge singlets present in string theory ;172] . 



which follows from > Hose- 

In view of the allowed range of values for M [c.f. Eq. ()3.67p ]. we conclude that 
our flaton field cannot be the GUT Higgs field investigated in Ref. |136j . We must 
remember however that in some other GUT models there are additional Higgs fields 
with much smaller vevs [33, iM* '10^ I198j so they are still good flaton candidates. 
The flaton field as a gauge singlet in string theory |172j remains as a viable option. 

Once we have found the required parameter space for A we must do the same for 
the other relevant parameter of the Lagrangian: the bare mass of the curvaton m„. 
The only bound on rricr is 

< i/osc ~ lO-^'M , (3.71) 

which is related to the fact that the oscillations of the curvaton around the minimum 
begin due to the sudden increment in the curvaton mass at the end of thermal 
inflation. That means, in view of Eq. ()3.67|) . that 

<10"^GeV. (3.72) 
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Such a small value for mg-, taking into account the soft supersymmetric contributions 
of order the gravitino mass for any scalar field which is not protected by a global 
symmetry, leads us to point a PNGB as a viable curvaton candidate [101 (HDl iHH| • 

Finally, we still need to understand the lower bound M > 10^^ GeV. To do that, 
we must study the solution to the moduli problem^. 



3.3.3 Solution to the moduli problem 

Among the unwanted relics that the inflationary epoch is not able to dilute are the 
moduh [nil Eni 1121 ISni ISI] • ModuH fields are fiaton fields with a vacuum expectation 
value of order the Planck mass. The decays of the fiaton and the curvaton fields 
increment the entropy, so that the big-bang moduli abundance, defined as that 
produced before thermal inflation and given by |137j 

^ lOrrip 

where $ is the vacuum expectation value of the moduli fields, gets suppressed by 
three factors. One is 

. ^, 9*iTpR) Tlp. 

^""-^:(^^' ^'-''^ 

due to the parametric resonance process [2S1 UHl IHH I212j following the end of the 
thermal inflation era, where the are the total internal particle degrees of freedom, 
TpR is the temperature just after the period of preheating, and Tc is the temperature 
at the end of thermal inflation; another is 



(27r745)(7,(Tp^)r|^ 



due to the fiaton decay, where is the temperature just after the decay^°, and 
j3 is the fraction of the total energy density left in the flatons by the parametric 
resonance process and the increment in the energy density of the curvaton {j3 < 1); 
the other is 

due to the curvaton decay, where Tjec is the associated reheating temperature which 
must be bigger than 1 MeV not to disturb the nucleosynthesis process^^. This 



^In the following subsection we correct one mistake in Ref. which led to a reduced parameter 
space for nia-- Conclusions are different of course, but they are now more positive than before. 

i^This is assuming that the fiaton has come to dominate the energy density just before decaying 
(see Fig. EIJ. 

i^We have assumed that does not change appreciably from the time when T ~ Tc to the 
time when T = T^. This is a good approximation since T^. 
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enhancement in the entropy depends on the temperature just after the flaton decay 

1013 GeV^ ^ ^ ^ 

which is obtained by setting ^ H and assuming that the flaton decay products 
thermahse promptly. Thus, the abundance of the big-bang moduli after thermal 
inflation is: 



2 / rji ^ , rji ^ 3 



> lO^GeVA-'ArSJl — I 1-^1(^1 X 



2 



The lower bound 



^ \mp) UMeVy Vm$/ 
f y/^ 1 ( mlM^ \ 1 /lO^GeV 

is obtained when taking into account the restriction n^/s < 10^^^ coming from 
nucleosynthesis . This is a weaker bound on A than those presented in Fig. 13.21 
Let's have a look at the thermal inflation moduli abundance defined as that 
produced during the preheating stage following the end of the thermal inflation era 

n^^ $^l^feVlOm|m|, <l>lVhT^Tdec 




/lO^GeVX f Vh \ 1 (lO^GeVV 

Here $7- corresponds to the vacuum expectation value of the thermal moduli fields. 
To suppress the thermal inflation moduli at the required level n$y/s < 10^^^ we 
require 

, ^ lO^GeV^ 2 

^ ^ -J^^i ■ (3.81) 

Again this is a weaker bound on A than those in Fig. 13.21 

The Eqs. fl3.78|l and ()3.80|1 give us information about the necessary conditions 
for the suppression of the big-bang and thermal inflation moduli, but they are based 
on the unknown parameters M and A. Since we still need to know if the range 
M < 10^^ GeV, required to obtain a low-energy scale inflation, is not forbidden by 
the requirements coming from the solution to the moduli problem, we must find a 
A-independent relation on M. This relation can be found noting that the increment 
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in the entropy due to the curvaton decay [c.f. Eq. ()3.76|) ] can be written in an 
alternative way: 

.1/4 



A. 



dec I 



SO the abundance of big-bang moduh after thermal inflation is: 



(3.82) 



> 102^GeV^M-V'(l-r]acc)'/'(;^)" 



X 



This means that 



M>10^2QgV, (3.84) 

to satisfy n^^/s < 10^^^. This is the lower bound on M we have used throughout 
this chapter. A similar treatment to the abundance of thermal inflation moduli 
[c.f. Eq. ()3.80|) ] leads to the bound M < 10^^ GeV, which is weaker than that ob- 
tained in Fig. 13.21 

Of course we might have considered the scenario where there are no moduli 
fields at all. Without the introduction of the moduli problem Eq. ()3.84j) becomes 
unnecessary. This does not help for the improvement of the required range of values 
for m„ but it does for A as the lower bound on M in Eq. ()3.84|) becomes replaced 
by M ^ 10^ GeV, which comes from the definition of the fiaton fields. In this way 
the range of values for M extends to smaller values well below 10^^ GeV until the 
coupling constant A eventually reaches the lower bound 10~^. 

The introduction of a period of thermal inflation into our curvaton scenario, 
sketched in Fig. 13.11 has helped us not only to lower the energy scale of the main 
inflationary epoch, but also to solve the moduli problem JOl OHl 1121 IMl ISZ| still 
present after ordinary inflation. The required parameter space for A has been plotted 
in Fig. 13. 2[ and the vacuum expectation value for the fiaton field has been showed 
to be in the range 10^^ GeV < M < 10^^ GeV. Our fiaton field, in view of the 
allowed range of values for M, should be one of the gauge singlets present in string 
theory |172j . The upper bound on mg- [c.f. Eq. f)3.72|) ]. m„ < 10^^ GeV, suggests 
the curvaton field could be a PNGB [101 EOl |Hni- This is because in the presence 
of supergravity all the scalar fields, whose masses are not protected by a global 
symmetry, acquire soft masses of the order of the gravitino mass if if < m3/2, and 
contributions contributions to the squared mass of order if if > m3/2 except 
during the radiation dominated era |128j . The smallness of the curvaton mass is 
in turn because of the very small value for iiosc- The parameter ifosc is directly 
proportional to M, so the bigger M is, the more possible to obtain a range of values 
for rricr compatible with the soft supersymmetric contributions. We, in the next 
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chapter, will look for a mechanism to improve the required range of values for the 
bare mass m„ and the coupling constant A in presence of the moduli problem jl89j . 

3.4 Conclusions 

We have presented a different type of curvaton scenario [ST], in which the scale of 
inflation can be much lower than H^, ~ 10^ GeV, which is the default lower bound 
for the standard curvaton model |126j . This scenario considers a curvaton, whose 
mass, being appropriately Higgsed, is substantially enlarged at a phase transition 
after the end of inflation ('heavy curvaton'). We have shown that this mechanism 
is indeed able to accommodate inflation scales as low as iJ* ~ 1 TeV or even lower. 

We have implemented the idea of a thermal inflation epoch, introduced in Refs. 
jl()8| llH7j to solve the moduli problem, as a second inflationary period nec- 
essary to lower the energy scale of the main inflationary stage. In our model, a 
flaton field x with bare mass coming from soft supersymmetric contributions and 
vacuum expectation value in the range 10^^ GeV < M < 10^^ GeV (i.e. one of the 
gauge singlets in string theory |172j ). is held at the origin of the scalar potential by 
finite-temperature effects. These effects are associated to the thermal background 
created by the main reheating epoch. When temperature falls below Vh thermal 
inflation begins. This period of thermal inflation lasts around ten e-folds until the 
temperature falls below liberating the flaton field to roll away towards the min- 
imum of the potential. The curvaton field is coupled to the fiaton one, through a 
coupling constant A in the range 10~^^ ^ ^ 10~^°, so its mass is largely increased 
at the end of thermal infiation. This increment is enough to lower the bound on 
to satisfactory levels, without sending the non-gaussianity constraint to the limit. 
However, the energy scale of the thermal infiation epoch is very small, requiring in 
turn a bare mass for the curvaton field of at most 10~^ — 10"^ GeV. Taking into 
account the soft supersymmetric contributions to mo-, the required smallness of m^j 
points toward using a PNGB curvaton |10l ISHl IHE] to achieve low-scale infiation. 

The type of mechanism that we presented is not completely compelling. It 
suffers from the problem that the mass of the curvaton before oscillation, as well as 
its coupling, have to be much smaller than one would expect. In the next chapter 
|189j it will be shown how this tuning problem can be at least alleviated. 
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Chapter 4 



Low scale inflation and the 
immediate heavy curvaton decay 

4.1 Introduction 

Low scale inflation is desirable in order to identify the inflaton field with one of the 
MSSM flat directions ^1 1^ or with one of the fields appearing in the SUSY break- 
ing sector, giving the inflaton a much deeper particle physics root. In contrast low 
scale inflation is not desirable because it makes very difficult the generation of the 
adiabatic perturbations by the inflaton, leading to multiple fine-tuning and model- 
building problems, unless the curvaton mechanism is invoked |138l I139| I159j (see 
also Refs. pUl I116| I156j ). With the aim of generating the curvature perturbation 
that gives origin to the large-scale structure in the observable universe, the curvaton 
mechanism has appeared as a nice and plausible option and a lot of research has 
been devoted to its study. Making the curvaton mechanism viable in a low energy 
inflationary framework would be the ideal situation but, unfortunately, the simplest 
curvaton model has shown to be incompatible with low enough values for the Hub- 
ble parameter during inflation |126j . Some general proposals to make the curvaton 
paradigm accommodate low scale inflation have recently appeared and specific mod- 
els have been studied too [KT |ll49^ll75j . In the previous chapter, a thermal inflation 
epoch was attached to the general curvaton mechanism making the curvaton field 
gain a huge increment in the mass at the end of the thermal inflationary period, trig- 
gering this way a period of curvaton oscillations, and lowering the main inflationary 
scale to satisfactory levels 'FT]. However, the parameters of the model required for 
this effect to take place showed to be extremely small to affect the reliability of the 
model. The purpose of this chapter is to study the same mechanism but in the case 
where the increment in the mass is so huge that the decay rate becomes bigger than 
the Hubble parameter and the curvaton decays immediately |189j . The results are 
very positive, offering a more natural parameter space. 
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4.2 Thermal inflation and the immediate heavy 
curvaton decay 

The thermal inflation model has been investigated before and found to be a very 
efficient mechanism to dilute the abundance of some unwanted relics, like the moduli 
fields, that the main inflationary epoch is not able to get rid of (see Refs. |136|I137] ). 
We will constrain the available parameter space for A and in the scalar potential 
of Eq. ()3.4H1 so that enough dilution of the moduli abundance is obtained. In 
Chapter IHl this was done for the case in which the fiaton-curvaton coupling term 
gives a huge contribution to the mass of the curvaton when the fiaton field is released 
and gets its vacuum expectation value M. In that case the effective curvaton mass 
TJifj may become bigger than the Hubble parameter giving birth to a period of 
curvaton oscillations and making the scale of the main inflationary period low enough 
{H^ ~ m^j2 ~ 10^ GeV) to think about the inflaton as a field associated to the SUSY 
breaking sector jl26t 1149] I175j . The evolution of the energy densities associated to 
the different fluids in that case are sketched in Fig. 13.11 

The purpose of this chapter is to analyse the scenario where there are no oscil- 
lations of the curvaton field. As it was pointed out in Ref. |175j both the curvaton 
decay rate and the associated lower bound [c.f. Eq. ()3.23j) ] 

r.>4> (4-1) 

TTLp 

are also increased when the fiaton field acquires its vacuum expectation value so that, 
if this increment is big enough for the curvaton decay rate to be bigger than the 
Hubble parameter, the curvaton field may decay immediately rather than oscillating 
for some time. Low scale inflation in this case is also possible to be attained [^I175j . 
but the lower bound on changes with respect to the case when the curvaton 
oscillatory process is triggered. The evolution of the energy densities associated to 
the different fluids in this case are sketched in Fig. 14.11 

In the scenario where curvaton oscillations are allowed, corresponding to Fo- < -ffpt, 
the lower bound on is [c.f. Eqs. fl!0^ . (IT^ . and (HUp ] 

> max{ /^/^ X 10^ GeV, //^^ x 10^^ GeV, 5^ x 10^ GeV } , (4.2) 

where / and 5, given by Eqs. ()3.13|) and ()3.27p . are less than 1, and with i^osc = -f^pt 
being the Hubble parameter at the end of the thermal inflation period (which also 
corresponds to the beginning of the curvaton oscillations). In contrast, the lower 
bound in the scenario where the curvaton field decays immediately, corresponding 
to F,, > i/pt, is [c.f. Eqs. fFO^ and (UHl] 

H, > max{ /^/^ X 10^ GeV, / x 10^^ GeV } , (4.3) 



which may challenge the WMAP constraint for the curvaton scenario j49] leading 
to excessive curvature perturbations from the inflaton field if / is not much smaller 
than unity. 
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A B C log a 

Figure 4.1: Evolution of the energy densities in the thermal inflation model where 
the curvaton field a decays immediately at the end of thermal inflation [189j . The 
continuous line corresponds to the global radiation energy density pr-o? the dashed 
dotted line corresponds to the global fiaton energy density p^^, and the dashed line 
corresponds to the global curvaton energy density po-o- The horizontal axis represents 
the expansion parameter a. From the left to A radiation dominates the energy 
density, although it decreases following p^^ oc a~^. At this stage the unperturbed 
components of the fiaton and curvaton fields x ^i^nd a are frozen at xo = and 
(To = 0"* making their energy densities constants. When p^o reaches Vh at A, thermal 
inflation begins. The thermal inflation period lasts until B when the temperature 
T becomes of the order of the fiaton mass m^. Thermal inflation is portrade by 
the dashed region. After thermal inflation ends, the parametric resonance process 
transforms a substantial fraction of p^ into pr |^ IHHl IHOl \'2-V2\ . The fiaton field is 
liberated by this time and begins oscillating around the minimum of its potential, 
behaving then as a matter fluid with p^^ cx a~^. The curvaton field increments 
suddenly its mass rria at as a result of the oscillations of x around the vacuum 
expectation value M. The increment is enough for the decay rate to overtake ifpt 
(the Hubble parameter at B) so that cr decays immediately. The curvaton energy 
density is transfered then completely to pr as it is the curvature perturbation too. By 
the time C, x already dominates the energy density before decaying into radiation. 
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The lower bound in Eq. ()4.2|) . for H^, ~ 10'^ GeV, was shown in the previous 
chapter to be satisfied for very small values for the flaton-curvaton coupling constant, 
A ~ 10~^^ — 10^^° (see Fig. 13. 2|) . and very small values for the bare mass of the 
curvaton field, mg- < 10~^ GeV [c.f. Eq. fl3.72|) ]. which suggests that the curvaton 
field could be a PNGB finHSOllHEl- This is, in any case, a quite negative result due 
to the required smallness of the parameters A and m^. However, when taking into 
account the lower bound in Eq. ()4.3|) . corresponding to the case when the decay 
rate V^j becomes bigger than ifpt, things change appreciably. 



4.2.1 The flaton-curvaton coupling constant A 

Thermal infiation ends when the thermal energy density is no longer dominant; 
thus, the Hubble parameter at the end of thermal infiation is associated to the 
energy density coming from the curvaton and the fiaton fields: 



3ml 3m| ' ^ ' 



so that 

ffpt~10"^^M. (4.5) 

Since the effective mass of the curvaton field after the end of thermal infiation, i.e., 
when X = average over oscillations of the flaton field and ctq = 0, is 

m„ = {ml + XM'^ f/^ ^ y/XM , (4.6) 

the parameter / [cf. Eq. fl3.13|) ] becomes 

f^E2i^lO-i^J__ (4.7) 

In view of the Eqs. ()4.3j) and ()4.7j) the smallest possible value for A, compatible 
with H^, ~ 10^ GeV, becomes A ~ 10~^°, which is very good because this already 
improves the results found in the previous chapter. Moreover, the effective flaton 
mass during thermal inflation = (m^ — Aa^)^/^ must be positive to trap the 
flaton fleld at the origin of the potential. Thus, Xa^ < m^, and the biggest possible 
value for A becomes [c.f. Eq. fl3.11|l ] 

which is already a small value but much bigger and more natural than that found 
in the case where curvaton oscillations are allowed. The lower bound on A vs ^dec 
is depicted in Fig. 14.21 Note that a small value for fidec, which is restricted to be 
f^dec > 0.01 in order to satisfy the WMAP constraints on non gaussianity |in2|ll38j . 
is desirable to obtain a higher value for A, so the biggest possible value A ~ 10~^ 
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Lower bound on the A coupling constant vs f2dec 



log A 




-1.8 -1.6 -1.4 -1.2 -1 -OS -0.6 -0.4 -0.2 



log Odec 

Figure 4.2: Lower bound on the fiaton-curvaton coupling constant A as a logarithmic 
plot. A more natural value for A requires a higher level of non gaussianity compatible 
with the WMAP constraints. 



is at the expense of a high level of non gaussianity. A smaller upper bound for 
A could be a possibility, according to Eq. ()4.8j) and Fig. 14. 2^ by increasing fldec- 
Nevertheless f^dec? in the scenario studied in this chapter, must satisfy ^Idcc < 1 to 
avoid a period of inflation driven by the curvaton field. 

Recalling, in the scenario where curvaton oscillations are allowed the coupling 
constant A is in the range (see Fig. IH.2j) 

10-22< A<10-i°, (4.9) 

whereas in the scenario where the curvaton decays immediately the range is 

10-i°< A< 10-^ (4.10) 

It is easy to see that the allowed range of values for A in Eq. ()4.1()j) , valid for the case 
where the curvaton field decays immediately at the end of the thermal inflation era, 
is complementary to the allowed range for A in Eq. ()4.9|) . valid when the curvaton 
has some time to oscillate before decaying (see Subsection I3.3.2|) . 

4.2.2 The bare curvaton mass ma 

The only bound on m„ is given by the fact that in the heavy curvaton scenario the 
bare mass must be smaller than the Hubble parameter at the end of the thermal 
inflation era, so that the sudden increment in the mass and the decay rate leads to 
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the immediate decay of the field avoiding in this case the oscillations. Thus, 

</7pt~10-^^M, (4.11) 

so we need to worry about the possible values for M. In the scenario where the 
curvaton field decays immediately the fiaton field is left immersed in a background of 
radiation, so it must decay before the time of nucleosynthesis in order not to disturb 
the abundances of the light elements. By setting ~ if we get the temperature 
just after the fiaton decay 

^ 10" GeV^-^, (4.12) 

which must be bigger than 1 MeV to satisfy the nucleosynthesis constraint. There- 
fore 

M<10^^GeV, (4.13) 

leading to an upper bound on the bare curvaton mass given by mg- < 1 GeV, which is 
again a more relaxed constraint than that found in the previous chapter for the case 
of an oscillating curvaton, but that still reduces the number of possible curvaton 
candidates, leaving essentially the PNGB j^Hl ISHl EEl- Recalling, in the scenario 
where curvaton oscillations are allowed the bare curvaton mass m„ is in the range 
[c.f. Eq. dSH] 

<10-^GeV, (4.14) 

whereas in the scenario where the curvaton decays immediately the range is 

m^<lGeV. (4.15) 

Some important constraints might come from the solution to the moduli problem 
and could limit the reliability of the Eqs. ()4.13j) and ()4.15|1 . Moduli fields are fiaton 
fields with a vacuum expectation value $ of order the Planck mass. The decay 
of the fiaton field increments the entropy density s, so that the big-bang moduli 
abundance, defined as that produced before thermal infiation and given by |137j 

^ ^ (4 16) 

„ in 3/2 1/2 ' \^-^^) 

s lUmp 

where m$ is the mass of the moduli fields, gets suppressed by three factors^ One is 

(4.17) 

due to the curvaton decay, where the are the total internal degrees of freedom, T„ 
is the temperature just after the curvaton decay, and Tc ~ is the temperature 
at the end of thermal infiation; another is 

. 9*{Tpr) Tj,^ 



^Eqs. H4.17|l and (|4.18|) correct a mistake in Ref. |189| . However Eqs. (|4.20l) and 14.2211 are not 

affected by that mistake and, therefore, the conclusions in Ref. |189| about the lower and upper 
bounds on M remain unchanged. 
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due to the parametric resonance process [231 UHl I212j following the end of the 
thermal inflation era, where Tpji is the temperature just after the period of preheat- 
ing; and the other is 

A (4 19) 

due to the flaton decay, where is the temperature just after the decay^, and P is 
the fraction of the total energy density left in the fiatons by the parametric resonance 
process (/3 < 1). Thus, the abundance of the big-bang moduli after thermal inflation 
is: 



> 10^ GeV^M" 



r\-/ 




which must be suppressed enough {n<^/s < 10^^^) so that the nucleosynthesis con- 
straints studied in Ref. jSSl are satisfied. This is easily achieved by imposing a lower 
bound on M: 

M>109GeV, (4.21) 

which does not affect the upper bounds on M and mg. in Eqs. ()4.13p and ()4.15p . 

We also have to take care about the abundance of the thermal inflation moduli, 
defined as that produced during the preheating stage following the end of the thermal 
inflation era: 



X 



> lO-^'GeWf^^Vf-^l X 

\mpj VlMeV 



p 




(4.22) 

Here $r corresponds to the vacuum expectation value of the thermal moduli fields. 
To suppress the thermal inflation moduli at the required level nq,^/s < 10~^^ we 
require 

M<10i^GeV, (4.23) 

which is precisely the same bound as in Eq. ()4.13|) . Recalling, the allowed range of 
values for the vacuum expectation value of the flaton fleld is 

10^ GeV <M< 10^*^ GeV , (4.24) 



^This is assuming for simplicity that the flaton has come to dominate the energy density just 
before decaying (see Fig. I4.f |l . 
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so that the moduh problem is solved and, in the best case, mo- ~ 1 GeV. The latter 
allowed range for M means that, unlike the case where the curvaton has some time 
to oscillate before decaying, the flaton field could be the GUT Higgs field studied in 
Ref. 

4.3 Some useful remarks 

Before concluding, we want to stress some points that can help to avoid possible 
confusion. The parameter space compatible with low scale inflation is a feature of 
the specific model studied, and we cannot say it is the same for all classes of models 
in the basis of Eqs. ()3.22|) . ()3.26p . and ()3.30|) . which provide just some general 
bounds. That is why specific models have been studied (see Refs. |^ I175j ). even 
when the general bounds were already known from Refs. \V2{S\ I149j . Although the 
claim, that the available parameter space is bigger for the immediate curvaton decay, 
was given before in Ref. |175j . we again cannot say that the available parameter 
space is the same for all classes of models in the basis of the bounds required to 
have low energy scale inflation. For example, from Eq. ()4.5|) . iJpt depends on M so 
there is no direct bound on it unless we know the bound on M ^. The bound on M 
comes in turn from the requirement that the fiaton decays before nucleosynthesis 
[c.f. Eqs. fl4.12|l and ()4.13|l ] and must be consistent with the adequate suppression 
of the thermal inflation moduli [c.f. Eqs. ()4.22j) and ()4.23p ]. These are, of course, 
features specific only to the model we are studying, and are therefore not present in 
Ref. HZSj. 

Naively, one would think that the bounds on A and mo- are found from that on 
-ffpt only through a mere change of variables. This is of course not true as the bound 
on ifpt is a very sensitive quantity that has to avoid disturbing the nucleosynthesis 
process and the adequate moduli abundance suppression. It is worth mentioning 
that the scenario discussed in this chapter differs appreciably from that studied in 
Chapter El due to the immediate curvaton decay, so that the conditions to satisfy the 
nucleosynthesis and thermal inflation moduli constraints are completely different^. 

Finally, the agreement between the bounds found in Ref. |175j (which are sup- 
posed to be general) and those found in this chapter is apparent and corresponds 
just to a mere coincidence. We justify this observation by noting that Eq. (6) in 
Ref. |175j is essentially the same as our Eq. ()3.1()j) . the latter being generalized to 
give Eq. ()3.22|1 . except for which in our Eq. ()3.1(ij) appears to be if dec- The 
expressions in the previous chapter were carefully derived so that the correct ex- 
pression is that given there [5T] . In contrast, Eq. (6) in Ref. |175j is just valid 
for the standard case where the curvaton field has some time to oscillate before 

^Notice that the bounds required to have low energy scale inflation [c.f. Eq. 14.3|l ] depend only 
on the ratio / = Hpt/rha-, and not exclusively on Hpt. 

^For example, the expressions for the big-bang and thermal inflation moduli abundances in 
Chapter|31[c.f. Eqs. H3.78(l and H3.80|l ] are different from those in this chapter [c.f. Eqs. H4.20(l and 
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decaying, so we can identify To- with ifdec- However, for the immediate decay case. 
To- > iJdec = -f^pt, which renders Eq. (6) in Ref. |175j invahd. Based on the previous 
discussion we claim that the bound H^i < 1 GeV, as are those on A and m^, is 
presented in this thesis for the first time in a correct way. 

4.4 Conclusions 

In this chapter we have investigated the required parameter space compatible with 
low scale inflation in the thermal inflation curvaton scenario where there are no 
oscillations of the curvaton field |189j . We have shown that the parameter space 
is greatly enhanced when the increment in the curvaton decay rate is big enough 
for the curvaton field to decay immediately at the end of the thermal inflation era. 
The best case corresponds to a flaton-curvaton coupling constant A ~ 10~^ and a 
bare curvaton mass mg- ~ 1 GeV, which are much bigger and more natural than 
the ranges 10^^^ ^ ^ 10^^° and mg- < 10^^ GeV found previously in Chapter IHl 
for the case where the curvaton oscillates for some time before decaying j3T]. In 
addition we have found lO^GeV < M < lO^^GeV for the vacuum expectation value 
M of the flaton field. Therefore, our flaton field as the GUT Higgs field discussed 
in Ref. jTB6ij is a viable option in this scenario. 
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Chapter 5 



Non-gaussianity from the 
second-order cosmological 
perturbation 

5.1 Introduction 

In chapters 121 El andlU we discussed some of the theoretical aspects of the origin of 
the large-scale structure in the Universe, emphasising the possibility to achieve low 
scale inflation in the curvaton scenario. We begin now the discussion of the statistical 
aspects, specifically the presence of non-gaussianities in the fields responsible for the 
origin of the curvature perturbation ( and/or in ( itself. Chapters El and IHl will deal 
with such an interesting subject. 

Cosmological scales leave the horizon during inflation and re-enter it after Big 
Bang Nucleosynthesis. Throughout the super-horizon era it is very useful to define 
a primordial cosmological curvature perturbation, which is conserved if and only 
if pressure throughout the Universe is a unique function of energy density (the 
adiabatic pressure condition) (see Subsection ITT^ [HI 1121 1121 II2I1 CZni 1^ • 
Observation directly constrains the curvature perturbation at the very end of the 
super-horizon era, a few Hubble times before cosmological scales start to enter the 
horizon, when it apparently sets the initial condition for the subsequent evolution of 
all cosmological perturbations. As discussed in Chapter |21 the observed curvature 
perturbation is almost Gaussian with an almost scale-invariant spectrum. 

Cosmological perturbation theory expands the exact equations in powers of the 
perturbations and keeps terms only up to the nth order. Since the observed cur- 
vature perturbation is of order 10~^, one might think that first-order perturbation 
theory will be adequate for all comparisons with observation. That may not be the 
case however, because the PLANCK satellite |171j and its successors may be sen- 
sitive to non-gaussianity of the curvature perturbation at the level of second-order 
perturbation theory |il03j . 

Several authors have treated the non-gaussianity of the primordial curvature 



75 



perturbation in the context of second-order perturbation theory. They have adopted 
different definitions of the curvature perturbation and obtained results for a variety 
of situations. In this chapter we revisit the calculations, using a single definition of 
the curvature perturbation which we denote by ( [1'31]. In some cases we disagree 
with the findings of the original authors. 

The outline of this chapter is the following: in Section we review two defini- 
tions of the curvature perturbation found in the literature, which are valid during 
and after inflation, and establish definite relationships between them; in section ESI 
we review a third curvature perturbation definition, which applies only during in- 
flation, and study it in models of inflation of the slow-roll variety; in Section E31 we 
describe the present framework for thinking about the origin and evolution of the 
curvature perturbation; in Section 15.51 we see how non-gaussianity is defined and 
constrained by observation; in Section 15.61 we study the initial non-gaussianity of 
the curvature perturbation, a few Hubble times after horizon exit; in Section IKTI we 
study its subsequent evolution according to some different models. The conclusions 
are summarised in Section ESI 

We shall denote unperturbed quantities by a subscript 0, and generally work with 
conformal time t]. Sometimes though we revert to physical time t. We shall adopt 
the convention that a generic perturbation g is split into a first- and second-order 
part according to the formula 

^ = ^1+2^2- (5.1) 



5.2 Two definitions of the curvature perturbation 
5.2.1 Preliminaries 

Cosmological perturbations describe small departures of the actual Universe, away 
from some perfect homogeneous and isotropic universe with the line element in 
Eq. ()2.2|) . For a generic perturbation it is convenient to make the Fourier expansion 

^7(x, V) = J d'k ^7k(r/)e^'^-^ , (5.2) 

where the spacetime coordinates are those of the unperturbed Universe. The inverse 
of the comoving wavenumber, , is often referred to as the scale. Except where oth- 
erwise stated, our discussion applies only to the super-horizon regime {k <^ aHi^f). 

When evaluating an observable quantity only a limited range of scales will be 
involved. The largest scale, relevant for the low multipoles of the Cosmic Microwave 
Background anisotropy, is k~^ ~ -f^today where -fftoday is the present Hubble param- 
eter. The smallest scale usually considered is the one enclosing matter with mass 
~ IO^Mq, which corresponds to k~^ ~ 10^^ Mpc ~ 10^^H^^\^^. The cosmological 
range of scales therefore extends over only six orders of magnitude or so. 

To define cosmological perturbations in general, one has to introduce in the 
perturbed Universe a coordinate system (t,x*), which defines a slicing of spacetime 



76 



(fixed t) and a threading (fixed x'). To define the curvature perturbation it is enough 
to define the slicing [127\ . 

5.2.2 Two definitions of the curvature perturbation 

In this chapter we take as our definition of ( the following expression for the spatial 
metric |SHl 1131131121111121111111^11^111112001 which applies non-perturbatively: 

Qij = a^{vi)^ije^^ . (5.3) 

Here 7jj has unit determinant, and the time-slicing is one of uniform energy density^. 

It has been shown under weak assumptions |127j that this defines C, uniquely, 
and that C, is conserved as long as the pressure is a unique function of energy density. 
Also, it has been shown that the uniform density slicing practically coincides with 
the comoving slicing (orthogonal to the flow of energy), and with the uniform Hubble 
slicing (corresponding to uniform proper expansion, that expansion being practically 
independent of the threading which defines it) |127j . The coincidence of these slicings 
is important since all three have been invoked by different authors. 

Since the matrix 7 has unit determinant it can be written 7 = Je'*, where / is 
the unit matrix and h is traceless jl27j . Assuming that the initial condition is set 
by inflation, h corresponds to a tensor perturbation (gravitational wave amplitude) 
which will be negligible unless the scale of inflation is very high. As we shall see later 
(see footnote^! in this chapter), the results we are going to present are valid even 
if h is not negligible, but to simplify the presentation we drop h from the equations. 
Accordingly, the space part of the metric in the super-horizon regime is supposed 
to be well approximated by 

g,j = a^{7])5ije^^ . (5.4) 
At first order, Eq. ()5.4|) corresponds to 

g,j = a\Ti)6,,{l + 2C). (5.5) 

Up to a sign, this is the definition of the first-order curvature perturbation adopted 
by all authors [c.f. Eqs. ()2.3|) and ()2.15|) ]. There is no universally agreed convention 
for the sign of C,. Ours coincides with the convention of most of the papers to which 
we refer, and we have checked carefully that the signs in our own set of equations 
are correct. 

At second order we have 

g,,=a\r^)6,,{l + 2C + 2e). (5.6) 

This is our definition of Q at second order |131j . 

^It is proved in Ref. |127j that this definition of C, coincides with that of Lyth and Wands 
|140j . provided that their shces of uniform coordinate expansion are taken to correspond to those 
on which the Hne element has the form Eq. H5.3(l without the factor e^^ (this makes the shces 
practicaUy flat if 7.y ~ 5ij). 
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Malik and Wands |145j instead defined ( by Eq. ()5.5p even at second order. 
Denoting their definition by a subscript MW, 



C^^ = C + C', (5.7) 

or equivalently 

Cs^w = C2 + 2 (Ci)' , (5.8) 

where Ci is the first-order quantity whose definition Eq. ()5.5|) is agreed by all authors. 

To make contact with calculations of the curvature perturbation during inflation, 
we need some gauge-invariant expressions for the curvature perturbation. As stated 
in Subsections I2I221 and I2I2IS1 'gauge- invariant' means that the definition is valid 
for any choice of the coordinate system which defines the slicing and threading^. 

We shall write gauge- invariant expressions in terms of ( and (^^w. First we 
consider a quantity ■i/'^^, defined even at second order by 

^,, = a2(r^)5,,(l-2V^^W). (5.9) 

This definition, which is written in analogy to Eq. ()5.5|1 . applies to a generic slicing. 
Analogously to Eq. ()5.4p we can consider a quantity ip, valid also in a generic slicing, 
defined by 

g,, = a\ri)S,je-^^ . (5.10) 

On uniform- density slices, ipi = ip^^ = —(i, ip^^ = — Cl^^, and 1^2 = —(2- We 
shall also need the energy density perturbation 6p, defined on the generic slicing, as 
well as the unperturbed energy density po- 

At first order, the gauge-invariant expression for ( has the well-known form 
[c.f. Eq. (ITTHll ] 

Ci = -^i-H^, (5.11) 
Po 

where TC = a' /a, and the unperturbed energy density satisfies pg = — 37i(po + -Po) 
with Pq being the unperturbed pressure. This expression obviously is correct for the 
uniform density slicing, and it is correct for all slicings because the changes in the 
first and second terms induced by a change in the slicing cancel [TTJ [T21 EH I124t I147j . 
At second order, Malik and Wands show that jl45j 

^Mw ^ -^Mw_^^^27^^^ + 2^(^^ + 27^^0 

Po Po Po Po 

which is, again and for the same reason as before, obviously correct for all the 
slices. Accordingly, from Eq. ()5.8j) . we can write a gauge invariant definition for our 

^In the unperturbed limit the slicing has to be the one on which all quantities are uniform and 
the the threading has to be orthogonal to it. 
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second-order (: ^ 



Po Po Po Po \ Po J \Wo 

where the relation 

^|;^^ = ^,-2M\ (5.14) 
coming from Eqs. and ()5.1()|1 . has been used. 

5.3 Slow-roll inflation and a third deflnition 

Now we speciahze to the era of slow-roll inflation |Hl I112[ I117| llHOj . We con- 
sider single-component inflation, during which the curvature perturbation ( is con- 
served, and mult i- component inflation during which it varies. After defining both 
paradigms, we give a third definition of the curvature perturbation which applies 
only during inflation. 



5.3.1 Single-component inflation 

In a single-component inflation model jll2| ll3Uj the inflaton trajectory is by defi- 
nition essentially unique. The inflaton field if parameterises the distance along the 
inflaton trajectory. In terms of the field variation, slow-roll inflation (see Subsection 
12.5.11) is characterised by the slow-roll conditions j lllf I112[ I13()j 



inf 



.-"inf. 



< 1 



Vo 



-f^infV^O 



< 1. 



(5.15) 
(5.16) 



The inflaton field can be taken to be canonically normalised, in which case these 
definitions are equivalent to conditions on the potential V 



dV 



2^2 ydifio, 

nip d'^V 



(5.17) 
(5.18) 



which, together with the slow-roll approximation, lead to the slow-roll behaviour 

dV 



d^po 



(5.19) 



■^This relation has recently been confirmed in Ref. jl()fi| (see also Ref. |1()7| ) using a nonlinear 
coordinate-free approach. 
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Even without the slow-roll approximation, slices of uniform Lp correspond to 
comoving slices because a spatial gradient of (y? would give non-vanishing momentum 
density. Since comoving slices coincide with slices of uniform energy density, the 
slices of uniform Lp coincide also with the latter. Also, since is a Lorentz scalar, its 
gauge transformation is the same as that of p. It follows [2175 ^^^^ replace 
p by <y9 in the above expressions: 



Ci 



S2 



inf 



..MW 



inf 



5^ 












inf , + 27i 


■)■ 


[ < 












■)■ 


( 







(5.20) 
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(5.21) 




(5.22) 



5.3.2 Multi-component inflation 

Now consider the case of multi-component inflation, where there is a family of in- 
equivalent inflationary trajectories lying in an A^-dimensional manifold of field space. 
If the relevant part of the manifold is not too big it will be a good approximation to 
take the fields to be canonically normalised. Then the inequivalent trajectories will 
be curved in field space"^. To define the trajectories one can choose a fixed basis in 
field space corresponding to fields 0i, ■ ■ ■ , 07v 

Assuming canonical normalisation, multi-component slow-roll inflation is char- 
acterised by the conditions 



V 



( dV ' 

\d(t>no, 



3-f^inf0no 



dV 



(5.23) 
(5.24) 
(5.25) 



The procedure of choosing a fixed basis is quite convenient for calculations, but 
a different procedure leads to a perhaps simpler theoretical description. This is 
to take if to parameterise the distance along the infiaton trajectories, just as in 
single-component inflation, but now with the proviso that uniform if corresponds 

■*More generally they will be non-geodesics, the geodesies being the trajectories which the 
background fields could follow if there was no potential term in the scalar Lagrangian |17S | . 



80 




V = const 

V = const 
V = const 



Figure 5.1: Two different procedures for defining the fields in two-component infla- 
tion. The fields are denoted by ip and '(9. (a) The field (p parameterises the distance 
along the inflaton trajectories, with uniform ip corresponding to the equipotential 
lines. The field parameterises the distance along the equipotentials. (b) The 
fields ip and i} are the components in a fixed orthonormal basis, aligned with the 
inflationary trajectory at a certain point in field space. The value of ip is now 
the displacement along the tangent vector and the value of ^9 is the displacement 
along the orthogonal vector. Working to second order in these displacements, the 
equipotentials no longer coincide with the lines of uniform ip. 



to uniform field potential (since we work in the slow-roll approximation, this means 
that the slices in field space of uniform p are orthogonal to the trajectories). Then, 
in the slow-roll approximation, slices of spacetime with uniform ip will again coincide 
with slices of uniform density (see Fig. 15.1b .). Since is a scalar, Eqs. ()5.2()j] and 
()5.21|) will then be valid. This is the simplest form of the gauge-invariant expression, 
though for a practical calculation it may be better to write it in terms of a fixed 
basis. 

There is a subtlety here. For the first-order case we could define p in a different 
way; around a given point on the unperturbed trajectory we could choose a fixed 
field basis, with one of the basis vectors pointing along the trajectory, and define p 
as the corresponding field component. Then we could choose p to be canonically 
normalised in the vicinity of the chosen point in field space. That would not work 
at second order though, because at that order it makes a difference whether p is 
the appropriate parameterisation of the distance along the trajectories (our adopted 
definition) or the distance along a tangent vector to the trajectory (the alternative 
definition) (see Fig. 15.1b ). Only our adopted one will make Eqs. fl5.21|) and ()5.22|) 
valid. 
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5.3.3 A third definition of the curvature perturbation 



The third definition in the hterature apphes only during infiation. It was given orig- 
inally by Acquaviva et. al. |3] for the single- component case, and the generalization 
to the multi-component case was noted by Rigopoulos |178j . We shall denote this 
definition by C^. 

The definition of Acquaviva et. al. and Rigopoulos is 

C2 = -W2 - -^inf—^ nj, ,9-1/2 nj ,,,,, ■ (5.26) 

V^O ^inf + ^^inf - riiniVolVo 

This is gauge-invariant by construction, with Lp defined as in Figure 1(a). 

It was pointed out by Vernizzi |217j (actually in the context of single-component 
inflation) that comparing this definition with Eq. ()5.2H1 gives simply 

_ /-MW ^^fniiClY /r oyN 

^2 ~ ^2 nj, I 9-172 -17 ,nl>l,r,l ' {O.ZI) 

^inf + ^^inf - '^infV^o/V'o 

In the limit of slow-roll the denominator of the last term becomes just 27if^f, and 
then 

= C2. (5.28) 

In other words, this third definition coincides with our adopted one in the slow-roll 
limit. 

Making use of the slow- roll parameters in Eqs. ()5.15j) and 1)5.1(11) . the expression 
in Eq. ()5.27|) gives to first-order in the slow-roll approximation 

C2^ = C2-(2£-r/^)(Ci)'. (5.29) 



5.4 The evolution of the curvature perturbation 

The simplest possibility for the origin of the observed curvature perturbation is that 
it comes from the vacuum fiuctuation of the infiaton field in a single-component 
model (see Section f2.5|) . More recently other possibilities were recognised and we 
summarise the situation now. Although the discussion is usually applied to the 
magnitude of the curvature perturbation, it applies equally to the non-gaussianity. 



5.4.1 Heavy, fight and ultra-fight fields 

On each scale the initial epoch, as far as classical perturbations are concerned, should 
be taken to be a few Hubble times after horizon exit during inflation. The reason is 
that all such perturbations are supposed to originate from the vacuum fluctuation 
of one or more light scalar fields, the fiuctuation on each scale being promoted to a 
classical perturbation around the time of horizon exit jH EHl 1^ I122| I124j . 
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Considering a fixed basis with canonical normalisation, a light field is roughly 
speaking one satisfying the flatness condition in Eq. ()5.24|) . The terminology is 
suggested by the important special case that the effective potential during inflation 
is quadratic. Then, a light fleld is roughly speaking that whose effective mass during 
inflation is less than the value if* of the Hubble parameter. More precisely, the 
condition that the vacuum fluctuation be promoted to a classical perturbation is 

m < ^H^ . (5.30) 

From now on we focus on the quadratic potential, and take this as the the deflnition 
of a light fleld. Conversely a heavy fleld may be deflned as one for which the condition 
in Eq. ()5.3m) is violated. 

During inflation light flelds slowly roll according to Eq. ()5.25|) (with the vacuum 
fluctuation superimposed) while the heavy flelds presumably are pinned down at an 
instantaneous minimum of the effective potential. As we have seen, multi-component 
inflation takes place in a subspace of fleld space. The flelds in this subspace are 
light, but their effective masses are sufficient to appreciably curve the infiationary 
trajectories. In the case of both mult i- component and single- component inflation, 
there could also be 'ultra-light' flelds, which do not appreciably curve the inflationary 
trajectory and which therefore have practically no effect on the dynamics of inflation. 

5.4.2 The evolution of the curvature perturbation 

To describe the behaviour of perturbations during the super-horizon era, without 
making too many detailed assumptions, it is useful to invoke the separate universe 
hypothesis jT21 \V27\ I14()| I197| 1222 ] after smoothing on a given comoving scale much 
bigger than the horizon^. According to this hypothesis the local evolution at each 
position is that of some unperturbed universe (separate universe). Of course the 
separate universe hypothesis can and should be checked where there is a sufficiently 
detailed model. However, it should be correct on cosmological scales for a very 
simple reason. The unperturbed Universe may be defined as the one around us, 
smoothed on a scale a bit bigger than the present Hubble distance. In other words, 
the separate universe hypothesis is certainly valid when applied to that scale. But 
the whole range of cosmological scales spans only a few orders of magnitude. This 
means that cosmological scales are likely to be huge compared with any scale that is 
relevant in the early Universe, and accordingly that the separate universe hypothesis 
should be valid when applied to cosmological scales even though it might fail on much 
smaller scales (this expectation was verified in a preheating example |113j to which 
we return later). 

^When considering linear equations, smoothing is equivalent to dropping short wavelengths 
fourier components. In the nonlinear situation the smoothing procedure could be in principle 
ambiguous. In a given situation one should state explicitly which quantities are being smoothed. 
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We are concerned with the curvature perturbation, which during the super- 
horizon era is conserved as long as the pressure is a unique function of the energy 
density (the adiabatic pressure condition) (see Subsection I2.2.3|) . The adiabatic 
pressure condition will be satisfied if and only if the separate universes are identical 
(at least as far as the relation between pressure and energy density is concerned) ^. 
The condition to have identical universes after a given epoch is that the specification 
of a single quantity at that epoch is sufficient to determine the entire subsequent 
evolution. 

In the case of single-component inflation, the initial condition may be supplied 
by the local value of the infiaton field, at the very beginning of the super-horizon 
era when it first becomes classical. Given the separate universe hypothesis, that is 
the only possibility if the infiaton is the only light field ever to play a significant 
dynamical role. This means that the curvature perturbation generated at horizon exit 
during single- component inflation will be equal to the one observed at the approach 
of horizon entry, provided that the inflaton is the only light field ever to play a 
dynamical role. 

If infiation is mult i- component, more than one field is by definition relevant 
during inflation. Then the curvature perturbation cannot be conserved during infla- 
tion. The variation of the curvature perturbation during multi-component inflation 
is caused by the vacuum fluctuation orthogonal to the unperturbed inflationary tra- 
jectory, which around the time of horizon exit kicks the trajectory onto a nearby 
one so that the local trajectory becomes position-dependent. After inflation is over, 
the curvature perturbation will be conserved if the local trajectories lead to prac- 
tically identical universes. In other words it will be conserved if the light (and 
ultra-light) flelds, orthogonal to the trajectory at the end of inflation, do not affect 
the subsequent evolution of the Universe. 

The curvature perturbation after inflation will vary if some light or ultra-light 
fleld, orthogonal to the trajectory at the end of inflation, affects the subsequent 
evolution of the Universe (to be precise, affects the pressure) [SHI EI!- As we shall 
describe in Section 15. 7| three types of scenario have been proposed for this post- 
inflationary variation of the curvature perturbation. 

5.5 Non-gaussianity 

5.5.1 Defining the non-gaussianity 

A gaussian perturbation is one whose Fourier components are uncorrelated . All 
of its statistical properties are deflned by its spectrum, and the spectrum Vg{k) = 

^Of course the identity will only hold after making an appropriate synchronization of the clocks 
at different positions. Having made that synchronization, horizon entry will occur at different 
times in different positions, which can be regarded as the origin of the curvature perturbation. 
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A^lk / aHi^^{)^s of generic perturbation is conveniently defined jll2[ I13U] by^ 

27r2 



ciS'ka) 



-6'ik,+k2)Vgik) 



the normalisation being chosen so that 



(/(x))^/;p.(.)f. 



(5.31) 



(5.32) 



On cosmological scales a few Hubble times before horizon entry, observation shows 
that the curvature perturbation is almost Gaussian with |y4^| ^ 5 x 10~^ [20fj. 

The simplest kind of non-gaussianity that the curvature perturbation could pos- 
sess is of the form 



c(x) = Cg(x) - |/nl (c|(x) - (CD) 



(5.33) 



where Q is Gaussian with {Q) = 0, and the non-linearity parameter /nl is inde- 
pendent of position. We will call this correlated non-gaussianity. Note that this 
definition assumes that {() = 0, which means that the zero Fourier mode (spatial 
average) is dropped. 

Following Maldacena |142j . we have inserted the prefactor —(3/5) so that in 
first-order perturbation theory our definition agrees with that of Komatsu and 
Spergel [103', which is generally the definition people use when comparing theory 
with observation. Working in first-order perturbation theory, these authors write 
$(x) = $3(x) + /nl ('^g(x) - i^l)), and their $ is equal to -3/5 times our ( ^. 

One of the most powerful observational signatures of non-gaussianity is a nonzero 
value for the three-point correlator, specified by the bispectrum B defined by 

(5.35) 



(Ck,Ck,Ck3) = {2n)-'/'Bik,, k2, ks)S%k^ + k2 + ks) . 
For correlated non-gaussianity (with the gaussian term dominating) 



B{ki, k2, fcs 



6. 

-g/NL 



Pc_{ki)P(^{k2) + cyclic permutations 



(5.36) 



^Technically the expectation values in this and the following expressions refer to an ensemble of 
universes but, because the stochastic properties of the perturbations are supposed to be invariant 
under translations, the expectation values can also be regarded as averages over the location of the 
observer who defines the origin of coordinates. 

®The actual quantity constrained by observational data is /^lj which is the non-linearity pa- 
rameter for the CMB temperature anisotropics: 



5r 

T 



5T 
T 



^(x)=(:^) (x) + /^L 



(5.34) 



At first order = 3/nl because 5T/Tq = (-1/5)C [c.f. Eqs. ITTJ and ^^]. However, to 
compare adequately the observational data with our /nl, we must calculate /^l in terms of /nl 
at second order (see e.g. Refs. [El UHl [131211 US)- 
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where P({k) = 27r^'P^(A;) / k^. For any kind of non-Gaussianity one may use the above 
expression to define a function /nl(^15 ^2, ^s)- 

Given a calculation of /nl using first-order perturbation theory, one expects in 
general that going to second order will change /nl by an amount of order 1 . On this 
basis, one expects that a first-order calculation is good enough if it yields |/nl| S> 1, 
but that otherwise a second-order calculation will be necessary. 

The definition Eq. ()5.36|) of /nl is made using our adopted definition of (. If ( 
in the definition is replaced by ^'^^ (with the zero Fourier mode dropped) then /nl 
should be replaced by 

/iT - /nl - ^ • (5.37) 

To obtain this expression we used Eq. ()5.7|1 and dropped terms higher than second 
order^. 

All of this assumes that the non-gaussian component of ( is fully correlated with 
the gaussian component. An alternative possibility fS^ that will be important for 
us is if ( has the form 

C(x) = C,,(x)-^/nl(CJ(x)-(CJ)) , (5.38) 

where (g and are uncorrected Gaussian perturbations, normalised to have equal 
spectra, and the parameter /nl is independent of position. We will call this uncor- 
related non-gaussianity. It can be shown [22^ that in this case, /nl as defined by 
Eq. ()5.36j) is given by 

5.5.2 Observational constraints on the non-gaussianity 

Taking /nl to denote the non-linearity parameter at the primordial era, let us con- 
sider the observational constraints. Detailed calculations have so far been made 
only with /nl independent of the wavenumbers, and only by using first-order per- 
turbation theory for the evolution of the cosmological perturbations after horizon 
entry. It is found |l()2j that present observation requires |/nl| ^ 10^ making the 
non-gaussian fraction at most of order 10~^. The use of first-order perturbation 
theory in this context is amply justified. Looking to the future though, it is found 
that the PLANCK satellite will either detect non-gaussianity or reduce the bound 
to |/nl| ^ 5 [m ll()3j . and that foreseeable future observations can reach a level 

i/nl| ~3 [niiiiis]. 

Although the use of first-order perturbation theory is not really justified for the 
latter estimates, we can safely conclude that it will be difficult for observation ever 
to detect a value |/nl| "C 1. That is a pity because, as we shall see, such a value is 

^Obviously the parameter /^l: which is the important one to make comparison with observa- 
tional data, does not depend on the chosen definition for /nl- 
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predicted by some theoretical scenarios. On the other hand, other scenarios predict 
|/nl| roughly of order 1. It will therefore be of great interest to have detailed second- 
order calculations, to establish precisely the level of sensitivity that can be achieved 
by future observations. A step in this direction has been taken in Ref. [20] (see 
also Refs. fH ^ I18 | 1^)- where a non-linear expression for the large-scale CMB 
anisotropy is given in terms of only the curvature perturbation (generalizing the 
first-order Sachs- Wolfe effect |191j ). 



5.6 The initial non-gaussianity 
5.6.1 Single-component inflation 

At first order, the curvature perturbation during single-component infiation is Gaus- 
sian. The amplitude of its time-independent spectrum is given by |112[ lllHOj [c.f. 

Eq, mm] 

and its spectral index = dlnV(^{k)/dlnk is given by 

= 2r]^ - 6e . (5.41) 

The squared amplitude of the spectrum r^t^ of the tensor perturbation, defined as a 
fraction of A^, is also given in terms of the slow-roll parameter e [c.f. Eq. ()2.127|) ]: 

ttc = I6e . (5.42) 

// the curvature perturbation does not evolve after single-component inflation is 
over observation constrains and r-r^, and hence the slow-roll parameters t]^ and 
e. A current bound |2()7j is —0.048 < < 0.016 and vtq < 0.46. The second bound 
gives e < 0.029, but barring an accurate cancellation the first bound gives e < 0.003. 
In most infiation models e is completely negligible and then the first bound gives 
—0.024 < Tj^p < 0.008 (irrespective of slow-roll infiation models, the upper bound in 
this expression holds generally, and the lower bound is badly violated only if there 
is an accurate cancellation). The bottom line of all this is that e and \r]^\ are both 
constrained to be < 10^^. 

Going to second order, Maldacena |142j has calculated the bispectrum during 
single-component infiation (see also Refs. 03 iSl IHS El dll IMl W^) His 

^"in Ref. |181j (see also Ref. J180') Rigopoulos et. al. calculated the three-point correlator in 
single-component slow-roll inflation using a stochastic approach. Their result agrees with Malda- 
cena's in the squeezed limit (where one of the scales crosses the horizon much earlier than 
the other two, ki <C A:2, fcs), but disagrees in the limit where the fc^'s form an equilateral triangle. 
Calcagni in Ref. extended this stochastic approach to calculate the non-gaussianity originated 
from a Dirac-Born-Infeld tachyonic inflaton and in braneworld scenarios, finding results identical 
to Maldacena's one. The three-point correlators calculated for both cases were found identical to 
that calculated in Ref. |181| . 
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result may be written in the form 



/nl = - [2r]^ -6e- 2ef{h, fcs, ^3)] , (5.43) 

with < / < 5/6. By virtue of the slow-roll conditions, |/nl| <^ 1 In other 
words, the curvature perturbation (, as we have defined it, is almost Gaussian during 
single-component inflation. 

From Eq. ()5.29|1 is also practically gaussian, but this quantity is defined 
only during inflation and therefore could not be considered as a replacement for (. 
More importantly, has significant non-gaussianity because, from Eq. ()5.37|) . it 
corresponds to ~ —5/3. 

One may ask why it is our ( and not (^"^ which is gaussian in the slow-roll limit. 
One feature that distinguishes our (, is that any part of it can be absorbed into the 
scale factor without altering the rest; indeed 

g,, = 6,,a\Ti)e^^^^^' = 6,fa\Ti)e^^ , (5.44) 

with a = ae^i (if we tried to do that with Cmw, the part of ( not absorbed would have 
to be re-scaled). This means that an extremely long- wavelength and possible large 
part of ( has no local significance. It also means, in the context of perturbation 
theory, that the first-order part of ( can be absorbed into the scale factor when 
discussing the second-order part. However, the gaussianity of ( does not seem to be 
related directly to this feature. Rather, it has to do with the gauge transformation, 
relating quantities ipA and V'b defined on different slicings. 
With our definition jl27j . the gauge transformation is 

x) - ^t, x) = -ANABit, x), (5.45) 

where ANab is the number of e-folds of expansion going from a slice -B to a slice A, 
both of them corresponding to time t In writing this expression we used physical 
time t instead of conformal time, the two related by dt = adrj. Along a comoving 
worldline, the number of e-folds of expansion is defined as N = J Hdr where H is 
the local Hubble parameter and dr is the proper time interval |127j . 

To understand the relevance of this result, take ■05 = and iI)a = —(■ The 
pressure is adiabatic during single-component inflation, which means that dt can be 

-'^^Near a maximum of the potential 'fast-roll' inflation |31II115| can take place with \ri^ \ somewhat 
bigger than 1. Maldacena's calculation does not apply to that case but, presumably, it gives initial 
non-gaussianity |/nl| ^ 1- Although the corresponding initial spectral index is far from 1, which 
means that the initial curvature perturbation produced by ip must be negligible, the precise initial 
value of /nl may in this case be important as long as another field (like the curvaton) be in charge 
of generating the observed curvature perturbation. 

^^This expression is valid even when the tensor perturbation is included jl27| . As a result, the 
gauge-invariant expressions mentioned earlier are still valid in that case, as are the results based 
on them including the present discussion. 
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identified with the proper time interval dr, and the proper expansion rate on shcing 
A is uniform |127j . As a result, to second order, 

C = i/i„f(t)At(t,x) + ^iji,f(t)(At(t,x))' 

~ i/i„fAt(t,x) + i|^(i7i,fAt(t,x))=^ 

^ ifi„fAt(t,x). (5.46) 

In the last line we made the slow-roll approximation, and from the second line we 
can see that the error in /nl caused by this approximation is precisely e. 

We also need the gauge transformation for the inflaton field (p in terms of At. 
Since the slices correspond to the same coordinate time, the unperturbed inflaton 
field can be taken to be the same on each of them which means that the gauge 
transformation for 6 if is 

5LpA{t, x) - 5(fB{t, x) = AipAB{t, x), (5.47) 

where Aip^B is the change in ip going from slice B to slice A. But slice A corresponds 
to uniform ip, which means that on slice B to second order 

i/inf(t)^^^^ = -i/inf(t)At(t,x)-^i7i„f(i)-(Ai(t,x))2 

ipo 2 ipo 

- -iJinfAt(t,x)-i-^(i7i,fAt(t,x))2 

- -/JinfAt(t,x), (5.48) 

where in the last line we used the slow-roll approximation. We can see that the 
fractional error caused by this approximation is v3o/-ffinf0o = e — rj^,. 

Combining Eqs. ()5.4fi|l and ()5.48j) we have in the slow-roll approximation 

C^-H„„W^^^. (5.49) 

with fractional error of order maxjr^,^, e} (this can also be seen directly from Eqs. ()5.2()j) 
and ()5.22|1 evaluated with ip = 0, but we give the above argument because it explains 
why the result is valid for ( as opposed to (^^^y 

The final and crucial step is to observe that in the slow-roll approximation is 
gaussian, with again a fractional error of order max{?7^,£:}. This was demonstrated 
by Maldacena jl42j but the basic reason is very simple. The non-gaussianity of 
ip comes either from third and higher derivatives of V (through the field equation 
in unperturbed spacetime) or else through the back-reaction (the perturbation of 
spacetime); but the first effect is small |112[ I13(J] by virtue of the flatness require- 
ments on the potential, and the second effect is small because ipo/H^^^ is small |112j . 
This explains why ( with our adopted definition is practically Gaussian by virtue of 
the slow-roll approximation. 
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5.6.2 Multi-component inflation 



The flatness and slow-roll conditions Eqs. ()5.23|) . ()5.24|1 . and ()5.25|1 ensure that 
the curvature of the inflationary trajectories is small during the few Hubble times 
around horizon exit, during which the quantum fluctuation is promoted to a classical 
perturbation. As a result, the initial curvature perturbation in first-order perturba- 
tion theory is still given by the amplitude in Eq. ()5.4U|) and the spectral index in 
Eq. ()5.41|) in terms of the field (p that we defined earlier. 

What about the initial non-gaussianity generated at second order? In the approx- 
imation that the curvature of the trajectories around horizon exit is completely neg- 
ligible, we can safely say that the initial non-gaussianity corresponds to |/nl| <^ 1- 
Confirming this expectation, Seery and Lidsey '2D0^ have calculated the three-point 
correlator of the perturbations in the fields involved in multi-component slow-roll 
inflation. Their result is given by 



(5.50) 



with 



ani 



d^^ 



Bijk{ki, k2, ^3 



ijk 



3 3 -|- cyclic permutations 



(5.51) 



Sjkfsiiki, k2, k^) + cyclic permutations in k and 



(5.52) 



As we will discuss in Chapter IHl where the 5N formalism is used to calculate the 
stochastic properties of C |lH2j , the contribution A/nl of the wavevector dependent 
parameter fij^ to the total /nl is in any case very small, being A/nl generically 
below {lb/2A)fsL^/W^ < 10"^ [141J where fsL is in the range 1/3 < fsL < 11/18. 



5.7 The evolution after horizon exit 
5.7.1 Single-component inflation and (2 

During single- component inflation the curvature perturbation (, as we have defined 
it, does not evolve. From its definition Eq. ()5.7p . the same is true of ^^w. 
In contrast C^, given by Eq. ()5.29|) . will have the slow variation |217j 

c^^ -(2^ - v^my ■ (5.53) 

^^The cyclic permutations in k and (j> in Eq. (|5.52() must be simultaneous, i.e. when exchanging 
indices i and j, for example, fci and ^2 must also be exchanged. Notice also that the calculation 
of Seery and Lidsey's is only valid when the magnitudes of the three wavevectors are roughly 
comparable, so that they exit the horizon at similar epochs. 
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This variation has no physical significance, being an artifact of the definition. 

Using a particular gauge, Acquaviva et. al. [3] have calculated in terms of 
first-order quantities ipi, 5ipi, and their derivatives, and they have displayed the 
result as an indefinite integral 

C^{t) = f A{t)dt + B{t). (5.54) 

Inserting an initial condition, valid a few Hubble times after horizon exit, this be- 
comes 

C2''(^) = C2itin{) + /* A{t)dt + B\l^^ . (5.55) 

In view of our discussion, it is clear that these equations will, if correctly evaluated, 
just reproduce the time dependence of Eq. ()5.53j) . 

The authors of Ref. jH] also present the respective equation for i^^, involving 
only first-order quantities, which is valid also before horizon entry. Contrary to the 
claim of the authors, this classical equation cannot by itself be used to calculate the 
initial value (more precisely, the stochastic properties of the initial value) of C^. In 
particular, it cannot by itself reproduce Maldacena's calculation of the bispectrum. 

It is true of course that in the Heisenberg picture the quantum operators satisfy 
the classical field equations. In first-order perturbation theory, where the equations 
are linear, this allows one to calculate the curvature perturbation without going 
to the trouble of calculating the second-order action |112j (at the nth order of 
perturbation theory the action has to be evaluated to order n-|- 1 if it is to be used). 
At second order in perturbation theory it remains to be seen whether the Heisenberg 
picture can provide a useful alternative to Maldacena's calculation, who adopted the 
interaction picture and calculated the action to third order. 

5.7.2 Multi-component inflation 

During mult i- component inflation the curvature perturbation by definition varies 
significantly along a generic trajectory, which means that non-gaussianity is gen- 
erated at some level. So far only a limited range of models has been investigated 
(iniEHlEniinSISZllinZllini. To keep the spectral tilt within observational bounds, 
the unperturbed trajectory in these models has to be specially chosen, but the choice 
might be justified by a suitable initial condition. 

We shall consider here a calculation by Enqvist and Vaihkonen in Refs. [6T| I13H 
I214j . Following the same line as Acquaviva et. al. [3], they study a two-component 
inflation model, in which the only important parts of the potential are 

V{^, d)=Vh + \mld^ + ]^m'y . (5.56) 

The masses are both supposed to be less than {3/2)H^, so that this is a two- 
component inflation model, and the above form of the potential is supposed to 
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hold for some number AA^ of e-folds after cosmological scales leave the horizon. 
They take the unperturbed inflation trajectory to have -^q = 0, and the idea is to 
calculate the amount of non-gaussianity generated after AA^ e-folds. Irrespective 
of any later evolution, this calculated non-gaussianity will represent the minimal 
observed one (unless non-gaussianity generated later happens to cancel it). 

It is supposed that the condition {^q = 0, as well as the ending of inflation, will 
come from a tree-level hybrid potential, 

V{^, ^)=V,- \ml^' + + ]^m'y + i/^V' • (5.57) 

Like the original authors though, we shall not investigate the extent to which 
Eq. ()5.57|1 can reproduce Eq. ()5.5fij) for at least some number of e-folds. We just 
focus on Eq. ()5.56p . with the assumption = for the unperturbed trajectory. 

Because ■t?o = 0, the unperturbed trajectory is straight, and at first order the 
curvature perturbation ^ is conserved. This is not the case though at second order. 
Adopting the definition C^, the authors of Ref. [HZj give an expression for similar 
to that in Eq. ()5.55|1 describing the evolution of the second-order curvature per- 
turbation on superhorizon scales^^. This equation, in the generalized longitudinal 
gauge, reads (from Eq. (67) in Ref. 



-2{5^,f + mliS^^f + {e- r^^)m,^,V-%{dkdHdid'5'd 



'1, 



~y-^d,d\dk5did^5^iy]dt + [ - V'^di{5did'5d 

(5.58) 

where is the inverse of the Laplacian operator, r]^ = m^/3H'^, and e is defined 

by 

which reduces to the e parameter in Eq. ()5.17|) for t = t^,, being the time when 
cosmological scales exit the horizon. 

Assuming that this expression is correct, we consider the non-gaussianity it may 
generate. Reviewing what it was done in Ref. |214j . we note first that at t = t^, 

SMQ-^MQ, (5-60) 



"'^'^The fields (p and d in Eq. (|5.56l) are supposed to be canonically normalised, which means that 
(p is not the field appearing in the Rigopoulos definition Eq. I|5.26(l of C^. Instead the authors of 
Ref. |()7| give an equivalent definition in terms of the canonically normalised fields. 
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which is a good approximation since at that time the amphtude of the spectrum 
of perturbations of any hght field (j) is ^ H^/2t[. Moreover, assuming slow-roll 
conditions we obtain 



e 



ee 



(5.61) 
(5.62) 
(5.63) 



where we have used = J^^ Hinidt, r]^ = m'^/3H^, and Eq. ()5.59|) . A similar 
expression for the evolution of S(pi is obtained by invoking the constancy of the 
first-order curvature perturbation (i. 

difiit) = SifiiQe-'^^^ . (5.64) 

Assuming that -ffmf, V-fiy ^'^^ almost constants in time, we end up with 



C2{t) - (2 it. 



sH^mp 



1 



dt + 



2r7^)7^][j, (5.65) 



where 



7^ = 3V-^di{dkd''6^id'6^h) + ^V-\di6^id'6^i 



(5.66) 



and we have used the equation of motion S'di + 3Hi^{6{}i + m'^d'di = to go from 
Eq. (IH3H|1 to Eq. 

The order of magnitude for C2^(t) — C2^(t*) is now easily estimated by means of 
the expressions in Eqs. ()5.6H) to ()5.64|) . and by neglecting the scale dependence of 
the non-local terms: 



^rrip ^ Jt, 



t r 



dt 



+ 



mH,^f\SA\' + mnf\S^i\' + V^H,^{\S^^\^] I}, (5.67) 



so that, using Eq. ()5.6()|) to write S-di in terms of (i, 



(5.68) 



It is unlikely that the exponential factor on the right hand side provides any signifi- 
cant enhancement to (2 if f produces most of the curvature perturbation. Therefore, 
the overall slow-roll factors give the actual magnitude. We have to remember that. 
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in this case, the right hand side is uncorrelated with the inflaton perturbation 6(j) 
which generates (f^. Thus, Eq. ()5.38p as opposed to Eq. ()5.33|) apphes, and the asso- 
ciated /nl would be ~ 10^^ 0{e^,ri^,ril), which is extremely small. If the observed 
( has a non-gaussian part (2 equal to Eq. ()5.f)8|) and a gaussian part generated 
mostly after inflation, one can obtain |/nl| > 1 by choosing i]^ > 0.26, 77,9 = r]^/2, 
N = 70, and (2 = lO'^C- 

As we will see in the next chapter, where the Enqvist and Vaihkonen model is 
studied by means of the non perturbative 6N formalism (see Subsection I6.5.3|) . the 
expression in Eq. ()5.(i8|l disagrees with the one calculated using the SN formalism 
through the appearance of non-local terms I 132|I133] . though the order of magnitude 
is similar |214j . We point out that the possible source of discrepancy is the use of a set 
of cosmological perturbation theory equations in Ref. [HZl based on those presented 
in Ref. [2]. A calculation made with another set of cosmological perturbation 
theory equations that do not involve non-local terms [144' reproduces exactly the 
result found using the SN formalism (see Subsection lfi.5.3|l . 

5.7.3 Preheating 

Now we turn to the possibility that significant non-gaussianity could be generated 
during preheating. Preheating is a stage of non-perturbative explosive resonant de- 
cay of scalar fields which might occur between the end of inflation and reheating 
|^l^l^l212j . the latter being taken to correspond to the decay of individual parti- 
cles which leads to more or less complete thermalisation of the Universe. Preheating 
typically produces marginally-relativistic particles, which decay before reheating. 

It was suggested a long time ago |2Il 122] that preheating might cause the cosmo- 
logical curvature perturbation to vary at the level of first-order perturbation theory, 
perhaps providing its main origin. More recently it has been suggested [001 EH 1121 
that preheating might cause the curvature perturbation to vary at second order, 
providing the main source of its non-gaussianity. 

If the separate universe hypothesis is correct, a variation of the curvature per- 
turbation during preheating can occur only in models of preheating which contain 
a non-infiaton field that is light during inflation. This is not the case for the usual 
preheating models that were considered in [00110111021, and accordingly one does not 
expect that significant non-gaussianity will be generated in those models^^. This is 
not in conflict with the findings of [001 IMl I02| because the curvature perturbation 
is not actually considered there. Instead the perturbation ip^"^ in the longitudi- 
nal gauge is considered, which is only indirectly related to ( by Eqs. ()5.7|) . ()5.11|) 
and ()5.12p We conjecture that non-gaussianity for the curvature perturbation 
on cosmological scales is not generated in the usual preheating models, but that 

^^The preheating model considered in |62| contains a field which may be heavy or light; we refer 
here to the part of the calculation that considers the former case. 

^^The slices of the longitudinal gauge are orthogonal to the threads of zero shear, and ?/;^^ on 
them is very different from the curvature perturbation ^. 
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instead the curvature perturbation remains constant on cosmological scales. This 
should of course be checked, in the same spirit that the constancy of the curvature 
perturbation was checked at the first-order level |113j . 

The situation is different for preheating models which contain a non-inflaton field 
that is light during inflation. At least three types of models have been proposed with 
that feature |S1 1^ 1^ 1^ llUUj . Except for [52] only the magnitude of the curvature 
perturbation has been considered, but in all three cases it might be that significant 
non-gaussianity is also generated. 

5.7.4 The curvaton scenario 

In the simplest version of the curvaton scenario jlH8| (see Subsection 

I2.4.2j) . the curvaton field a is ultra-light during inflation, weakly coupled, and has 
no significant evolution until it starts to oscillate during some radiation-dominated 
era. Until this oscillation gets under way, the curvature perturbation is supposed 
to be negligible (compared with its final observed value). The potential during the 
oscillation is taken to be quadratic, which will be a good approximation after a few 
Hubble times even if it fails initially. The curvature perturbation is generated during 
the oscillation, and is supposed to be conserved after the curvaton decays. Here 
we give a generally-valid formula for the non-gaussianity in the curvaton scenario, 
extending somewhat the earlier calculations. 

The local energy density of the curvaton field is given by [c.f. Eq. (j2.5Hj) ] 



where cra(?7, x) represents the amplitude of the oscillations and nia is the effective 
mass. It is proportional to 0(77, x)~'^ where a is the locally-defined scale factor. 
This means that the perturbation dpa/pao is conserved if the slicing is chosen so 
that the expansion going from one slicing to the next is uniform |14Uj . The fiat 
slicing corresponding to ip^^ = has this property |127| ll4Uj and accordingly 5prj 
is defined on that slicing (see Subsection I2.4.2|l . 

Assuming that the fractional perturbation is small (which we shall see is de- 
manded by observation) it is given by 



where we have extended to second order the Eq. ()2.65j) . We first assume that cr(x) 
has no evolution between inflation and the onset of oscillation. Then daa/oao "will 
be equal to its value just after horizon exit, which we saw earlier will be practically 
gaussian. 

The total density perturbation is given by 



(5.69) 




(5.70) 




(5.71) 
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where Q = Pao/Ptotaio oc a is the fraction of energy density contributed by the 
curvaton. Adopting the sudden-decay approximation, the constant curvature per- 
turbation obtaining after the curvaton decays is given by Eqs. fl5.11|) and ()5.13|1 . 
evaluated just before curvaton decay and with ip = 0. In performing that calcula- 
tion, the exact expression Eq. ()5.7()|1 can, without loss of generality, be identified 
with the first-order part Sp^j^/p^g, the second- and higher-order parts being set at 
zero. 

Adopting the first-order curvature perturbation in Eq. ()5.1H1 . one finds |138j 
non-gaussianity coming from the second term of Eq. ()5.7()j) , 

/nl = -|:, (5.72) 

with 

r = , (5.73) 

evaluated just before decay. Going to the second-order expression one finds [T2j 
additional non-gaussianity. The final non-linearity parameter /nl = /nl^ + 5/3 
is given by 

If Q just before the curvaton decay (f2dec) is much less than 1 {Qdec 1) then 
/nl is strongly negative and the present bound on it requires Qdcc ^ 0.01 (combined 
with the typical value ^ 5 x 10~^, this requires 6pa/pao <^ 1 as advertised). If 
instead fl^ec = 1 to good accuracy, then /nl = +5/4. Either of these possibilities 
may be regarded as generic whereas the intermediate possibility (|/nl| ~ 1 but 
/nl 7^ 5/4) requires a special value of f2dec just a bit less than 1. 

Finally, we consider the case that a evolves between horizon exit and the era 
when the sinusoidal oscillation begins. If (Xa (the amplitude of oscillation at the 
latter era) is some function g{a^:) of the value a few Hubble times after horizon exit, 
then 

1 

6(ra = g'6(T, + -g"{6(r,Y , (5.75) 

where the prime means derivative with respect to a*. Repeating the above calcula- 
tion one finds 

5 5 5 / gg"' 

— — r IH ; 

3 6 4r V g'- 

The final term is the first-order result (given originally in |126j ). the middle term is 
the second-order correction found in ^^I, and the first term converts from to 
/nl- 

5.7.5 The inhomogeneous reheating scenario 

The final scenario that has been suggested for the origin of the curvature per- 
turbation is its generation during some spatially inhomogeneous reheating process 



/nl = 77 + - - I 1 + ^ • (5-76) 
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|S3 Ea EH (see also Refs. [HSl EHl CSOl CHI IHH 1221 )• Before a reheating process 
the cosmic fluid is dominated by matter (non-relativistic particles, or small scalar 
field oscillations which are equivalent to particles) which then decay into thermalised 
radiation. At least one reheating process, presumably, has to occur to give the initial 
condition for Big Bang Nucleosynthesis, but there might be more than one. 

The inhomogeneous reheating scenario in its simplest form supposes that the 
curvature perturbation is negligible before the relevant reheating process, and con- 
stant afterwards. The inhomogeneous reheating corresponds to spatial fluctuations 
in the decay rate of the inflaton field to ordinary matter, which lead to fiuctuations 
in the reheating temperature. The coupling of the infiaton to normal matter is de- 
termined by the vacuum expectation values of scalar fields in the theory. If those 
fields are fight they will fiuctuate leading to density perturbations through the de- 
scribed mechanism. Inhomogeneities in the infiaton decay rate lead to a a spatially 
varying value (a perturbation) of the local Hubble parameter ifreh(x) at the decay 
epoch (or equivalently of the local energy density). 

In contrast with the curvaton scenario, where the form p^- can reasonably be taken 
as Po- oc 0"^, the inhomogeneous reheating scenario does not suggest any particular 
form for /7i.eh(x)- Depending on the form, the inhomogeneous reheating scenario 
presumably can produce a wide range of values for /nl I224j . 

5.8 Conclusions 

We have examined a number of scenarios for the production of a non-gaussian pri- 
mordial curvature perturbation, presenting the results with a unified notation jlHlj . 
These are the single-component infiation, mult i- component infiation, preheating, 
curvaton, and inhomogeneous reheating scenarios. Although the trispectrum may 
give a competitive observational signal IHHl I215j . we have focused only on 
the bispectrum which is characterised by the parameter /nl. In all cases our treat- 
ment is based on existing ones, though we do not always agree with the original 
authors. 

Table 5.1: Non-gaussianity according to different scenarios for the creation of the 
curvature perturbation. For the simplest curvaton scenario, /nl = +5/4 is a 
favoured value. 

Scenario |/nl| < 1 |/nl| - 1 /nl < -1 /nl > 1 

Single-component infiation yes no no no 

Multi-component infiation likely possible possible possible 
Simplest curvaton scenario unlikely likely likely no 



The preheating and inhomogeneous reheating scenarios cover a range of pos- 
sibilities, which have not been fully explored but which can presumably allow a 
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wide range for /nl- The same is true of multi-component inflation, except that 
extremely large values comparable with the current bound |/nl| ^ 10^ seem rela- 
tively unlikely. In contrast, the simplest curvaton scenario can produce a strongly 
negative value (even violating the current bound). However, in the important spe- 
cial case where the curvaton dominates the energy density before it decays, it gives 
precisely /nl = +5/4. Finally, for the single-component inflation case, Maldacena's 
calculation combined with current constraints on the spectral tilt show that it has 
magnitude less than 10^^. These result are summarised in the Table l5?Tl 

In the near future, results from WMAP |22Hj or elsewhere may detect a value 
|/nl| 3> 1. If that does not happen, then PLANCK |171j or a successor will either 
detect a value |/nl| ~ 1, or place a bound |/nl| < 1- The precise level at which this 
will be possible has yet to be determined because it would require a second-order 
calculation of all relevant observational signatures. The example of the simplest cur- 
vaton scenario, where /nl = +5/4 is a favoured value, shows that such a calculation 
and the eventual observations will be well worthwhile. 
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Chapter 6 



The inflationary prediction for 
primordial non-gaussianity 

6.1 Introduction 

In this chapter we present for the first time a powerful method to calculate the nor- 
malisation /nl of the bispectrum in slow-roll inflation, by means of the knowledge 
of the evolution of a family of unperturbed universes |132j . The wavevector depen- 
dence of /nl will be, in general, negligible |14H I2()()j compared with the (possibly 
big) contribution coming from the evolution of the unperturbed universes. This 
method will be applied to selected examples. In particular we will see how the level 
of non-gaussianity in the curvaton scenario (see Subsection I5.7.4p and the second 
order curvature perturbation in the hybrid model of Enqvist and Vaihkonen |144j 
are successfully reproduced. 

The primordial curvature perturbation of the Universe, is already present a few 
Hubble times before cosmological scales start to enter the horizon |1121 ll3Uj . Its 
time-independent value at that stage seems to set the initial condition for the sub- 
sequent evolution of all cosmological perturbations. As a result, observation probes 
the stochastic properties of (, which is found to be almost gaussian with an almost 
scale-invariant spectrum. 

According to present ideas ( is supposed to originate from the vacuum fluctua- 
tions during inflation of one or more light scalar fields, which on each scale are pro- 
moted to classical perturbations around the time of horizon exit [11 1751 IHHll22|n^lj . 
One takes inflation to be almost exponential (quasi de Sitter spacetime) correspond- 
ing to a practically constant Hubble parameter i^inf, and the effective masses of the 
fields to be much less than H^,. This ensures that the fields are almost massless and 
live in almost unperturbed quasi de Sitter spacetime, making their perturbations 
indeed almost gaussian and scale invariant (see Subsections 12.4.1) and This 
automatically makes ( almost scale invariant, (see Subsections 12.4.21 and and 
can (though not automatically |13H I1H8] ) make it also almost gaussian. 

All of this is of intense interest at the present time, because observation over 
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the next few years will rule out most existing scenarios for the generation of C, by 
detecting or bounding the scale dependence and non-gaussianity of (. We will now 
describe a general procedure for calculating the level of non-gaussianity, by means 
of the SN formalism [T2Z1 CnE] (see also Refs. dm EEl EH]). 



6.2 Defining the curvature perturbation 



Perturbations of the observable Universe are defined with respect to an unperturbed 
reference universe, which is homogeneous and isotropic (a FRW universe) (see Sec- 
tion Its line element may be written as 



ds'^ = —dt^ + o?{t)5ijdx^dx^ 



(6.1) 



defining the unperturbed scale factor a(t), time t, and the Cartesian spatial coordi- 
nates X. 

The curvature perturbation is only of interest after the universe has been smoothed 
on some scale {k/a)~^ much bigger than the horizon H~^. To define it, one takes the 
fixed-t slices of spacetime to have uniform energy density, and the fixed-x worldlines 
to be comoving. The spatial metric is [SHI iSl El 1123 CSIl II121 CHIl CHSl 
IMl] [c.f. Eq. Q] 



gij = a'(t)e2^(*'")7,,-(t,x) = a'(t, x)7,,(t, x) 



(6.2) 



In this expression, '~fij{t,x.) has unit determinant, so that a volume of the Universe 
bounded by fixed comoving spatial coordinates is proportional to the locally defined 
scale factor d^{t,x). In the inflationary scenario the factor just accounts for 
the tensor perturbation, but its form is irrelevant here (see Subsections 15.2.21 and 
I5.6.H) . According to this definition, ( is the perturbation in In a. Only the spatial 
variation of ( is significant, and to make contact with observation we can work with 
its Fourier components in a box a bit bigger than the observable Universe, setting 
the zero mode equal to zero so that ( has vanishing spatial average. 

One can also consider a slicing whose metric has the form in Eq. ()6.2|) without 
the ( factor, which we call the fiat slicing. Starting from any initial fiat slice at time 
tini, let us define the amount of expansion 



iV(t,x) 



In 



d{t) 



to a final slice of uniform energy density. Then [llUt I127t I196| I2im I211j 



C(t,x) = 5iV = iV(t,x)-iVo(t) 



where 



Noit) = In 



■ a{t) 



(6.3) 



(6.4) 



(6.5) 
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is the unperturbed amount of expansion. 

To make use of the above formahsm we assume that in the superhorizon regime 
{aH ^ /c), the evolution of the Universe at each position (the local evolution) is 
well approximated by the evolution of some unperturbed universe [121 1127| I14()| I197[ 
I222j . This 'separate universe' assumption will presumably be correct on cosmological 
scales because these scales are so big |14Uj . 

By virtue of the separate universe assumption, N{t, x) is the amount of expansion 
in some unperturbed universe, allowing C, to be evaluated knowing the evolution of 
a family of such universes. For a given content of the Universe it can be checked 
using the gradient expansion I127| I179| method, but we do not wish 

to assume a specific content. 

The separate universe assumption leads also to local energy conservation. Indeed, 
using the uniform density slicing, and remembering that a determines the expansion, 

pit) = -3H[p{t) + Pit, x)] = -?>(h + c) [pit) + P(t, x)] , (6.6) 

where H = a/a, p is the energy density, and P is the pressure. During any era when 
P is a unique function of p (the adiabatic condition), P is uniform on the chosen 
slicing; then ( vanishes (because it is uniform and its spatial average vanishes) 
so that ( is conserved (see Subsection I2.2.3|l . This consequence of the separate 
universe assumption was first recognised in full generality in Refs. |127[I179] (see also 
Ref. |193j for the case of inflation, Refs. |14UH222] for the case of linear perturbation 
theory, and Refs. |lU6t llU7j for a coordinate-free treatment). 



6.3 The inflationary prediction 

The evolution of the observable Universe, smoothed on the shortest cosmological 
scale, is supposed to be determined by the values of one or more light scalar fields 
when that scale first emerges from the quantum regime a few Hubble times after 
horizon exit. Defined on a flat slicing, each field (pi at this epoch will be of the form 

0i(x) = + (50i(x). 

Because quasi exponential inflation is assumed, and only light fields are consid- 
ered, it is a good approximation to take the S(j)i to be almost massless fields living 
in unperturbed quasi de Sitter spacetime |112| 1130] . In these circumstances the 
perturbations 6(f)i generated from the vacuum are almost gaussian, with an almost 
flat spectrum whose amplitude is |j^ EZj 

M^^^- (6-7) 

Now we invoke the separate universe assumption, and choose the homogeneous 
quantities 0jp to correspond to the unperturbed universe. Then Eq. ()6.4p for ( 
becomes 

at, x) = N{p{t), 01 (x), 02(x), ■ ■ ■) - N{p{t), 02o, ■ ■ ■) • (6.8) 
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In this expression, the expansion is evaluated in an unperturbed universe, from 
an epoch when the fields have assigned values to one when the energy density has an 
assigned value p. This expression |lin| I127| I196[ 12 10^ I211j allows one to propagate 
forward the stochastic properties of C to the epoch when it becomes observable, 
given those of the initial field perturbations. 

Since the observed curvature perturbation is almost gaussian, it must be given 
to good accuracy by one or more of the linear terms 

C(t,x)^^iV,(t)50,(x), (6.9) 

i 

where we use the notation 

dN 

N, = TT-, (6.10) 

Nij = TT-jr-: (6.11) 

with the field perturbations being almost gaussian. Here we include for the first 
time the quadratic terms |132j 

C(t,x) = Y,N,{t)S^i + ^Y.N,^{t)S<f>Mj ■ (6.12) 

i ij 

They may be entirely responsible for any observed non-gaussianity if the field per- 
turbations are gaussian to sufficient accuracy^. 



6.4 Non-gaussianity 
6.4.1 The bispectrum 

The stochastic properties of the perturbations are specified through expectation 
values which, according to the infiationary paradigm, are taken with respect to the 
time independent (Heisenberg picture) quantum state of the Universe (to be precise, 
the quantum state of the universe before it somehow collapses to give the observed 
Universe). Focusing on (, we consider Fourier components, 

^^--y^C(i,x)exp(?k-x) . (6.13) 

The stochastic properties of a gaussian perturbation are specified entirely by the 
spectrum Vc^{k) = A^i^k / aHi^f)"''^ , defined through 

^k^c.) = Pdk)SHki + k2) , (6.14) 

^Here and elsewhere, we are not displaying a homogeneous term needed to make the spatial 
average of C vanish. 



102 



and 

Vdk)^^,Pdk). (6.15) 

From Eqs. dHUj) and 

^?=(§)'i:iV/. (6-16) 

Non-gaussianity is defined through higher correlations. We consider only the 
three-point correlation^. It defines the bispectrum 5^ through (see Subsection l5.5.1|) 

(CkiCk.Ck3) = i2n)-'/'Bdh, h, h)S%k, + k2 + kg) . (6.17) 
Its normalisation is specified by a parameter /nl according to jEl I103| I142j 
6 

= —-fTsshiki, k2, k3)[P(^{ki)P(^{k2) + cyclic perturbations] . (6.18) 



In first-order cosmological perturbation theory the gauge-invariant gravitational po- 
tential $g during matter domination before horizon entry is $g = — (5/3)C, and our 
definition of /nl coincides with the definition |l()3j 

Bq>g = 2/NL(fci, k2, A;3)[P$^(/ci)P$g(A;2) + cychc perturbations] . (6.19) 

At second-order these definitions of /nl differ (see for example Refs. [HI UHl HZl UHl 
[inil2ni- See also footnotes El and H in Chapter E)). 

We shall take (A;) and /nl to be evaluated when cosmological scales approach 
the horizon and ( becomes observable. Observation gives \A^\ ^ 5 x 10~^ j2U7j . and 
|/nl| ^ 100 |102j . Absent of a detection, this will eventually come down to roughly 

i/nlI < 1 m. 

Ignoring any non-gaussianity of the 6 (pi, our formula in Eq. ()6.12j) makes /nl 
practically independent of the wavenumbers. Indeed, generalising the result found 
in Ref. we have calculated (see Appendix 

-fNL= \_ '^^^f +ln(fcL)^ ^'^^_''^ '^,3 ■ (6-20) 



5 



2 



In deriving this expression we used the amplitude of the spectrum A^^^ ^ H^./2t[ 
of the field perturbations, and used Eq. ()6.16p to eliminate H^, in favour of A^. As 
discussed in Ref. [32 , the logarithm can be taken to be of order 1, because it involves 
the size of a typical scale under consideration, relative to the size L of the region 
within which the stochastic properties are specified. Except for the logarithm, /nl 
is scale independent if the field perturbations are gaussian. 

^The four-point correlation may give a competitive observational signature and can be calculated 
in a similar fashion 132 ES HH ESI • 
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If only one 5 (pi is relevant, Eq. ()6.12|) becomes 



C(t,x) = iV,50, + liV,,(5</),)2, (6.21) 



and because the first term dominates, Eq. ()6.20|) becomes 



(6.22) 



In this case, /nl may equivalently be defined |lU3j by writing 

C = Cg-^/NLC|, (6.23) 

where Cg is gaussian. 

To include the possible non-gaussianity of the 5 (pi, we define the bispectra Bijk 
of the dimensionless field perturbations {27t /H^)S(j)i and their normalisation fijk, in 
exactly the same way that we defined and /nl [cf. Eqs. ()5.50|) and ()5.51|) ]. 
These bispectra add the following contribution to /nl in Eq. ()6.20|) (see Appendix 
EI): 

A/nl = . (6.24) 



The fijk, generated directly from the vacuum fluctuation, will depend strongly on 
the wavenumbers [cf. Eq. ()5.52j) ]. 



6.4.2 Cosmological perturbation theory 

In the superhorizon regime the non-linear theory jl27j that we have used is a com- 
plete description. The basic expression in Eq. ()6.8j) is non-perturbative, giving ((t, x) 
in terms of the initial fields and the expansion of a family of unperturbed universes. 
The second-order expansion in Eq. ()6.12j) is a matter of convenience. As we shall 
see it seems to be adequate in practice, but Eq. ()6.8|) would still be apphcable if the 
expansion converged slowly or not at all. 

Cosmological perturbation theory (CPT) is completely different. It is applicable 
both inside and outside the horizon, being at each instant a power series in the 
perturbations of the metric and the stress-energy tensor, together with whatever 
variables are needed to completely specify the latter and close the system of equa- 
tions. During inflation these variables are the components of the infiaton, while 
afterwards they may involve oscillating fields and a description of the particle con- 
tent. First-order CPT is usually adequate and can describe non-gaussianity at the 
level |/nl| 3> 1 which has to be generated by the second-order term in Eq. ()6.12j) . 
Second-order CPT^ is generally needed only to handle non-gaussianity at the level 
|/nl| ~ 1- 

^For a full set of CPT equations see e.g. Ref. |144| . 



104 



Quantised CPT is needed to calculate the stochastic properties of the initial field 
perturbations which are the input for our calculation. The slow-roll spectrum 
in Eq. ()6.7|1 comes from the first-order calculation. The bispectrum is a second-order 
effect and has, in the context of slow-roll infiation, been calculated in Refs. |142[ I2()()j 
[c.f. Eqs. ()5.5()|) - ()5.52|l ]. It is shown elsewhere jl41j that A/nl is in this case 
negligible compared to 1, being generically below {15/2A)fsL\/fTC^ ^ 10~^ with 
1/3 < fsL < 11/18 [c.f. Eqs. ()5.52|) and ()6.24|) ]"^. Higher order correlators have not 
been calculated yet and would give an additional contribution to Eq. ()6.2n|l . which 
presumably is also negligible. A few single field non slow-roll models P HHl HHHl EHHI 
have been investigated where it is found that A/nl could be much bigger than 1. 
From now on we take the 50j to be gaussian. 

In the regime aH ^ k, perturbation theory must be compatible with Eq. ()6.12j) . 
In particular, the non-local terms, present at second order for a generic perturbation, 
must be absent for (. Finally, CPT is needed to evolve the perturbations after 
horizon entry, but that is not our concern here. In the following, we apply our 
formalism to calculate /nl in various cases and compare it with the CPT result 
where that is known. 

6.5 The SN formalism in some multi-component 
models 

6.5.1 A two- component inflation model 

As a first use of Eq. ()6.20|) we consider the two-component infiation model of Kadota 
and Stewart |9(H l9T] , estimating for the first time the non-gaussianity which it pre- 
dicts. The model works with a complex field $, which is supposed to be a modulus 
with a point of enhanced symmetry at the origin. The scalar potential is given by 

V = Vh- m2|$|2 + -Am^[<l>'^ + <l>*^] + + 1)AW|$|^ , (6.25) 
3 2 

with A being fixed so that the vacuum energy vanishes at the minimum of the 
potential, and u = {1,2,3, ...}. Writing 

$ = |<l>|e*^ (6.26) 

the tree level potential has a maximum at $o = and depends on both |$| and 6. 
A one-loop correction turns the maximum into a crater and infiation occurs while ^ 
is rolling away from the rim of the crater (see Fig. 16.11) . The curvature perturbation 

^In Ref. 183J (see also Refs. [831 1182) ') Rigopoulos et. al. calculated the three-point correlator 
of C in general multi-component inflationary models using a stochastic approach. Their result is 
quite puzzling since the wavevector dependence seems to be significant even well after horizon exit. 
Reconciliation between the approach followed in Ref. |188j and ours |182j is desirable jl84j . 
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Figure 6.1: Equipotential layers and trajectory of the field $ = + i^i)/\/2 in 
the Kadota and Stewart's model of Refs. jHOlini]. Figure taken from Ref. PT| . 



is supposed to be constant after the end of slow-roll inflation. For 6*0 ^ 6c, with 6c 
being a parameter of the model, it is found that [HI] 



which is too small ever to be observed. 
6.5.2 The curvaton model 

In the curvaton model |1H8| I1H9| I159j (see also Refs. [011 1116^ I156j and Subsection 
I2.4.2|) the curvature perturbation ( grows, from a negligible value in an initially 
radiation dominated epoch, due to the oscillations of a weakly coupled light field a 
(the curvaton) around the minimum of its quadratic potential 




(6.27) 



Through the first term of Eq. ()6.2Up 



6_o 

6c 



(6.28) 



y,(t,x) = -m>2(t,x). 



(6.29) 
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N{pdec, Pose, (y*) = ^ In ( = ^ In 




(6.31) 



where mo- is the curvaton effective mass. Due to the oscillations, the initially negli- 
gible curvaton energy density redshifts as 

Po(t, x) ^ \mlal{t, x) oc a'''{t, x) , (6.30) 

where Oa represents the amplitude of the oscillations. Meanwhile the radiation 
energy density pr redshifts as a""^. Soon after the curvaton decay, the standard Hot 
Big-Bang is recovered and C, is assumed to be conserved until horizon reentry. 

To calculate /nl using Eq. ()6.20|) we first realise that (the unperturbed value 
of cr a few Hubble times after horizon exit) is the only relevant quantity since the 
curvature perturbation produced by the infiaton, and imprinted in the radiation 
fluid during the reheating process, is supposed to be negligible. Thus, Eq. ()6.22p 
applies. Second, we can redefine N as the number of e-folds from the beginning 
of the sinusoidal oscillations to the curvaton decay. This is because the number of 
e-folds from the end of inflation to the beginning of the oscillations is completely 
unperturbed as the radiation energy density dominates during that time. Thus, 
is now a function of three variables 

where g = cXosc is the amplitude at the beginning of the sinusoidal oscillations as 
a function of a^. Here the curvaton energy density just before the curvaton decay 
Po-dec is expressed in terms of the total energy density pdec at that time, the total 
energy density at the beginning of the sinusoidal oscillations posc, and g by 

P^d.. = lral[gia.)rl P'- . (6.32) 

^ \ Pose / 

After evaluating d/da^ = g'd/dg, at fixed pdec and posc, we obtain 

2 a' 

Na, = -r- , (6.33) 
"J 9 

where 

r = ^-^^ , (6.34) 

being pr^^^ the radiation energy density just before the curvaton decay, giving 

H^r q' 

in agreement with first-order cosmological perturbation theory in the sudden decay 
approximation [1381 1139t 1159] (see Subsection I2.4.2|) . Differentiating again we find 
from Eq. (j6:22|l 

which nicely agrees with the already calculated /nl using first- and second-order 
perturbation theory (see Refs. pT j I126| |13H I138j and Subsection 15. 7. 4|1 . 
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6.5.3 Another two-component model 

Finally we consider the two-component inflation model of Ref. (see also Refs. |13H 
I214j and Subsection 15 . 7. 2j) . For at least some number N of e-folds after cosmological 
scales leave the horizon, the potential in Eq. ()5.5fij) is written as 

with the flrst term dominating, and rj^ and r^.^ being the slow-roll rj parameters. The 
idea is to use Eq. ()6.2()|) to calculate the non-gaussianity after the e-folds which, 
barring cancellations, will place a lower limit on the observed non-gaussianity. 

The slow-roll equations give the background fleld values ipo{N) and 'do{N) after 
N e-folds, in terms of those obtaining just after horizon exit 

ifoiN) = ^oeM-Nv^), (6.38) 
MN) = ^oexpi-Nri^). (6.39) 



This gives 



K,(iV, ^0, ^o) = vjl + Iv.^e-'''^^ + Iv^^e-'^'A , (6.40) 



2 rrip 2 rrip 



and allows us to calculate the derivatives of N with respect to (fo and at flxed 
V. Focusing on the case 'do = considered in Ref. [HZj, we flnd 

C = ^ - ^ f ^ V + ^e2^(^--) (^]\ (6.41) 

in agreement with the second-order perturbation calculation of Ref. [144^. If the 
observed ( has a non-gaussian part Co- equal to the last term of Eq. ()6.4H) and a 
gaussian part generated mostly a/ter inflation, one can obtain |/nl| > 1 by choosing 
r]^ > 0.26, r]^ = r/^/2, N = 70, and (a = lO'^C- 

Our calculated expression for the coeflicient of is in disagreement with the 
one found in Ref. |214j [c.f. Eq. ()5.68|) ] which uses a set of CPT equations based on 
those presented in Ref. After converting the variable used there |1311 1217] to 
our ( (see Subsection 15 . 3 . 3|) . these equations give ( in terms of flrst-order quantities, 
but they contain non-local terms involving the inverse Laplacian |133j . Comparison 
with our non-linear expression in Eq. ()6.8j) shows that such terms must cancel if 
correctly evaluated. 

^The initial calculations of the structural form of ( in this model using CPT [^3 \l'M\ were in 
gross conflict with Eq. H6.41|l . This is because the time evolutions of ip and "i? outside the horizon 
were not considered. The sources of discrepancy were recognised in Ref. |214| . showing that the 
actual order of magnitude is in agreement with Eq. (|5.41|) except for the presence of non-local 
terms (see Subsection l5.7.2|l . 
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6.6 Conclusions 



The 6N formalism is a non-perturbative approach to calculate the curvature per- 
turbation ( at all orders, in terms only of background quantities that describe the 
evolution of a family of unperturbed universes. Such a formalism was originally 
introduced to calculate the spectrum of ( at first order 196, 2101 l^llj (see also 
Ref. |llUj ) in multi-component inflationary models. Now, with the increasing inter- 
est in the non-gaussian features of ( in both single- and multi-component inflationary 
models, the formalism has been extended to calculate the curvature perturbation 
at second order (2 and the normalisation /nl of the bispectrum |132j . The 6N for- 
malism relies on the separate universe assumption, which says that on superhorizon 
scales the Universe behaves locally as if it were unperturbed I127| 11401 ^W7\ \'2'2'2\ . 
and on the intrinsic gaussianity of the fields 0j involved. The quantities (2 and /nl 
are easily given in terms of the first and second derivatives of the unperturbed num- 
ber of e-folds N, from an epoch when the fields have assigned values to one when 
the energy density has an assigned value p, with respect to the unperturbed fields 
(pig . The possible intrinsic non-gaussianity of the fields (pi would lead to an additional 
contribution to /nl, highly wavevector dependent pOti] but in any case negligible 
compared to 1 |141j . The 6N formalism reproduces |1H2[I2()()'] the well known results 
in single-component inflation |142j , in the curvaton scenario ^3 1126| ll^H I138j , and 
in the 'hybrid' model of Enqvist and Vaihkonen |144j . In addition, the Kadota and 
Stewart's modular inflation model [HHl 1^ has served as an example of the power 
of this formalism. The 6N formalism is an interesting alternative to cosmological 
perturbation theory where, order by order, the relevant expressions to calculate (2 
and /nl tend to be more and more complicated (see for example Refs. [SI IFTT j I144j ). 
Nevertheless, it is true that the latter is valid on all the scales, while the former 
is only valid in the superhorizon regime. Fortunately, the subhorizon effects are 
negligible |14H I200j making the 6N formalism reliable and completely independent 
of cosmological perturbation theory. 
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Chapter 7 
Conclusions 



Prior to the standard Hot Big-Bang, a period of accelerated expansion seems to 
have been crucial [7^]. Not only does this period solve the classical problems of 
the Big-Bang cosmological model, namely the horizon, flatness, and unwanted relics 
problems, but it also amplifies the fluctuations in the light scalar fields (pi living 
in the Friedmann-Robertson- Walker spacetime [El EH El EHl 11121 CESl ITRKllTm] . 
This inflationary process serves also for the scalar field fluctuations to become clas- 
sical soon after horizon exit [H |7S1 El ll^^l I124j . giving birth to the primordial 
perturbations in the energy density that generate the temperature anisotropies in 
the cosmic microwave background radiation, and, through gravitational collapse, 
the large-scale structure observed today. Such primordial perturbations can be de- 
fined perturbatively on a homogeneous and isotropic background, but the freedom 
to choose the perturbed coordinate system makes them gauge dependent. To char- 
acterize adequately the primordial perturbations, we introduce the gauge-invariant 
curvature perturbation ( ^l] which represents the intrinsic spatial curvature on 
slices of uniform energy density (or slices with zero flow of energy). While the pres- 
sure is a unique function of the energy density, ( is conserved at all orders |127|ll79j . 
which makes this an ideal quantity. In this thesis we have explored some of the the- 
oretical and statistical aspects of the origin of the large-scale structure, such as the 
two most known scenarios to generate ( (the inflaton and the curvaton scenario), 
the required inflationary energy scale in those scenarios, and the non-gaussianity 
associated to (. 

In general ( depends on the perturbations in the scalar fields (pi during the 
inflationary period, whose spectra Vs^pX^) are generically the same for all kinds of 
quasi exponential expansion either the respective field dominates the energy density 
or not. The only appreciable difference is in the way the scale dependence is given in 
terms of the Hubble parameter H^, a few Hubble times after horizon exit and the mass 
m^- of the respective field. However, the spectrum of (, Vi^{k), varies significantly 
among the different possible scenarios for the origin of the large-scale structure, 
although in all cases it is almost gaussian and scale invariant. We have studied in 
Chapter |2]two different scenarios for the origin of (, the inflaton |5l lll2[ITT71ll3nj and 
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the curvaton scenario |138|ll39lll59j (see also Refs. [BHllll6lfl56j ). pointing out their 
different signatures and connecting them with the amphtude of gravitational waves 
produced in each scenario. Both scenarios have their advantages and drawbacks. 
For example, the single- component inflaton scenario presents a consistency relation 
that relates the amplitude of 'Pc_{k) with the amplitude of the gravitational waves 
spectrum Vrik) Pl lUSllTTT|lll2t[T^ . The possible detection of gravitational waves 
is consistent with this scenario and would serve as a smoking gun if the consistency 
relation is satisfied. The drawback is that the generation of C, by the inflaton field 99 
(the field which drives inflation) imposes severe constraints on the theoretical model 
building [101 llfilj . In contrast, the curvaton scenario circumvents the latter 
problem since the generation of C is assigned to a weakly coupled scalar field a (the 
curvaton) different to the inflaton. The drawbacks are that there is no consistency 
relation in this case, and that the scenario is inconsistent with a detectable level of 
gravitational waves |138j . 

There exists several well motivated inflationary models that locate the infla- 
ton field within a particle physics framework ^1251 1130] . Most of these models are 
however unrealistic due to the strong constraints mentioned above. In particular, 
the generation of adiabatic perturbations is almost inconsistent with the low in- 
fiationary energy scale, given by if*, required to identify the inflaton with one of 
the many flelds present in supersymmetry. The curvaton scenario comes to rescue 
these models, allowing for a much lower inflationary energy scale. But how low may 
this energy scale be in the curvaton scenario?. In Chapter El we have discussed the 
lower bound on H,, in the simplest curvaton setup, showing that it is high enough 
|126j {H^ > 10^ GeV) to fail at rescuing the inflationary models whose potentials are 
generated by some mechanism of gravity-mediated supersymmetry breaking where 
if* ~ 10'^ GeV is required. The general conditions to obtain low scale inflation in the 
curvaton scenario have been given in that chapter [ST], in terms of three quantities 
e, /, and 5, that parameterize respectively the evolution of the curvaton fleld from 
the time of horizon exit to the beginning of its oscillations, the effective curvaton 
mass ificr at the end of a phase transition with respect to the Hubble parameter 
at the same time, and the time of the phase transition with respect to the time of 
horizon exit. In Chapters El and |3] we have invoked the 'heavy curvaton' picture 
|126j , deflned as the setup where a suddenly increases its mass at the end of a phase 
transition much later than inflation |149j . The mass increment is given by the cou- 
pling of a (parameterized by the constant A) with another fleld which acquires a 
large vacuum expectation value at the end of the same period. Thus, the smallness 
of the parameters / and 5 are exploited to allow for an inflationary scale consistent 
with gravity- mediated susy breaking^. In the cases presented in these two chapters 
the phase transition is associated with the end of a second (thermal) inflationary 
period, which was originally introduced to dilute the abundances of unwanted relics 
that the flrst (main) inflationary period is not able to do (for instance the moduli 

"'^The parameter e is in these cases unmodified because it is assumed that the unperturbed 
component of a is frozen throughout inflation and until osciUations begin. 



Ill 



fields) |1361 1137j . Thermal inflation is driven by the confinement of a second (flaton) 
field X the origin of the potential due to the thermal effects from the radiation 
background left by the inflaton decay ^Hj- The eventual rolling of Xi towards the 
minimum of its potential, ends thermal inflation and triggers an increment in the 
bare mass of a through the coupling of the latter with x- Solving the moduli 
problem while satisfying adequately all the constraints in the first model discussed 
in Chapter El where the curvaton oscillates for some time before decaying, restricts 
the parameter space to a region where the two important parameters, A and are 
required to be very small (10~^^ ^ ^ % 10~^° and m„ < 10~^ GeV) jSI]. It is likely 
then that the curvaton field is a pseudo Nambu-Goldstone boson. In contrast, as 
discussed in the second model presented in Chapter |31 if the increment in the mass 
of the curvaton field is so high that the decay rate overtakes the Hubble parameter, 
leading to the immediate decay of a, the parameter space is less constrained, with 
more natural values for A and rricr (10~^° ^ ^ 10~^ and mo- < 1 GeV) |189j . 

Although ( is found to be almost perfectly gaussian in most of the models that 
give account of its origin, as required by observation, the possibly present small 
non-gaussianity is being scrutinized by present experiments like the WMAP satel- 
lite |223j and will be the focus of future experiments like the PLANCK satellite 
|171j . The first statistical significant quantity that gives us information about the 
level of non-gaussianity is the bispectrum, which corresponds to the three-point cor- 
relator of C ^103j . Its normalisation is given by the parameter /nl, which has been 
found to be |/nl| < 10^ |lU2j but that will eventually go down to |/nl| ~ 1 in the 
forthcoming years unless there is an earlier detection |lU3j . In view of the difficulty 
at discriminating between models by means of only the spectral index and/or consis- 
tency relations, the detection of non-gaussianity and a precise determination of /nl 
would be useful tools to serve this purpose. For example, in the single-component 
inflaton scenario the level of non-gaussianity is very small, being |/nl| or the order 
of the slow-roll parameters |142j . In contrast, in the curvaton scenario |/nl| could be 
much higher, being /nl negative, or, if a has already dominated the energy density 
before decaying, the precise value for /nl would be /nl = +5/4 fTf l \1'M\ I138j . In 
Chapters El and ini we have addressed this subject, following a perturbative approach 
in the former, and a non-perturbative one in the latter. In Chapter El we have 
presented in an unified way the different definitions of the second-order curvature 
perturbation ^2 present in the literature I13H 11401 1145t 1178] . The translation 
rules to go from one definition to another have been explicitly given, and would help 
to avoid the possible confusion when confronting different papers and results which 
use different definitions for (2 |131j . We have examined the predictions for /nl com- 
ing from the single-component inflation model, the multi-component one, and the 
curvaton scenario, discussing also the respective predictions in preheating and the 
inhomogeneous reheating scenario. Although multi-component inflation, preheating, 
and the inhomogeneous reheating scenarios do not predict a definite value or set of 
values for /nl, being very dependent on the specific model, the single-component 
inflaton scenario as well as the curvaton one do give definite predictions for /nl 
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|131j . Multi-component slow-roll inflation is studied in Chapter IHl following the non- 
perturbative 6N formalism |11UI 11271 1196| 12101 I^Hj - that allows us to calculate ( 
at all orders only by knowing the evolution of a family of unperturbed universes. 
Immediate applications of this formalism are the calculation of the spectral index 
jl96j . and /nl !132j . in a general mult i- component slow- roll inflation model. The 
formalism relies on the separate universe assumption, which says that the local evo- 
lution of the Universe on superhorizon scales is the same as that of an unperturbed 
universe 

Ha inn [UnilinZl 1222] ■ The normalisation parameter /nl is easily given m 
terms of flrst and second derivatives of the number of e-folds, from an initially flat 
slice to a flnal uniform energy density slice, with respect to the flelds living in the 
former. The intrinsic non-gaussianity in the flelds (pi relevant for the evolution of 
the family of universes may be also taken into account, but its contribution, which 
is highly wavevector dependent |142t l2U0j . is in any case negligible compared to 
1 jl41j . Comparison of this formalism with cosmological perturbation theory has 
been made, in the case of the curvaton scenario [132J and the 'hybrid' [HZl model of 
Enqvist and Vaihkonen jl32t I144[ I214j . According to Ref . [144|| the results for the 
latter model following the two approaches agree, and moreover they refute the claim 
about the possible presence of non-local terms in cosmological perturbation theory 

pmiHTiiTra]. 
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Appendix A 

/nl from the 6N formalism 



In this appendix we derive explicitly the expression in Eq. ()fi.2()|l which gives the 
normalisation /nl of the bispectrum of ( defined in Eqs. ()fi.l7j) and (jfi.lSj) . 

We begin by writing ( as given by Eq. ()6.12|) including the homogeneous term 
which makes the spatial average of ( vanish: 



(A.l) 



i ij ij 

The corresponding mode function is then written as 

(A.2) 

Making use of the above formula, the product of three mode functions Ck is 
therefore 



CkiCkaCkg = E^.»^>i^.fc'^ki0i5k20i5k30fc 
ijk 



ijkl 



5ki0i5k30j 



^ki0i(5k2 0j 



(27r)3/2 



(27r)3/2 

5k2+p0fc'5p0r 

3/2'^ki+p0fc5p0i* 



5k3+p0fc^: 



ijkl 



+ 4i0i^k30j^^(k2)((50fc50/) 

+ Sk24>i6k^(f)j6^(ki){6(j)kS(j)i) 



+ 



E N,N,,N 



Im 



ijklm 



d?p d?p' 
(27r)3/2 (27r)3/2 
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5ki</'i42+p<^i'5p(/)^(5k3+p'0z'5, 



P V^m 



+ 



(27r)3/2 '^ki0i(^k2+p0i'^p0fc 

(27r)3/^ '^k3 0i(^k2+p0;'^p0m 

(27r) V2 '^k3 0»'^ki+p0;(^p0^ 
(27r)3/^ '^k20i(^ki+p0j(^p0fc 

+ ^(2^)' E iViiVifeiV;m[5'(ki)5'(k2)((5<^j#fc)(#/#„)(5k3<^, 

+ 53(ki)5^(k3) (50/0,) (50,50„)5k2 0^ 
ijklmn 

I (2^^'^^+p'^'^p'^^* / (^;:^^'^2+p'0feV0; / ^^;:^^k3+p"0m 

ijfeimn 

5^(kl)(50i50j) y" j^^5i^^+p(l)kSp(p^ j ^^^43+p'0m5p'0n 
+ 53(k2)(50fc(50,) y ^^^^4i+p0i(^p0* y" ^^^1^3/2 ^k3+p'0mVC 



^2^^3/2 '^k3+p0m5p<i 
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(27r)3/2 

ij'fc/rrtTi 

(A.3) 

The next step is to take the average of the latter expression. In doing so, we 
make use of the following decompositions jH]: 

(5ki0i5k2</'j)(^k3 0fe'^k4 0«) + (5ki0i5k3 0fc)('^k20i44 0«) + (4i0i5k4 0«)(420i43<^fc) , 

(A.4) 



(4i 01^2 4>j^k3 4>Ai 0/(^k5 0m (^ke 0n) = 

(4i0i5k20i) (5k3 0fc5k4 0i) (5k50m5k6 0n) + {Sk^(t>iSk2<f>j) (^a^fc^s^m) (5k40i5k6 0n) 
+ (4i0i420j)('^k3 0fc46 0n) (44 0^45 0m) + (^ki 0i43 0fc) (42 0j^k4 0i ) ("^kg 0m5k6 0n) 
+ (5ki0j^k3 0fc)(^k20i45 0m)('5k4 0z460n) + ((^ki 0j'5k3 0fc) (^k20j^k6 0n) ((^k4 0i(^k5 0m) 
+ (5ki0i5k4 0i)(^k20j430fc)(5k50m460n) + (^ki 0i5k4 0i ) (^k20j45 0m) (43 0fc5k6 0n) 
+ (4i0j44 0«)(420j^k60n)(5k3 0fc450m) + (^ki 0i45 0m) (420i<^k3 0fc) ((^k4 0i(^k6 0n) 
+ (5ki0j^k5 0m)(^k20i44 0/)('^k30fc46 0n) + ((^ki 0j'5k5 0m) (5k20i'^k6 0n) ('^k3 0fc44 0i) 
+ (5ki0i5k6 0n)(420i5k3 0fc)(44 0«45 0m) + (^ki 0i5k6 0n) (5k2 0i44 0i) (43 0fc'^k5 0m) 
+ (4i0i46 0n) (420i45 0m) (43 0/0^4 0i) > 

(A.5) 

where we have neglected connected n-point correlators with n > 3, which nobody 
has calculated yet but that are presumably very small |2UUj , and products of n-point 
correlators with m-point correlators, where n > 3 and m > 2, which we believe give 
a much smaller contribution than that coming from the three-point correlator |135j . 

From Eqs. ()A.3j) . ()A.4jl . and ()A.5|1 . and after doing some algebra, we obtain the 
three point correlator function of the Ck mode functions: 



(CkiCk2Ck3) = E^.i^.i^.fc(^ki0i5k20i43 0fc) 

ijk 




+ cyclic permutations 



ki + pPp3|k3 - p|3 |ki + p|3p3|k2 - p 



X 
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The integrals in the last line of the previous expression are calculated following the 
arguments presented in Ref. giving as a result 



1 



+ 



1 



|ki + pPp3|k3 _ p|3 |ki + pPp3|]^2 - p 



87rln(A;L) 



cychc 



(A.7) 

where k ^ represents a typical scale under consideration and L is the size of the 
region within which the stochastic properties are specified. That is why \n{kL) can 
be taken of order 1. 

Finally, making use of the /nl definition in Eqs. ()6.17|) and ()6.18|1 : 



(Ck,Ck.Ck3) = (27r)-^/V(ki + k2 + k 



6 

--/nl(A;i, k2, h)[P(^{ki)Pi^{k2) + cyclic] 



(A.8) 

the results found in Eqs. ()A.6|) and ()A.7|1 . and the expression in Eq. ()fi.lfi|l for the 
amplitude of the spectrum of C, in terms of the first derivatives of A^: 



Ai 



we arrive to the desired expression 



UN': 



2 

i 1 



(A.9) 



3. 



" f +ln(fcL) 



ki 



(A.IO) 



with an additional contribution 

T^ijk N^jNjN^kftjkjki, fca, ^3) 



A/] 



NL 



3/2 



(A.U) 



The latter expression comes from the first line in Eq. ()A.6j) containing the three- 
point correlator (5ki0i5k2 0i^k3 0fc)- This three-point correlator is normalised in such 
a way that the bispectra Bij^ of the dimensionless perturbations {2tt / H^)6(f)i is 
written in the same way as in Eqs. ()6.17|) and ()6.18|) with /nl(^i, ^2, ^3) being 
replaced by fijk{ki, ^2, /cs)- 
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